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ABSTRACT 


This manual describes "how to use" the computer code, HITCAN (High Temperature 
Composite ANalyzer) . HITCAN is a general purpose computer program for predicting 
nonlinear global structural and local stress -strain response of arbitrarily 
oriented, multilayered high temperature metal matrix composite structures. This 
code combines composite mechanics and laminate theory with an internal data base 
for material properties of the constituents (matrix, fiber and interphase). The 
thermo -mechanical properties of the constituents are considered to be nonlinear ly 
dependent on several parameters including temperature, stress and stress rate. 
The computation procedure for the analysis of the composite structures uses the 
finite element method. HITCAN is written in FORTRAN 77 computer language and at 
present has been configured and executed on the NASA Lewis Research Center CRAY 
XMP and YMP computers. 

This manual describes HITCAN* s capabilities and limitations followed by input/ 
execution/ output descriptions and example problems. The input is described in 
detail including (1) geometry modeling, (2) types of finite elements, (3) types 
of analysis, (4) material data, (5) types of loading, (6) boundary conditions, 
(7) output control, (8) program options, and (9) data bank. 
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CHAPTER 1 

INTRODUCTION 


The potential use of High Temperature Metal Matrix Composite (HTMMC) materials 
in propulsion systems has already been recognized. The advantages of HTMMC 
materials are high operational temperatures, high specific moduli and strengths, 
tailorable properties, dimensional stability, and hygral resistance. The thermo- 
mechanical properties of components made from HTMMC materials exhibit a nonlinear 
dependence on parameters such as temperature, stress, and stress rate. Since 
comprehensive experimental investigations are prohibitive in cost, it is 
advantageous to have computational schemes which can simulate the nonlinear 
response of components made from HTMMC materials. 

Research related to various aspects of HTMMC materials and structures has been 
conducted at NASA Lewis Research Center (LeRC) for several years. This work has 
focused on high temperature material behavior, constitutive law development, MMC 
experimental mechanics, mathematical modeling, and nonlinear structural analysis 
and simulation. Building upon this research effort, a High Temperature 
Composites ANalyzer (HITCAN) , has been developed. 

HITCAN is a general purpose computer program for predicting global structural and 
local stress-strain response of arbitrarily oriented, multilayered high 
temperature metal matrix composite structures both at the constituent (fiber, 
matrix, and interphase) and the structural level. The thermo -mechanical 
properties of the constituent materials are considered to be nonlinearly 
dependent on several parameters including temperature, stress, and stress rate. 
The computational procedure uses the finite element method, which employs an 
incremental direct iteration procedure to solve the nonlinear equations. A 
schematic of this approach is shown in Figure 1.1. 

HITCAN includes : 

- a dedicated mesh generator, adapted from COBSTRAN (Reference 1); 

- capability for simulating nonlinear behavior at all levels of 
composite material, adapted from METCAN (Reference. 2); 

- finite element structural analysis, adapted from MHOST (Reference 3). 

All three computer programs, COBSTRAN, METCAN, and MHOST were developed in-house 
at NASA and are used as modules. This makes HITCAN a modular stand-alone code, 
independent of commercial codes. 

HITCAN, written in the FORTRAN 77 language, has been configured and executed on 
the LeRC CRAY XMP and YMP computer systems. The code is made up of approximately 
16,000 lines. The companion codes residing in the HITCAN library, COBSTRAN, 
METCAN, and MHOST consist of approximately 7000, 10000, and 51000 lines, 
respectively. 
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Figure 1.1: HITCAN: An Integrated Approach for High Temperature Composite Structural Analysis 
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After a brief description of the capabilities and limitations in the HITCAN 
computer code, a step-by-step outline of the procedure necessary to utilize 
HITCAN, i. e., the preparation of the input data and the creation of a databank 
of constituent material properties and parameters , is given. Chapter 3 describes 
the input data file preparation. The shell scripts which are required to compile 
and run HITCAN on the CRAY X-MP and Y-MP are described in Chapter 4. A 
description of the output is given in Chapter 5. Finally, in Chapter 6, the 
input for three example problems is explained. 

The potential user of HITCAN is reminded that the program is in a state of on- 
going development and the methodology which HITCAN comprises is of an 
evolutionary nature. 

The structure of the input/output file structures in HITCAN is shown in Figure 

1 . 2 . 
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CHAPTER 2 

CAPABILITIES AND LIMITATIONS 


HITCAN is capable of predicting global structural and local stress -strain 
response of multilayered HTMMC structures exhibiting nonlinear material behavior. 
Each layer of the composite can be constructed of different materials and can be 
arbitrarily oriented. The constitutive model employed in METCAN is specifically 
designed for HTMMC, therefore it is recommended that HITCAN be used only for 
metal matrix composites. 

The current version of the code is based on a rectilinear coordinate system. 
Arbitrary shaped geometries can be modeled using interpolators included in the 
mesh generation module of the code. 

At the present time, the following analyses are available in HITCAN: 

- Incremental static analysis with nonlinear anisotropic material behavior 

- Dynamic analysis using direct time integration 

- Modal analysis (free vibration) 

- Buckling analysis (first critical buckling load) 

The element library includes 3 four-noded elements, i.e. plate, plane stress, and 
plane strain, and 1 eight noded element, i.e. a 3D solid. The current mesh 
generation capability of the code allows modeling of solid structures using any 
of the 4 types of elements. The user may also input a finite element model 
directly. The code is capable of handling a variety of boundary conditions, 
loadings (centrifugal, concentrated, distributed, pressure, temperature, static, 
transient, cyclic, and impact), and various types of structures (such as beam, 
plate, ring, curved panel, and built-up structures). A list of HITCAN's analysis 
capabilities can be found in Table 2.1. 

The limitations of the code are : 

- Formulation assumes small displacement and small strain theory; 

- Elements of different types cannot be combined; 

- Hollow structures can be modeled using the plate element only; 

- The finite element model generated by HITCAN can have a nonuniform mesh 
only along the x-axis; 

- If the curvature is large, the mesh will not be uniform. 

HITCAN presents analysis results at the global, laminate, and ply levels. 
Results include displacements , reactions , stresses , and modes shapes . The code 
also has the capability to generate post-processing files for PATRAN. 
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~~ Type of Elements 

Type of Analysis'^ - 

Plate 

Plane 

Strain 

Plane 

Stress 

8-Node 

Solid 

Static 

tested 

- 

- 

tested 

Buckling (a) 

tested 

- 

- 

- 

Load Stepping 

tested 

- 

- 

tested 

Modal (Natural Vibration Modes) (b) 

tested 

- 

- 

- 

Time-domain 

- 

- 

- 

- 

Loading 

Mechanical 

tested 



tested 

Thermal 

tested 



tested 

Cyclic 

- 

- 

- 

- 

Impact 

- 

- 

- 

- 

Constitutive Models (c) 
P = Constant 

tested 


• 

tested 

P = f ( T ) (temperature dependence) 

tested 

! 

i 

tested 

P = f la) (stress dependence) 

tested 

- 


tested 

P = Ha) (stress rate dependence) 

tested 


- 

tested 

P a f(t) (creep) 

- 



- 

P = f( Tcf , a a ) (combination) 

tested 

_ 


tested 

P = f( Ta, a .a , t ) (creep combination) 

- 

_ 

_ 

- 

Fiber Degradation 

tested 


_ 

tested 

Fabrication-induced Stresses 

tested 



tested 

Ply Orientations 
Arbitrary 

tested 

- 

- 

tested 


(a) Tested 1 buckling mode (c) Constitutive models: Notation 

(b) Tested 4 vibration modes E : P. 1 ?? 61 *' 65 ?'■ Stress 

T: Temperature q\ Stressrate 

t : Time 


Table 2.1: HITCAN Capabilities for Composite Materials 
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CHAPTER 3 


INPUT FORMAT 

A single HITCAN input file includes all the data necessary for the selection of 
analysis options, parameters, mesh generation, composite material type and 
construction, loading, initial conditions, boundary conditions, and print 
options. However, the values of material properties at the constituent (fiber, 
matrix, and interphase) level are not entered in the HITCAN input file, they 
reside in a separate file, labeled "data bank”. This file contains material 
property data for each of the composite systems used in the analysis. Each 
composite system is identified with a material name which is entered in the 
HITCAN input file. The user can edit the "data bank" file to define material 
property data for any material required. Additional information on the "data 
bank" file can be found in Section 3.8. 

The input file consists of two blocks. The first block contains the title and 
the program option cards. The second block consists of card groups. Figure 3.1 
illustrates the two blocks. The program option cards either control the flow 
through the program or activate various card groups. There are twenty-eight 
program option cards, four of these control the flow through the program. They 
are HPLATE, S3DS0LID, and SPLATE option cards which set the type of finite 
element model to be generated by HITCAN, and the READ IN MODEL option which 
enables the user to input into HITCAN a finite element model consisting of eight 
node solid elements. One of these four cards must be included in the "Program 
Option Cards" block of the input deck. If the program option card HPLATE is 
specified then the following program option cards must be included in the input 
file 

PLATE 

PLYORDER 

ENDOPTION 

If SPLATE is used, the program option cards required in the input file are 

PLATE or STRESS or STRAIN 

PLYORDER 

ENDOPTION 

When the S3DS0LID is specified, then 

BRICK 

PLYORDER 

ENDOPTION 


must be used. 
Chapter 3 
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To separate the two blocks the ENDOPTION card is used. This program option card 
is placed at the end of the "Program Option Cards", as shown in Figure 3.1. 
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Figure 3.1: Organization of Input File in HITCAN 
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I 


A card group is a set of one or more cards. Some of the card groups are 
dependent on either a program option card being specified or a parameter defined 
in a previous card group. These optional card groups and their respective 
dependencies are explained within each card group in Section 3.3. If an optional 
card group is not activated by a previous program option card or user selected 
parameter, the optional card group must not be included in the input file. The 
input file should , however . be maintained in numerical order by card group 
number . 

The input file can be broken into 7 functions. These are shown in Figure 3.2. 
In sections 3.1 to 3.7, the program option cards and card groups used in each one 
of these functions will be described. By depicting the input file in this 
manner, the user can quickly assemble an input file. In section 3.9 a summary 
of all the program option cards is given. 

Note that in both the input file and the "data bank" file English units are 
required. 
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3.1 GEOMETRIC MODELING 

HITCAN generates a finite element mesh of the structure, based on coordinates of 
a few representative points . This is accomplished by interpolating nodal values 
of geometrical coordinates using a cubic spline. The interpolating function 
requires that the surface geometry be represented by a continuous and s ingle - 
valued function. Although HITCAN is a stand-alone computer program, it does 
create a PATRAN neutral file, so that the user can view the mesh generated. 
Figure 3.3 shows the different types of structures that can be modeled by HITCAN. 
There are 3 distinct mesh generation schemes in HITCAN, as categorized below. 
Any 1 of the 3 mesh generation schemes is activated by choosing the appropriate 
program option card, shown in the parentheses below. 

(1) For solid structures: use plate, plane stress, or plane strain element 
(SPLATE) . 

(2) For hollow structures (panel structure): use plate element (HPLATE) . 

(3) For solid structures: use 3D solid element (S3DS0LID) . 

The coordinate system must be rectilinear and right-handed. The center of the 
coordinate system can be placed anywhere. However, in choosing the directions 
of the three axes, the user should keep it in mind that the mesh can be 
nonuniform only along the x-axis. 

Also, HITCAN provides the user with the ability to enter a finite element model 
created by a program other than HITCAN. At the present time, this option can 
only be used with the eight-node solid element. This option can be activated 
with the program option card READ IN MODEL. 

The following table summarizes the different mesh generation schemes and their 
program option cards and card groups. 


TYPE OF STRUCTURE 

PROGRAM OPTION CARD 

CARD GROUPS 

Solid structure using 



plane stress, plane strain. 

SPLATE 

5 

or plate element 


27 

Hollow structures 

HPLATE 

2 



9 



24 

Solid structures 



using eight-node 

S3DSOLID 

4 

solid element 


26 

Read in a finite 


3 

element model 

READ IN MODEL 

w 

25 
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CYLINDRICAL SHELL ROOF 
use: SPLATE Element 


BLADE 

use: SPLATE Element 



ONE QUARTER OF A RING HOLLOW STRUCTURE 

use: SPLATE Element (PANEL STRUCTURE) 

use: HPLATE Element 


Figure 3.3: Typical Structures That Can Be Modeled Using HITCAN 
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S PLATE MODEL OPTION 

The S PLATE model option allows modeling of structures with arbitrary shapes in 
x and y directions, but a through- the* thickness plate type shape in the z- 
direction. An example of the type of structure that can be model with HITCAN is 
a solid curved panel shown in Figure 3.4. The finite element mesh is generated 
in the following 4 steps. 

Step 1: A right-handed rectilinear coordinate system (x,y,z) is 

defined, placing the center of the coordinate system at a 
convenient point. 

Step 2: To obtain the surface geometry of a structure, it must be 

divided into several (y,z) sections along the x-axis. The 
number of cross sections selected along the x-axis depends upon 
the curvature of the structure along the x-axis. The mesh 
generator can fit any curve up to a third degree polymonial, 
using a cubic spline. Hence, the x-axis must be divided 
into enough sections such that each section can be modeled 
using a cubic spline. 

Step 3: Each (y,z) cross section along the x-axis consists of a set of 

points. The number of (y,z) points needed to define a specific 
cross section are selected so that the curve between the two 
adjacent points can be modeled using a cubic spline 
interpolator. A different number of (y,z) points can be 
selected for different cross sections, along the x-axis. 

Step 4: Once the structure is divided into cross sections , along the x- 

axis in steps 2 and 3 above, the user has already created a 
coarse finite element mesh. To further subdivide the x and 
(y,z) sections for obtaining the desired number of elements, 
the user needs only to input the number of intermediate nodes. 
The nodal coordinates for the intermediate nodes of the finite 
element mesh are automatically interpolated. 
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Figure 3.4: A Shell Structure Modeled Using The SPLATE Option in HITCAN 
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Note that (y,z) coordinates in each cross section can be input in 2 different 
ways , as follows : 

(1) For each y point, two values of the z-coordinate , for the upper and 
the lower surface, must be input (see Figure 3.5). The thickness and 
the mid- surface of the structure are automatically calculated by the 
code . This option is activated by entering a value of 1 or 2 for the 
variable IGRD. The geometry modeling data are input through 
variables X, Y, ZU, and ZL for each cross section. 

(2) For each y point, one value of the z-coordinate, for the mid-plane 
of the structure and the structure thickness must be input (see 
figure 3.5). The upper and the lower surfaces of the structure are 
automatically calculated by the code. This option is activated by 
inputting a value of 3 or 4 for the variable IGRD. The geometry 
modeling data are input through the variables X, Y, Z, and TB for 
each cross section. 
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Figure 3.5: Input Options For The SPLATE Model Option 
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Note that by setting IGRD equal to 1 or 4, HITCAN will divide the model into IU-1 
elements along the x-axis and JU-1 elements along the y-axis using a cubic spline 
interpolating function. With IGRD having a value of 2 or 3, the input points are 
assumed to be the nodes and the variables X, Y, ZL, ZU, Z, and TB are nodal 
quantities. HITCAN will automatically generate the element connectivities. 
These variables are illustrated in the Figure 3.6 on the following page. 

Program Option card: SPLAT E MODEL 

Card Groups: 5 and 27 

Variable 

Columns Format Name Entry 

1st card in card group 5 

1-4 14 NSECT The number of input cross sections. A 

sufficient number must be used so that the 
surface geometry can be properly represented 
by a cubic spline. Note, if the surface is 
linear only two points are needed. The 
coordinates of the points are entered below, 
using X, Y, Z, ZL, ZU, and TB. 

2nd card in card group 5 

1-4 14 IGRD Sets the input format. If IGRD is set to 1 

or 4, a finite element mesh will be 
generated. Setting IGRD equal to 2 or 3 the 
input points will be assumed to be nodes of 
a finite element mesh, the element 
connectivity will then be automatically 
created. 

5-8 14 IU For IGRD equal to 1 or 4, the number of nodes 

in the finite element model along the x-axis. 
With IGRD set to 2 or 3, IU is the number of 
nodes that are input along the x-axis; i. e. , 
IU is equal to NSECT. 

9-12 14 JU For IGRD equal to 1 or 4, the number of nodes 

in the finite element model along the y-axis . 
With IGRD set to 2 or 3, JU is the number of 
nodes that are input along the y-axis; i. e. , 
JU is equal to MSECT. 
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Variable 

Columns Format Name Entry 

3rd card in card group 5 

1-8 F8.4 XBEG The initial x coord, of structure, use if 

IGRD - 1 or 4 

9-16 F8.4 XEND The final x coord, of structure, use if IGRD 

- 1 or 4 

1st card in card group 27 

1-80 2014 MSECT Number of input points at each cross section. 

If IGRD is equal to 2 or 3, then the number 
of input points must be the same for each 
cross section. 


2nd card in card group 27 

1-4 14 NXDIV Number of elements between two consecutive 

output sections 


3rd card 

in card group 

27 






1-8 

F8.4 

X 

X coordinate 

of 

an 

input 

point 

9-16 

F8.4 

Y 

Y coordinate 

of 

an 

input 

point 

17-24 

F8.4 

Z or ZU 

If IGRD - 1 

or 

2, 

use 

ZU. ZU is the Z 


coordinate of the upper surface. If IGRD - 
3 or 4, use Z. Z is the Z coordinate of the 
mid-plane. 


25-32 F8.4 TB or ZL If IGRD - 1 or 2, use ZL. ZL is the Z 

coordinate of the lower surface. If IGRD - 
3 or 4, use TB. Tb is the wall thickness. 

One block of data for each input cross section. Each block will contain MSECT(J) 
cards of card #3 and 1 card of card #2. The coordinates of each input point will 
be on one card. The total number of cards will be MSECT(l) + MSECT(2) +. . .+ 
MSECT(I) + NSECT, where I-l.NSECT. 
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EXAMPLE: 


In this example, the input required for a plate of length 4 in. , of width 2 in. , 
and a thickness of 0.1 in. is given. In the finite element mesh there will be 
4 elements along the x-axis and 4 elements along the y-axis. 


Card GrouD 5 





1 

2 

3 

4 5 6 

1 

O 

. .0 

. . .0 

0 

0 0 0 

Z 

4 

4 4 




0. 

0 4.1 

0 



Card Group 27 





1 

2 

3 

4 5 6 

1 

. . .0 

0 

0 

0 0 0 

2 

3 

2 





0.0 

-1.0 

0.0 

0.1 


0.0 

1.0 

0.0 

0.1 

0 

4.0 

-1.0 

0.0 

0.1 


4.0 

1.0 

0.0 

0.1 

Columns 

; Field 

Name 

Value 

Description 

1-4 

NSECT 


i 

There is 2 input sections. 

1-4 

IGRD 


4 

Since IGRD equals - 4, a finite element 
model will be created. 

5-8 

IU 


4 

Number of nodes along the x-axis. 

9-12 

JU 


4 

Number of nodes along the y-axis. 

1-8 

XBEG 


0.0 

The initial x coordinate of the plate. 

9-16 

XEND 


4.0 

The final x coordinate of the plate. 

1-4 

NXDIV 


3 

Number of elements between the 1st and 
2nd output sections. 

1-10 

X 


0.0 

X coordinate of input point #1 of input 
section #1. 

11-20 

Y 


-1.0 

Y coordinate of input point #1 of input 
section #1. 

21-30 

Z 


0.0 

Z coordinate of input point #1 of input 
section #1. 

31-40 

TB 


0.1 

The thickness of input point #1 of 
input section #1. 

1-10 

X 


0.0 

X coordinate of input point #2 of input 
section #1. 

11-20 

Y 


1.0 

Y coordinate of input point #2 of input 
section #1. 
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Columns 

Field Name 

Value 

21-30 

z 

0.0 

31-40 

TB 

0.1 

1-4 

NXDIV 

0 

1-10 

X 

4.0 

11-20 

Y 

-1.0 

21-30 

Z 

0.0 

31-40 

TB 

0.1 

1-10 

X 

4.0 

11-20 

Y 

1.0 

21-30 

Z 

0.0 

31-40 

TB 

0.1 


Description 

Z coordinate of input point #2 of input 
section #1. 

The thickness of input point #2 of 
input section #1. 

X coordinate of input point #1 of input 
section #2. 

Y coordinate of input point #1 of input 
section #2. 

Z coordinate of input point #1 of input 
section #2. 

The thickness of input point #1 of 
input section #2. 

X coordinate of input point #2 of input 
section #2. 

Y coordinate of input point #2 of input 
section # 2 . 

Z coordinate of input point #2 
of input section #2. 

The thickness of input point #2 of 
input section #2. 
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HPLATE MODEL OPTION 

The HPLATE model option allows modeling of hollow sandwich type structures with 
different shaped top and bottom plates, joined by multiple different sized plate 
type spars. Figure 3.7 shows an example of this type of structure. The top and 
bottom plates can have arbitrary shapes in x and y directions . The spars are 
defined as plates with faces in the x,z plane and thickness in the y-direction, 
joining the top and bottom plates, at user-specified points. Note if the spars 
in the structure are evenly spaced, then the program option card PANEL must also 
be specified. 

HITCAN assumes the nodal input points are on the outer surface of the shell. 
HITCAN automatically corrects the grid point positions by moving them to the mid- 
wall position in a direction normal to the surface. The program option card 
will suppress this mid-wall correction and will retain the grid points on the 
external profile of the shell. 

The finite element mesh is generated in the following steps . 

Step 1: A right-handed rectilinear coordinate system (x,y,z) is 

defined, placing the center of the coordinate system at a 
convenient point. 

Step 2: The basic shape of the top and bottom plates is divided into 

cross sections ((y,z) planes), definable with cubic splines. 
The number of cross sections selected along the x-axis depends 
upon the curvature of the structure along the x-axis. The 
curvature along the x-axis can be different for the top and 
bottom plates of the sandwich structure. This must 
be kept in mind when selecting the number of cross sections. 

Step 3: Each cross section ((y,z) plane) consist of input points. The 

number of (y,z) input points needed to define a specific cross 
section are again selected so that the curve between the two 
adjacent points can be modeled using a cubic spline 
interpolator. A different number of (y,z) input points can be 
selected for different cross sections and for the top and 
bottom surfaces. 

Step 4: The desired nodal point coordinates are automatically 

interpolated, by defining the number of elements desired 
between each output section and the number of elements desired 
along the y-axis. The spars are divided into same number of 
elements along the x-axis as the top and the bottom plates. 
The spars can only consist of one element in the y and z 
directions . 

Figure 3.8 illustrates several of the variables used in the HPLATE model option. 
Program Option card: 
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Figure 3.7: A Panel Structure Modeled Using The HP LATE Option in HITCAN 
Chapter 3 March, 1992 
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Figure 3,8: Variables Used In The HPLATE Model Option 
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( 


Program Option Cards: HPLAT E MODEL OPTION PROFILE 

PROFILE 

PANEL 


Card Groups: 2, 9, and 

Variable 


Columns 

Format 

Name 

1st card 

in card 

group 

2 

1-4 

14 


NSECT 

5-8 

14 


NXSPAR 

9-12 

14 


NSPAR 

13-16 

14 


NYSPC 

1st card 

in card 

group 

9 

1-80 

2014 


LSECT(l) 

1-80 

2014 


LSECT(l) 

1st card 

of card 

group 

24 

1-80 

2014 


NSPDES 

2nd card 

of card 

group 

24 

1-4 

14 


NXSPC 

5-10 

F6.2 


XH 


I 


24 


Entry 


Number of output sections . 

Number of output sections containing spars 
Number of spars in an output section. Note 
that each output section must have the same 
number of spars. 

Number of elements desired along the y-axis. 
The top and bottom plates will have the same 
number of elements . If PANEL is specified 
the number of elements will be SPC+1. 


Number of input cross sections for the top 
surface. A sufficient number must be used 
so that the surface geometry can be properly 
represented by a cubic spline. 

Number of input cross sections. 


Spar ply designation numbers. 


Number of elements in the FE model 
between output sections 
X coordinate of an output section 


I 
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3rd card in card group 24 


1-10 

F10.4 

SY(1) 

11-20 

F10.4 

SY(2) 

21-30 

F10.4 

SY(1) 

31-40 

F10.4 

SY(2) 

41-50 

F10.4 

SY(1) 

51-60 

F10.4 

SY(2) 

61-70 

F10.4 

SY(1) 

71-80 

F10.4 

SY(2) 

Cards 2 and 

3 are 

repeated NSECT+1 

4th card of 

card group 24 

1-80 

2014 

MSECT 


5th card 

of card group 

24 

1-80 

2014 

MSECT 

6th card 

in card group 

24 

1-8 

F8.4 

X 

9-16 

F8.4 

Y 

16-24 

F8.4 

Z 

25-32 

F8.4 

THK 


Y coordinate of spar as measured from the 
left edge of the structure 

Wall thickness of spar 

Y coordinate of spar as measured from the 
left edge of the structure 

Wall thickness of spar 

Y coordinate of spar as measured from the 
left edge of the structure 

Wall thickness of spar 

Y coordinate of spar as measured from the 
left edge of the structure 

Wall thickness of spar 

times 


An array containing the number of input 
points for each input cross section. This 
is for the top surface. 


An array containing the number of input 
points for each input cross section. This 
is for the top surface. 


X coordinate of an input point on the top 
surface 

Y coordinate of an input point on the top 
surface 

Z coordinate of an input point on the top 
surface 

Wall thickness of an input point on the top 
surface 
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This card is repeated MSECT(l.l) + MSECT(2,1), +. . .+ MSECT(I,1), 
where I - 1,LSECT(1). 


7th card in card group 24 


1-8 

F8.4 

X 

X coordinate 
surface 

of an 

input point 

on 

the 

top 

9-16 

F8.4 

Y 

Y coordinate 
surface 

of an 

input point 

on 

the 

top 

16-24 

F8.4 

Z 

Z coordinate 
surface 

of an 

input point 

on 

the 

top 

25-32 

F8.4 

THK 

Wall thickness of an input point 
surface 

on 

the 

top 


This card is repeated MSECT(1,2) + MSECT(2,2), +. . .+ MSECT(I,2), 
where I - 1,LSECT(2). 
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EXAMPLE: 


In this example, the input is given for a panel. The length of the panel is 0.5 
in., the width is 0.2 in., and the height of the panel is 0.075 in. The panel 
has 2 spars one on each end. Each spar has a thickness of 0.02 in. 

Card Group 2 

1 2 3 4 5 6 


1 0 0 0 0 0 0 


2 2 

Card Group 

2 7 

9 




1 

2 


3 

4 5 6 

1 0, 

2 2 

Card Group 

0. . . 

24 


. . .0. . . 

0 0 0 

1 

2 


3 

4 5 6 

1 0, 

1 1 
4 .0 

0. . . 

* * * 

. . .0. . . 

0 0 0 

.0 

.5 

.02 


.2 

.02 

.0 

2 2 

2 2 

.02 


.2 

.02 

.0 

-.1 


.04 

.02 

.0 

.1 


.04 

.02 

.5 

-.1 


.04 

.02 

.5 

.1 


.04 

.02 

.0 

- .1 


.035 

.01 

.0 

.1 


.035 

.01 

.5 

-.1 


.035 

.01 

.5 

.1 


.035 

.01 

Columns 

Field Name 

Value 

Description 

1-4 

NSECT 

i 


Number of output sections is 1. 

5-8 

NXSPAR 

2 


Both output sections are to contain spars . 

9-12 

NSPAR 

2 


Number of spars is 2 . 

13-16 

NYSPC 

7 


Number of elements along the Y-axis 

1-4 

LSECT(l) 

2 


The top surface is described by 2 input 
sections. 

5-8 

LSECT(2) 

2 


The bottom surface is described by 2 input 
sections . 

1-4 

Chapter 3 
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Columns 

Field Name 

Value 

Description 

5-8 

NSPDES(2) 

1 

Spar #2 is described by ply designation #2. 

1-4 

NXSPC 

4 

Number of elements along the x-axis between 
output sections #1 and #2. 

5-10 

XH 

0.0 

The x coordinate of the 1st output section. 

1-10 

SY(1) 

0.0 

Y coordinate of spar #1 . 

11-20 

SY(2) 

0.02 

Wall thickness of spar #1. 

21-30 

SY(1) 

0.2 

Y coordinate of spar #2. 

31-40 

SY(2) 

0.02 

Wall thickness of spar #2. 

1-4 

NXSPC 

0 


5-10 

XH 

0.0 

The x coordinate of the 2nd output section. 

1-10 

SY(1) 

0.0 

Y coordinate of spar #1. 

11-20 

SY(2) 

0.02 

Wall thickness of spar #1. 

21-30 

SY(1) 

0.2 

Y coordinate of spar #2. 

31-40 

SY(2) 

0.02 

Wall thickness of spar #2. 

1-4 

MSECT(l) 

2 

Two inputs described the 1st input section 
on the top surface. 

5-8 

MSECT(2) 

2 

Two inputs described the 1st input section 
on the top surface. 

1-4 

MSECT(l) 

2 

Two inputs described the 1st input section 
on the bottom surface. 

5-8 

MSECT(2) 

2 

Two inputs described the 1st input section 
on the bottom surface. 

1-10 

X 

0.0 

X coordinate of input point #1 at input 
section #1 on the top surface. 

11-20 

Y 

-.1 

Y coordinate of input point #1 at input 
section #1 on the top surface. 

21-30 

Z 

0.04 

Z coordinate of input point #1 at input 
section #1 on the top surface. 

31-40 

THK 

0.02 

Wall thickness of input point #1 at input 
section #1 on the top surface. 

1-10 

X 

0.5 

X coordinate of input point #2 at input 
section #1 on the top surface. 

11-20 

Y 

0.1 

Y coordinate of input point #2 at input 
section #1 on the top surface. 

21-30 

Z 

0.04 

Z coordinate of input point #2 at input 
section #1 on the top surface. 

31-40 

THK 

0.02 

Wall thickness of input point #1 at input 
section #2 on the top surface. 

1-10 

X 

0.5 

X coordinate of input point #1 at input 
section #2 on the top surface. 

11-20 

Y 

-.1 

Y coordinate of input point #1 at input 
section #2 on the top surface. 

21-30 

Z 

0.04 

Z coordinate of input point #1 at input 
section #2 on the top surface. 

31-40 

THK 

0.02 

Wall thickness of input point #1 at input 
section #2 on the top surface. 
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Columns 

1-10 

Field Name 
X 

Value 

0.0 

Description 

X coordinate of input point #2 
at input section #2 on the top 
surface . 



11-20 

Y 

0.1 

Y coordinate of input point #2 
section #2 on the top surface. 

at 

input 

21-30 

Z 

0.04 

Z coordinate of input point #2 at input 
section #2 on the top surface. 

31-40 

THK 

0.02 

Wall thickness of input point #2 
section 42 on the top surface. 

at 

input 

1-10 

X 

0.0 

x coordinate of input point #1 at 
section #1 on the bottom surface. 

input 

11-20 

Y 

-.1 

Y coordinate of input point 4 1 at input 
section #1 on the bottom surface. 

21-30 

Z 

0.035 

Z coordinate of input point #1 
section #1 on the bottom surface. 

at 

input 

31-40 

THK 

0.02 

Wall thickness of input point #1 
section #1 on the bottom surface. 

at 

input 

1-10 

X 

0.0 

X coordinate of input point 42 
section #1 on the bottom surface. 

at 

input 

11-20 

Y 

0.1 

Y coordinate of input point 42 
section #1 on the bottom surface. 

at 

input 

21-30 

Z 

0.035 

Z coordinate of input point 42 
section #1 on the bottom surface. 

at 

input 

31-40 

THK 

0.01 

Wall thickness of input point #1 
section #2 on the bottom surface. 

at 

input 

1-10 

X 

0.5 

X coordinate of inpt point #1 
section 42 on the bottom surface. 

at 

input 

11-20 

Y 

-.1 

Y coordinate of input point #1 
section 42 on the bottom surface. 

at 

input 

21-30 

Z 

0.035 

Z coordinate of input point #1 
section #2 on the bottom surface. 

at 

input 

31-40 

THK 

0.01 

Wall thickness of input point #1 
secttion 42 on the bottom surface 

at 

input 

1-10 

X 

0.5 

X coordinate of input point 42 
section 42 on the bottom surface. 

at 

input 

11-20 

Y 

0.1 

Y coordinate of input point 42 
section 42 on the bottom surface. 

at 

input 

21-30 

Z 

0.035 

Z coordinate of input point 42 
section 42 on the bottom surface. 

at 

input 

31-40 

THK 

0.01 

Wall thickness of input point 42 
section 42 on the bottom surface. 

at 

input 
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i 


S3DSOLID MODEL OPTION 

This model option enables the user to create a finite element model of a solid 
structure using 3D solid elements. Figure 3.9 shows an example of the type of 
structure that can be modeled using this model option. The finite element mesh 
is generated in the following steps. 

Step 1: A right-handed rectilinear coordinate system (x,y,z) is 

defined, placing the center of the coordinate system at a 
convenient point on the structure. 

Step 2: The basic shape of the structure is divided into input planes. 

As shown in the Figure 3.10, a quarter of a ring that was 
broken into two input planes. One input plane represents the 
inside surface of the ring, the other the outside surface of 
the ring. The number of input planes used will depend upon the 
curvature of the structure along the z-axis. A sufficient 
number must be provided so that a cubic spline can be used. 

Step 3: Each input plane is separately divided into several sets of 

input points. The number of sets used must be sufficient so 
that the curvature in the input plane can be represented by a 
cubic spline. The number of sets of points for each input plane 
can be different; likewise, the number of points in each set 
can be different. 

Step 4: The desired nodal point coordinates are automatically 

interpolated, by defining the number of elements desired 
between each output section, the number of elements along the 
y-axis, and the number of elements along the z-axis. 
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Figure 3.9: A Ring Modeled Using The S3DSOLID Option In HITCAN 
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Program Option card: S3DS OLID MODEL OPTION 

Card Groups: 4, and 26 


Columns 

Format 

Variable 

Name 

1st card of 

card group 4 


1-4 

14 

NIPL 


5-8 

14 

NOSC 

9-12 

14 

NEYY 

13-16 

14 

NETT 


Entry 

Number of input planes . A sufficient 
number must be used so that the 
curvature can be properly represented 
by a cubic spline. 

Number of output sections . 

Number of elements along the 
y-axis in the finite element model. 
Number of elements along the z-axis in 
the finite element model. 


1st card of card group 26 

1-4 14 NXSPC Number of elements between output 

sections along the x-axis. 

5-10 F6.2 X X coordinate of an output section. 

The above card is repeated NOSC+1 
times . 


2nd card of card group 26 

1-80 2014 LSECT 


3rd card of card group 26 

1-80 2014 MSECT 


Card #3 is repeated NIPL times . 

4th card of card group 26 

1-8 F8.4 X 

9-16 F8.4 Y 

17-24 F8.4 Z 


An array containing the number of sets 
of points in an input plane. A 
sufficient number must be used so that 
the curvature lying in the input plane 
can be properly represented. 


An array containing the number of input 
points for each set of input points . 


X coordinate of 
Y coordinate of 
Z coordinate of 


an input point, 
an input point, 
an input point. 


This card is repeated as follows. One block of data for each input plane, i. e. , 
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NIPL blocks of data. Each block will contain LSECT sets of points, with each set 
containing MSECT input points. Thus for input plane #2, card #4 is repeated 
MSECT(l) + MSECT(2) +. . .+ MSECT(I) times, where I - 1,LSECT(2). 


I 
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EXAMPLE: 

The input for a sector of a ring is shown in this example. The inside radius of 
the ring is 2.875 in. and the outside radius of the ring is 3.475. 

Card Group 4 


1 2 3 4 5 6 

1 0 0 0 0 0 0 

2 3 3 6 

Card Group 26 

1 2 3 4 5 6 

1 0 0 0 0 0 0 

2 .0 

.6 


2 

2 




2 

2 




2 

2 





.0 

.0 

2.875 



.0 

.1 

2.873 



.6 

.0 

2.875 



.6 

.1 

2.873 



.0 

.0 

3.475 



.0 

.121 

3.473 



.6 

.0 

3.475 



.6 

.121 

3.473 


Columns 

Field Name 

Value 

Description 

1-4 


NIPL 

2 

Number of input planes. 

5-8 


NOSC 

3 

Number of output sections. 

9-12 


NEYY 

3 

Number of elements along the Y-axis. 

13-16 


NETT 

6 

Number of elements through the thickness. 

1-4 


NXSPC 

2 

Number of elements between output sections 
#1 and #2. 

5-10 


X 

0.0 

X coordinate of the 1st output section. 

1-4 


NXSPC 

0 


5-10 


X 

0.6 

X coordinate of the 2nd output section. 

1-4 


LSECT(l) 

2 

Number of input sections on input plane #1. 

5-8 


LSECT(2) 

2 

Number of input sections on input plane #2 . 

1-4 


MSECT(l) 

2 

Number of input points in the 1st set of 
input points of input plane #1. 

5-8 


MSECT(2) 

2 

Number of input points in the 2nd set of 


input points of input plane #1. 
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I 


Columns 

1-4 

Field Name 
MSECT(l) 

Value 

2 

5-8 

MSECT(2) 

2 

1-8 

X 

0.0 

9-16 

Y 

0.0 

17-24 

Z 

2.875 

1-8 

X 

0.0 

9-16 

Y 

0.1 

17-24 

Z 

2.873 

1-8 

X 

0.6 

9-16 

Y 

0.0 

17-24 

Z 

2.875 

1-8 

X 

0.6 

9-16 . 

Y 

0.1 

17-24 

Z 

2.873 

1-8 

X 

0.0 

9-16 

Y 

0.0 

17-24 

Z 

3.475 

1-8 

X 

0.0 

9-16 

Y 

0.121 

17-24 

Z 

3.473 

1-8 

X 

0.6 

9-16 

Y 

0.0 

17-24 

Z 

3.475 


Chapter 3 


Description 
Number of input 
input points of 
Number of input 
input points of 
X coordinate of 
input points on 

Y coordinate of 
input points on 
Z coordinate of 
input points on 
X coordinate of 
input points on 

Y coordinate of 
input points on 
Z coordinate of 
input points on 
X coordinate of 
input points on 

Y coordinate of 
input points on 
Z coordinate of 
input points on 
X coordinate of 
input points on 

Y coordinate of 
input points on 
Z coordinate of 
input points on 
X coordinate of 
input points on 

Y coordinate of 
input points on 
Z coordinate of 
input points on 
X coordinate of 
input points on 

Y coordinate of 
input points on 
Z coordinate of 
input points on 
X coordinate of 
input points on 

Y coordinate of 
input points on 
Z coordinate of 
input points on 


points in the 1st set of 
input plane # 2 . 
points in the 2nd set of 
input plane #2. 
point # 1 in the 1st set of 
plane #1. 

point #1 in the 1st set of 
plane #1. 

point # 1 in the 1st set of 
plane #1. 

point #2 in the 1st set of 
plane #1. 

point #2 in the 1st set of 
plane #1. 

point #2 in the 1st set of 
plane #1. 

point #1 in the 2nd set of 
plane #1. 

point #1 in the 2nd set of 
plane #1. 

point #1 in the 2nd set of 
plane #1. 

point #2 in the 2nd set of 
plane #1. 

point #2 in the 2nd set of 
plane #1. 

point #2 in the 2nd set of 
plane #1. 

point #1 in the 1st set of 
plane #2. 

point #1 in the 1st set of 
plane # 2 . 

point #1 in the 1st set of 
plane #2. 

point #2 in the 1st set of 
plane #2. 

point #2 in the 1st set of 
plane #2. 

point #2 in the 1st set of 
plane #2. 

point //I in the 2nd set of 
plane #2. 

point #1 in the 2nd set of 
plane #2. 

point #1 in the 2nd set of 
plane #2. 
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Columns 

Field Name 

Value 

1-8 

X 

0.6 

9-16 

Y 

0.121 

17-24 

2 

3.473 


Description 

X coordinate of point #2 in the 2nd set of 
input points on plane #2. 

Y coordinate of point #2 in the 2nd set of 
input points on plane #2. 

2 coordinate of point #2 in the 2nd set of 
input points on plane #2. 
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READ IN MODEL 

This option allows the user to input a finite element model using eight-node 
solid elements. The user must provide the element connectivities and the nodal 
coordinates. When this option is chosen, the user must also specify the BRICK 
program option card. 


Program option cards : READ IN MODEL and BRICK 


Card groups: 


3 and 25 


Columns Format 


Variable 

Name Entry 


1st card 

in card 

group #3 

1-4 

14 

MAXNP 

5-8 

14 

NETOT 


Number of nodes in the finite element model. 
Number of elements in the finite element 
model. 


1st card in card 

group #25 

1-4 

14 

IELE 

5-36 

814 

KEI 


Element number. 

Element connectivity. See description of the 
solid element for the correct order of the 
nodes in the element connectivity. 


Card #1 is repeated NETOT times. 
2nd card in card group #25 


1-5 

15 

NNUM 

Node number. 

6-10 

11-20 

5X 

F10.4 

X 

X coordinate. 

21-30 

F10.4 

Y 

Y coordinate. 

31-40 

F10.4 

Z 

Z coordinate. 

Card #2 

is repeated 

MAXNP times 



I 


l 
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EXAMPLE: 

Card Group 3 

1 2 3 4 5 6 

1 0 0 0 0 0 0 

8 1 


Card Grout) 25 

1 2 
i n n 

3 

0 

4 5 

0 0 . . 

6 

0 

1 1 


2 3 4 

5 6 7 

8 


1 


1.0 

1.0 

0.0 


2 


0.0 

1.0 

0.0 


3 


0.0 

1.0 

1.0 


4 


1.0 

1.0 

1.0 


5 


1.0 

0.0 

0.0 


6 


0.0 

0.0 

0.0 


7 


0.0 

0.0 

1.0 


8 


1.0 

0.0 

1.0 


Columns 


Field Name 

Value 

Description 


1-4 


MAXNP 

8 

Number of nodes 

is 8. 

5-8 


NETOT 

1 

Number of elements is 1. 

1-4 


IELE 

1 

Designation number of element #1. 

5-8 


KEI(l) 

1 

1st node in the 

element connectivity. 

9-12 


KEI(l) 

2 

2nd node in the 

element connectivity. 

13-16 


KEI(l) 

3 

3rd node in the 

element connectivity. 

17-20 


KEI(l) 

4 

4th node in the 

element connectivity. 

21-24 


KEI(l) 

5 

5th node in the 

element connectivity. 

25-28 


KEI(l) 

6 

6th node in the 

element connectivity. 

29-32 


KEI(l) 

7 

7th node in the 

element connectivity. 

32-36 


KEI(l) 

8 

8th node in the 

element connectivity. 

1-5 


NNUM 

1 

Node number for 

node #1. 

11-20 


X 

1.0 

X coordinate of 

node #1. 

21-30 


Y 

1.0 

Y coordinate of 

node #1. 

31-40 


Z 

0.0 

Y coordinate of 

node #1. 

1-5 


NNUM 

2 

Node number for 

node #2. 

11-20 


X 

0.0 

X coordinate of 

node it 2. 

21-30 


Y 

1.0 

Y coordinate of 

node it 2. 

31-40 


Z 

0.0 

Y coordinate of 

node #2. 

1-5 


NNUM 

3 

Node number for 

node #3. 

11-20 


X 

0.0 

X coordinate of 

node #3. 

21-30 


Y 

1.0 

Y coordinate of 

node #3. 

31-40 


Z 

1.0 

Y coordinate of 

node #3 . 

1-5 


NNUM 

4 

Node number for 

node #4. 
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Columns Field Name Value Description 


11-20 

X 

1.0 

21-30 

Y 

1.0 

31-40 

Z 

1.0 

1-5 

NNUM 

5 

11-20 

X 

1.0 

21-30 

Y 

0.0 

31-40 

Z 

0.0 

1-5 

NNUM 

6 

11-20 

X 

0.0 

21-30 

X 

0.0 

31-40 

Z 

0.0 

1-5 

NNUM 

7 

11-20 

X 

0.0 

21-30 

Y 

0.0 

31-40 

Z 

1.0 

1-5 

NNUM 

8 

11-20 

X 

1.0 

21-30 

Y 

0.0 

31-40 

Z 

1.0 


X coordinate of node #4. 

Y coordinate of node # 4. 

Y coordinate of node # 4. 
Node number for node #5. 
X coordinate of node #5. 

Y coordinate of node #5. 

Y coordinate of node #5 . 
Node number for node #6. 
X coordinate of node #6. 

Y coordinate of node #6 . 

Y coordinate of node #6. 
Node number for node #7. 
X coordinate of node #7. 

Y coordinate of node #7. 

Y coordinate of node #7 . 
Node number for node # 8. 
X coordinate of node #8. 

Y coordinate of node #8 . 

Y coordinate of node #8. 
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3.2 ELEMENT TYPE 

There are four types of elements available in HITCAN, a plane stress element, a 
plane strain element, a plate element, and eight-node solid element. These 
elements are activated by program option cards. The table below lists the 
available elements and their corresponding program option cards. 


ELEMENT TYPE 

PROGRAM OPTION 
CARD 

Eght-Node Solid 

BRICK 

Four-Node Plate 

PLATE 

Plane Stress 

STRAIN 

Plane Strain 

STRESS 
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< 

BRIC K 

When this option is specified, a 8 -node isoparametric 3D solid element will be 
used. This option card is to be used in conjunction with the S3DS0LID and the 
READ IN MODEL option cards. This element has 3 translational degrees of freedom 
identified by u*, Uy, and u z , as shown in Figure 3.11. 



Figure 3.11: Eight-Node Solid Element 
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l 

PLAT E 

This option selects the four-node isoparametic plate element both for the solid 
( SPLATE model option ) and the hollow ( HPLATE model option ) structures . This 
element, which was derived from the Reissner-Mindlin theory for plates and 
shells, has 6 degrees -of -freedom at each node. These degrees -of- freedom are 
identified by u^, Uy, u z , H*, Hy, and H z . This element is shown in Figure 3.12. 



Figure 3.12: Plate Element 
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STRA IN 

This option specifies, that a four-node isoparametic plane strain element will 
be used in the analysis. Each node within this element has 2 degrees -of -freedom, 
identified by u* and Uy. This element is shown in Figure 3.13. 



Figure 3.13: Plane Strain Element 
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STRE SS 

This option specifies, that a four-node isoparametic plane stress element will 
be used in the analysis. Each node within this element has 2 degrees -of -freedom, 
identified by u* and Uy. This element is shown in Figure 3.14. 


u 


X 



EDGE TRACTION 
LOADING 


Y 


Figure 3.14: Plane Stress Element 
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33 TYPE OF ANALYSIS 

There are four types of analyses available in HITCAN: 

- Incremental static analysis with nonlinear anisotropic material behavior 

- Dynamic analysis using direct time integration 

- Modal analysis (free vibration) 

- Buckling analysis (first critical buckling load). 

The code will assume an incremental static analysis, unless the program option 
card DYNAMIC is specified. Note that for the static analysis no program option 
card or card groups are required. If a dynamic analysis is chosen by the user, 
four additional sets of data can be entered. They are: 

- Initial accelerations at selected nodes 

- Damping coefficients 

- Initial displacements at selected nodes 

- Initial velocities at selected nodes. 

One of the most features of HITCAN is it's capability to generate a restart file. 
This feature is useful in two ways . 

By using this feature the user can make several small runs instead of one large 
run, thus reducing turn around time. Secondly, if HITCAN fails to converge, a 
restart file is automatically generated. This enables the user to then continue 
the analysis at a smaller load increment. 

The table below lists the analyses and the restart feature with their 
corresponding program option card and card groups. 


TYPE OF ANALYSIS 

OPTIONS 

PROGRAM OPTIONS CARDS 

CARD GROUPS 

Static Dynamic 


DYNAMIC 

13 


Initial 

ACCELERATION 

22 


Acceleration 


39 


Damping 

DAMP 

14 


Initial 

DISPLACEMENT 

20 


Displacements 


37 


Initial 

VELOCITY 

21 


Velocity 


38 

Budding 


BUCKLE 

11 




44 

Modal 


MODAL 

11 




43 

Restart 


RESTART 

10 
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DYNAMIC ANALYSIS 


This option activates the direct time integration dynamic analysis. The 
integration scheme used is the Newmark-Beta method. For a description of this 
method, see Reference 7. 


Program option cards: DYNAM IC 

Card group : 13 

Variable 

Columns Format Name Entry 


1st card 
1-4 

14 

INCDYN 

Number of load increments between updating 

EXAMPLE: 

Card Group 13 
i 

2 3 

of the material properties . 
4 5 6 

i 

.0 

. .0 0. . . 

0 0 0 

i 

Columns 

Field 

Name Value 

Description 


1-4 INCDYN 1 The material properties will be updated every 

load increment. 
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INITIAL ACCELERATION 

This option enables the user to specify initial acceleration, at user-selected 
nodal points, for the direct time integration analysis. If this option is not 
specified, the initial accelerations are set to 0. 

Program option card: ACCEL ERATION 

Card groups: 22 and 39 

Variable 

Col limns Format Name Entry 

1st card of card group 22 

1-4 14 NACC Number of acceleration data sets . Each data 

set can specify acceleration for all degrees 
of freedom at several nodal points with equal 
increment in their node numbers . 

1st card of card group 39 

1-4 14 IBEG The first node number in this set of nodes. 

5-8 14 I END The last node in this set of nodes. 

9-12 14 INCR Increment at which the node numbers are 

incremented . 

2nd card of card group 39 

I- 10 F10.4 ACELIN(l) Initial acceleration for the 1st degree-of- 

freedom for this set of nodes . 

II- 20 F10.4 ACELIN(2) Initial acceleration for the 2nd degree-of- 

freedom for this set of nodes . 

21-30 F10.4 ACELIN(3) Initial acceleration for the 3rd degree-of- 

freedom for this set of nodes. 

31-40 F10.4 ACELIN(4) Initial acceleration for the 4th degree-of- 

freedom for this set of nodes. 

41-50 F10.4 ACELIN(5) Initial acceleration for the 5th degree-of- 

freedom for this set of nodes. 

51-60 F10.4 ACELIN(6) Initial acceleration for the 6th degree-of- 

freedom for this set of nodes. 

These two cards are repeated NACC times. 
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EXAMPLE: 

Card Group 22 

1 2 3 4 5 6 

1 0 0 0 0 0 0 

1 


Card Group 39 


1 

1 0 

2 

0. . . 


3 

0. . . 

4 5 6 

0 0 0 

4 7 

10.0 

1 

0.0 

0. 

0 

0.0 0.0 0.0 

Columns 

Field Name 

Value 

Description 

1-4 

NACC 

i 


Number of initial acceleration data 

1-4 

IBEG 

4 


sets is 1. 

The 1st node in this series of nodes is 1. 

5-8 

I END 

7 


Last node in this series of nodes is 7. 

9-12 

INCR 

1 


Increment at which the node numbers are 

1-10 

ACELIN(l) 

10.0 


incremented is 1. 

The 1st degree-of- freedom has an initial 

11-20 

ACELIN ( 2 ) 

0.0 


acceleration of 10.0 in. /sec. 2 

The 2nd degree -of -freedom has an initial 

21-30 

ACELIN ( 3 ) 

0.0 


acceleration of 0.0 in. /sec. 2 

The 3rd degree-of- freedom has an initial 

31-40 

ACELIN (4) 

0.0 


acceleration of 0.0 in. /sec. 2 

The 4th degree-of -freedom has an initial 

41-50 

ACELIN ( 5 ) 

0.0 


acceleration of 0.0 in. /sec. 2 

The 5 th degree-of -freedom has an initial 

51-60 

ACELIN ( 6 ) 

0.0 


acceleration of 0.0 in. /sec. 2 

The 6 th degree-of -freedom has an initial 


acceleration of 0.0 in. /sec. 2 
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I 


This option defines damping for the direct time integration analysis. At the 
present time, only Rayleigh clamping is available in HITCAN. Thus damping matrix 
[C] is of the form 


[C] - a[M] + 


where a is the damping coefficient 
coefficient of the stiffness matrix 

of the mass matrix [M] and B 

[KJ. 

is the damping 

Program option card: 

DAMP 



Card group: 

14 



Columns Format 

Variable 

Name 

Entry 


1st card 




1-8 F8.4 

9-16 F8.4 

DAMPMS 

DAMPST 

Damping coefficient for the mass matrix. 
Damping coefficient for the stiffness matrix. 


I 


March, 1992 | 


Chapter 3 


INPUT46 



HITCAN User's Manual - Version 1.0 


INITIAL DISPLACEMENT 


This option enables the user to specify initial displacement, at user-selected 
nodal points, for the direct time integration analysis. If this option is not 
specified, the initial displacements are set to 0. 


Program option card: DISP LACEMENT 

Card groups: 20 and 37 


Variable 

Columns Format Name Entry 

1st card of card group 20 

1-4 14 NDIS Number of displacement data sets. Each data 

set can specify displacement for all degrees 
of freedom at several nodal points with equal 
increment in their node numbers. 


1st card of card group 37 


1-4 

14 

IBEG 

The first node number in this 

set of nodes. 

5-8 

14 

I END 

The last node in this set of 

nodes . 

9-12 

14 

INCR 

Increment at which the node 
incremented . 

numbers are 


2nd card of card group 37 


1-10 

F10.4 

11-20 

F10.4 

21-30 

F10.4 

31-40 

F10.4 

41-50 

F10.4 

51-60 

F10.4 


DISPIN(l) 

DISPIN(2) 

DISPIN(3) 

DISPIN(4) 

DISPIN(5) 

DISPIN(6) 


Initial displacement for the 1st degree-of- 
freedom for this set of nodes. 

Initial displacement for the 2nd degree -of - 
freedom for this set of nodes. 

Initial displacement for the 3rd degree-of- 
freedom for this set of nodes . 

Initial displacement for the 4th degree-of- 
freedom for this set of nodes. 

Initial displacement for the 5th degree-of- 
freedom for this set of nodes . 

Initial displacement for the 6th degree-of- 
freedom for this set of nodes. 


These two cards are repeated NDIS times. 
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I 


4 5 6 

0 0 0 


4 5 6 

. .0 0 0 

0.0 0.0 0.0 


Description 


Number of initial acceleration data sets 
is 1 . 

The 1st node in this series of nodes is 1. 
Last node in this series of nodes is 14. 
Increment at which the node numbers are 
incremented is 2. 

The 1st degree-of-freedom has an initial 
displacement of 0.0 in. 

The 2nd degree-of-freedom has an initial 
displacement of 0.0 in. 

The 3rd degree-of-freedom has an initial 
displacement of 5.0 in. 

The 4th degree-of-freedom has an initial 
displacement of 0.0 in. 

The 5th degree-of-freedom has an initial 
displacement of 0.0 in. 

The 6th degree-of-freedom has an initial 
displacement of 0.0 in. 
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INITIAL VELOCITY 

This option enables the user to specify initial velocities, at user-selected 
nodal points, for the direct time integration analysis. If this option is not 


specified, i 

the initial 

velocities 

are set to 0. 

Program option card: 

VELOCITIES 

Card groups 


21 and 38 
Variable 

Columns 

Format 

Name 

Entrv 

1st card of 

card group 

21 


1-4 

14 

NVEL 

Number of velocity data sets . Each data set 
can specify the velocity for all degrees of 
freedom at several nodal points with equal 
increment in their node numbers. 

1st card of 

card group 

38 


1-4 

14 

IBEG 

The first node number in this set of nodes. 

5-8 

14 

IEND 

The last node in this set of nodes. 

9-12 

14 

INCR 

Increment at which the node numbers are 
incremented. 

2nd card of 

card group 

38 


1-10 

F10.4 

VEL0IN(1) 

Initial velocity for the 1st degree-of- 
freedom for this set of nodes. 

11-20 

F10.4 

VELOIN ( 2 ) 

Initial velocity for the 2nd degree-of- 
freedom for this set of nodes. 

21-30 

F10.4 

VELOIN ( 3 ) 

Initial velocity for the 3rd degree-of- 
freedom for this set of nodes. 

31-40 

F10.4 

VELOIN (4) 

Initial velocity for the 4th degree-of- 
freedom for this set of nodes. 

41-50 

F10.4 

VEL0IN(5) 

Initial velocity for the 5th degree-of- 
freedom for this set of nodes. 

51-60 

F10.4 

VEL0IN(6) 

Initial velocity for the 6th degree-of- 
freedom for this set of nodes. 

These two cards are repeated NVEL 

times . 
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EXAMPLE: 

Card Group 21 

1 2 3 4 5 6 

1 0 0 0 0 0 0 

1 


Card Group 38 


1 

, . . 0 . . . 
7 1 

1.0 


2 

, 0 . 


3 

. 0 . 


4 

. 0 . 


5 

. 0 . 


6 

.0 


0.0 


0.0 


0.0 


0.0 


0.0 


Columns 

Field Name 

Value 

Description 


1-4 

NVEL 

1 

Number of initial velocity 

data sets is 1. 

1-4 

IBEG 

4 

The 1st node in this series 

of nodes is 1. 

5-8 

I END 

7 

Last node in this series of 

nodes is 7. 

9-12 

INCR 

1 

Increment at which the node numbers are 
incremented is 1 . 

1-10 

VELOIN(l) 

10.0 

The 1st degree -of- freedom 
velocity of 1.0 in. /sec. 

has an initial 

11-20 

VELOIN ( 2 ) 

0.0 

The 2nd degree -of- freedom 
velocity of 0.0 in. /sec. 

has an initial 

21-30 

VEL0IN(3) 

0.0 

The 3rd degree-of- freedom 
velocity of 0.0 in. /sec. 

has an initial 

31-40 

VELOIN (4) 

0.0 

The 4th degree -of- freedom 
velocity of 0.0 in. /sec. 

has an initial 

Columns 

Field Name 

Value 

Description 


41-50 

VEL0IN(5) 

0.0 

The 5th degree -of- freedom 
velocity of 0.0 in. /sec. 

has an initial 

51-60 

VELOIN ( 6 ) 

0.0 

The 6th degree -of- freedom 
velocity of 0.0 in. /sec. 

has an initial 
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BUCKLING ANALYSIS 

This option activates the buckling analysis. To determine the critical buckling 
load, MHOST uses the subspace iteration method. For a description of this method, 
see Reference 7. When this option is used, a buckling analysis is performed at 
the initial load and the times specified in the array TIMEMB. 

Program option card: BUCK LE 

Card groups: 11 and 44 

Variable 

Columns Format Name Entry 

1st card of card group 11 

14 NEIGV Number of critical buckling modes to be 

extracted. At the present time this value 
should be set to 1. 

14 NSUBD Number of subspace dimensions for the 

eigenvalue extraction. A sufficient 
dimension is reserved as a default. 

14 MHITER Maximum number of iterations allowed in MHOST 

for the subspace iteration. The default for 
the maximum number of iterations is 5. 

2nd card of card group 11 

1-8 F8.2 RESID The allowable tolerance in MHOST. The default 

value is 5.0. 

1st card of card group 43 

1-10 F10.6 TIMEMB Times at which a buckling analysis is 

desired. A maximum of 8 times are allowed. 
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EXAMPLE: 

Card Group 11 

1 2 3 4 5 6 

1 0 0 0 0 0 0 

1 0 10 

5.0 


Card Group 43 

1 2 3 4 5 6 

1 0 0 0 0 0 0 

10.0 


Columns 

Field Name 

Value 

Description 


1-4 

NEIGV 

i 

Must be set to 1. 


5-8 

NSUBD 

0 

The default value is used. 


9-12 

MHITER 

10 

The allowable number of iterations 
is set to 10. 

in MHOST 

1-8 

RES ID 

5.0 

The allowable tolerance in MHOST is 

set 5.0. 

1-10 

TIMEMB 

10.0 

A buckling analysis is desired at 

10 sec. 
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MODAL ANALYSIS 


When this option is used HITCAN will perform a free vibration dynamic analysis 
to determine the frequencies and the mode shapes. To determine the frequencies 
and the mode shapes, MHOST uses the subspace iteration method. For a description 
of this method, see Reference 7. 


Program option card: MODAL 

Card groups: 11 and 43 


Variable 


Columns 


Format 

Name 

Entry 

1st card 

of 

card group 

11 


1-4 


14 

NEIGV 

Number of modes to be extracted. At the 
present time this value should be set to 1. 

5-8 


14 

NSUBD 

Number of subspace dimensions for the 
eigenvalue extraction. A sufficient 
dimension is reserved as a default. 

9-12 


14 

MHITER 

Maximum number of iterations allowed in MHOST 
for the subspace iteration. The default for 
the maximum number of iterations is 5. 

2nd card 

of 

card group 

11 


1-8 


F8.2 

RESID 

The allowable tolerance in MHOST. The 

default value is 5.0. 

1st card 

of 

card group 

43 


1-10 


F10.6 

TIMEMS 

Times at which a modes are to be extracted. 


A maximum of 8 times are allowed. 
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< 

EXAMPLE: 

Card Group 11 

1 2 3 4 5 6 

1 0 0 0 0 0 0 

1 0 10 

5.0 


Card Group 43 

1 2 3 4 5 6 

1 0 0 0 0 0 0 

5.0 


Columns 

Field Name 

Value 

Description 

1-4 

NEIGV 

i 

Must be set to 1. 

5-8 

NSUBD 

0 

The default value is used. 

9-12 

MHITER 

10 

The allowable number of iterations in MHOST 




is set to 10. 

1-8 

RESID 

5.0 

The allowable tolerance in MHOST is set 5.0. 

1-10 

TIMEMB 

10.0 

A buckling analysis is desired at 5 sec. 
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RESTART 

This feature enables the user to conduct an analysis in several runs, an useful 
option for large problems. When the RESTART option is specified, a restart file 
created in a previous run is input. This file contains the necessary information 
to continue the analysis, including, the load step number, increment number, ply 
stresses, microstresses, microstress rates, material failure flags, and nodal 
displacements. For a restart run, the input is the same as before, except that 
RESTART is now specified in the program option cards. 

A restart file is created by specifying the variable MSTART. MSTART is the number 
of increments to be preformed in a particular run. After the analysis has 
progressed through MSTART increments, a restart file is created and the run is 
terminated. Note that a restart file will also be created when the maximum 
number of allowable iterations ( the variable MITER ) in HITCAN is exceeded. 
Note that card group #10 must be entered whether or not the program option card 
RESTART is specified. 

Program option card: RESTA RT 

Card group : 10 

Variable 

Col unrns Format Name Entry 


1st card in card group 10 


1-4 

14 

NTISTP 

Number of load steps. 

5-8 

14 

NMECH 

Number of mechanical cycles used in METCAN 
to account for cyclic damage (used for 
fatigue analysis) . 

9-12 

14 

NTHER 

Number of thermal cycles used in METCAN to 
account for cyclic damage (used for fatigue 
analysis). 

13-16 

14 

LINC 

Number of load increments between load steps , 
see following figure. 

17-20 

14 

MSTART 

Write a restart file after this many load 
increments . 

21-24 

14 

MITER 

Maximum number of iterations allowed for 
global convergence per load increment. 



Variable 


Columns 

Format 

Name 

Entry 

2nd card 




1-8 

F8.4 

TOL 

Allowable global tolerance on convergence 

Chapter 3 



March, 1992 


INPUTS 5 


EXAMPLE: 


HITCAN User's Manual - Version 1.0 


I 


Card Group 10 

1 2 3 4 5 6 

1 0 0 0 0 0 0 

1 1 1 1 2 10 

1.0 


Columns 

Field Name 

Value 

Description 


1-4 

NTISTP 

i 

Number of load steps is 1. 


5-8 

NMECH 

i 

Number of mechanical cycles is 1. 


9-12 

NTHER 

i 

Number of thermal cycles is 1. 


13-16 

LINC 

i 

Number of load increments is 1. 


17-20 

MSTART 

2 

Write a restart file after 2 load increments. 

21-24 

MITER 

10 

Number of allowable iterations for 
convergence. 

global 

1-8 

TOL 

1.0 

Tolerance on the global convergence 

is 1.0. 


< 
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3.4 MATERIAL DATA 

In this section, the necessary program option cards and card groups used in 
building the composite model from constituent properties are described. The user 
builds the composite model by defining the plies over the surface of the 
structure and through the thickness. This produces an integrated laminated model 
of the entire structure. 

The user has the flexibility of selecting the number of plies, the ply thickness, 
the fiber volume ratio, the void volume ratio, and the fiber orientation in 
building either a symmetric or an unsymmetric ply layup. When PLYORDER is used 
in conjunction with the program option card, UNSYMMETRICAL, the user can generate 
an unsymmetric ply layup. However, when the HPLATE model option is specified, 
only symmetric ply layups can be used. 

Each ply is referred to by a designation number determined by the order in which 
it is listed in the input. A different ply designation number must be assigned 
if any one or more of the following variables; material properties, ply 
thickness, fiber volume ratio, void volume ratio, and fiber orientation, are 
different. However, a ply used more than once, but with same value to all of 
these variables, can be assigned only one ply designation number. 

For a symmetric ply layup, the ply order is assumed to be from the bottom surface 
of the structure to the mid- thickness line moving in the positive y or z 
directions . This order is then reversed and the ply layup is then continued from 
the mid- thickness line to the opposite surface of the wall. If the ply layup is 
unsymmetric, i. e. , the program option card UNSYMMETRICAL is specified, the ply 
layup is input in two parts. The first part is the ply layup from the bottom 
surface of the structure to the mid-thickness line moving in the positive y or 
z directions. The second part is the ply order from the opposite surface to the 
mid- thickness line moving in the negative y or z directions. 

The ply orientation angle, measured in degrees, is defined from the HITCAN global 
x-axis to be positive in the positive x-y quadrant for the SPLATE, HPLATE, 
S3DS0LID, and READ IN MODEL options. For the spars of the HPLATE model option 
the angle is positive in the positive x-z quadrant. Figure 3.15 illustrates the 
composite geometry relative to the HITCAN coordinate system. Note at the present 
time only one ply description is allowed for each spar. A sufficient number of 
these plies will be generated to fill the spar thickness. The orientation of the 
fibers in these plies will be along the x-axis. 

To account for the existence of a discrete interphase between the fiber and the 
matrix, the user is allowed to enter the interphase thickness as a fraction of 
the fiber diameter. This is done by specifying the program option card 
INTERPHASE and entering value of the variable PINTER in card group 23. 
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In creating this section of input, the card groups 1 and 28 must always be 
specified. These two card groups contain the necessary data to describe the 
plys. The program option cards PLYORDER and UNSYMMETRICAL can be used only with 
the model options HPLATE, SPLATE, and S3DS0LID. For the model option READ IN 
MODEL, the user must use card groups 12 and 31. The table below summarizes the 
program option cards and card groups used in this section. 



PROGRAM 
OPTION CARDS 

CARD GROUPS 

COMMENTS 

Ply Description 


1 

This data is needed 




for all model options. 

Data 


28 


Symmetrical 


6 

Use for SPLATE, 


PLYORDER 


HPLATE, & S3DS0LID 

Ply Order 


29 

model options. 

Unsymmetrical 

PLYORDER 

7 

Use for SPLATE, 




& S3DSOLID 

Ply Order 

UNSYMMETRICAL 

30 

model options. 

Ply Layup 






8 

Use for the 

Specified at 


12 

model option 

Each Node 


31 

READ IN model. 




Used only if an 

Interphase 

INTERFACE 

23 

interface between 




the fiber and 

Data 



matrix is to be 




modelled. 
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Figure 3.15: Fiber Composite Geometry Relative To HITCANCoordinate System 
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PLY DESCRIPTION CARDS 

These cards contain the necessary data to describe each unique ply in the 
structure. These card groups are required for all model options. 


Program option card: 

None 


Card group: 

1 and 28 



Variable 


Columns 

Format 

Name 

Entrv 

1st card 

in card group 

1 


1-4 

14 

NDES 

Number of ply descriptions. 

1st card 

in card group 

28 


1-10 

F10.4 

PERT(l) 

Initial thickness (percent of thickness at 




each input point) . 

11-20 

F10.4 

PERT ( 2 ) 

Final thickness (percent of thickness at each 




input point) . 

1-30 

F10.4 

PERT ( 3 ) 

Initial x coordinate (percent length) . 

1-40 

F10.4 

PERT(4) 

Final x coordinate (percent length) . 

41-50 

F10.4 

PERT ( 5 ) 

Initial y coordinate (percent width) . 

51-60 

F10.4 

PERT ( 6 ) 

Initial y coordinate (percent width) . 

2nd card 

in card group 

28 


1-4 

A4 

CODES (1) 

Type of fiber. 

5 

IX 



6-9 

A4 

CODES (2) 

Type of matrix. 

10 

IX 



11-20 

F10.4 

CODES (3) 

Ply thickness. 

21-30 

F10.4 

CODES (4) 

Void volume ratio. 

31-40 

F10.4 

CODES (5) 

Fiber volume ratio. 

41-50 

F10.4 

CODES (6) 

Ply orientation angle. 

These two cards are required for 

each ply description, thus these cards are 


repeated NDES times. Ply designation numbers are assigned in the order in which 
the ply descriptions are listed. 


Chapter 3 


March, 1992 


INPUT60 



HITCAN User's Manual - Version 1.0 


EXAMPLE: 


Card Group 1 


1 

i n 

eg c 

3 

0. . . 

4 5 6 

0 0 0 

1 

Card Group 28 



1 

2 

3 

4 5 6 

1 0 

. .0. . . 

0. . . 

0 0 0 

0.0 

100.0 

0.0 

100. 0.0 100.0 

SICA TI15 

0.2 

0.0 

0.5 45.0 

Columns 

Field Name 

Value 

Description 

1-4 

NDES 

i 

The number of plys to be described is 1. 

1-10 

PERT ( 1 ) 

0.0 

The initial thickness is 0.0X of the 
thickness at each input point. 

11-20 

PERT ( 2 ) 

100.0 

The final thickness is 100. 0Z of the 
thickness at each input point. 

21-30 

PERT ( 3 ) 

0.0 

The initial x coordinate of the ply is at 
0.0X of the length at each input point. 

31-40 

PERT(4) 

100.0 

The final x coordinate of the ply is at 
100. OX of the length at each input point. 

41-50 

PERT ( 5 ) 

0.0 

The initial y coordinate of the ply is at 
0.0X of the width at each input point. 

51-60 

PERT ( 6 ) 

100.0 

The final y coordinate of the ply is at 
100. OX of the width at each input point. 

1-4 

CODES (1) 

SICA 

The fiber material is designated in the 
databank as SICA. 

6-9 

CODES (2) 

TI15 

The matrix material is designated in the 
databank as T115 . 

11-20 

CODES (3) 

0.2 

The ply thickness is 0.2 in. 

21-30 

CODES (4) 

0.0 

The void volume ratio is 0.0. 

31-40 

CODES (5) 

0.5 

The fiber volume ratio is 0.5. 

41-50 
Chapter 3 

CODES (6) 

45.0 

The fiber orientation angle is 45 
degrees . 
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PLY LAYUP FOR A SYMMETRICAL PLYORDER 

This program option card allows the user to create a symmetrical ply layup. When 
PLYORDER is used in conjunction with the program option card, UNSYMMETRICAL, the 
user can generate an unsymmetric ply layup. However, when the HPLATE model 
option is specified, only symmetric ply layups can be used. 

Program option card: PLYOR DER 

Card group: 6 and 29 


Variable 

Columns Format Name Entry 

1st card in card group 6 

1-4 14 MAXPLY Number of plies specified for the half- 

thickness at the point of maximum wall 
thickness in the structure for the model 
options HPLATE, SPLATE, and S3DS0LID. 

1st card in card group 29 

1-80 2014 MPLY An array containing the ply order , using ply 

designation numbers, for one half of the wall 
thickness. The order is from the bottom 
surface to the mid- thickness line moving in 
the positive y or z directions. 


EXAMPLE: 

Card Group 6 

1 2 3 4 5 6 

1 0 0 0 0 0 0 

2 


Card Group 29 

1 2 3 4 5 6 

1 0 0 0 0 0 0 

1 1 
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Columns 

Field Name 

Value 

1-4 

MAXPLY 

2 

1-4 

MPLY(l) 

1 

5-8 

MPLY ( 2 ) 

1 


Description 

There are two plys from the bottom 
surface to the mid- thickness line. 
The 1st ply up from the bottom surface 
has the designation number of 1. 

The 2nd ply up from the bottom surface 
has the designation number of 1. 


Chapter 3 


March, 19912 


INPUT63 


HITCAN User's Manual - Version 1.0 


PLY LAYUP FOR AN UNSYMMETRICAL PLYORDER 

This option specifies an unsymmetric ply layup. This option is to be used in 
conjunction with the PLYORDER program option card. 

Program option card: UNSYM METRICAL PLYORDER and 

PLYOR DER 

Card group: 7 and 30 


Variable 

Columns Format Wants Entry 

1st card of card group 7 

1-4 14 LMAX Number of plies specified for the half- 

thickness at the point of maximum wall 
thickness in the structure for the model 
options SPLATE and S3DS0LID. 

2nd card of card group 30 

1-80 2014 NPLY An array containing the ply order, using ply 

designation numbers, for one half of the wall 
thickness . The order is from the top surface 
to the mid- thickness line moving in the 
negative y or z directions. 


EXAMPLE: 

Card Group 7 

1 2 3 4 5 6 

1 0 0 0 0 0 0 

2 


Card Group 30 

1 2 3 4 5 6 

1 0 0 0 0 0 0 

1 2 
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Columns 

Field Name 

Value 

1-4 

MAXPLY 

2 

1-4 

MPLY(l) 

1 

5-8 

MPLY(2) 

2 


Description 

There are two plys from the bottom 
surface to the mid- thickness line. 
The 1st ply up from the mid- thickness 
line has the designation number of 1. 
The 2nd ply up from the mid- thickness 
line has the designation number of 1. 
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PLY LAYUP FOR READ IN MODEL OPTION 

This option enables the user to enter a ply layup for the READ IN MODEL option. 

Program option card: None 

Card group: 8, 12, and 31 


Variable 


Columns 

Format 

Name 

Entrv 

1st card 

of card 

group 

8 


1-4 

14 


MAXPLY 

Maximum number of plies at any one node. 

1st card 

in card 

group 

12 


1-4 

14 


NPLSET 

Number of data sets describing the ply layup 





at each node. 

1st card 

in card 

group 

31 


1-4 

14 


IBEG 

The first node number in this set of nodes. 

5-8 

14 


I END 

The last node number in this set of nodes. 

9-12 

14 


INCR 

Increment at which the node numbers are 





incremented . 

13-16 

14 


NOP 

Number of plies at these nodes. 

2nd card 

of card 

group 

31 


1-80 

2014 


MPLY 

An array containing the ply order, using ply 


designation numbers , for one half of the wall 
thickness. The order is from the bottom 


surface to the mid- thickness line moving in 
the positive y or z directions. 


These 2 cards are repeated NPLSET times . 
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EXAMPLE: 


Card Groun 8 





1 2 

3 

4 

5 6 

1 

3 

.0 0 

0 

. .0. 

0 0 

Card Group 12 





1 2 

3 

4 

5 6 

i 

i 

.0 0 

0 

. .0. 

0 0 

Card Grout) 31 





1 2 

3 

4 

5 6 

i 

.0 0 

0 

. .0. 

0 0 

1 16 

1 3 




1 2 

1 




Columns 

Field Name 

Value 


Description 

1-4 

MAXPLY 

2 


The maximum number of plies is 3. 

1-4 

NPLSET 

1 


One data set is used to describe the 
ply layup or each node in the model. 

1-4 

IBEG 

1 


The 1st node in this series of nodes 
is 1. 

5-8 

IEND 

16 


The last node in this series of nodes 
is 16. 

9-12 

INCR 

1 


The increment between node numbers 
is 1 . 

13-16 

NOP 

3 


The number of plies to be read in 
is 3. 

1-4 

MPLY(l) 

1 


The 1st ply up from the bottom surface 
line has the designation number of 1. 

5-8 

MPLY ( 2 ) 

2 


The 2nd ply up from the bottom surface 
has the designation number of 2. 

9-12 

MPLY ( 3 ) 

1 


The 3rd ply up from the bottom surface 
has the designation number of 1. 
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i 

INTERPHASE DATA 


This option enables the user to specify the existence of a discrete interphase 
between the fiber and the matrix. The thickness of the interphase is specified 
as a fraction of the fiber diameter, i.e., Pinter = (D 0 - D)/D 0 , where D 0 
represents the initial fiber diameter and D the fiber diameter after degradation. 

Program option card: INTE RPHASE 

Card group: 23 

Variable 

Columns Format Name Entry 

1st card of card group 23 

1-8 F8.4 PINTER Thickness of the interphase as a fraction of 

fiber diameter. 


EXAMPLE: 


Card Group 23 

1 2 3 4 5 6 

1 0 0 0 0 0 0 

0.1 


I 


Columns 

Field Name 

Value 

Description 

1-8 

PINTER 

.i 

The thickness of the interphase is . 1 


x fiber diameter. 
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3.5 TYPES OF LOADING 

There are four types of loading available in HITCAN. They are: 

Centrifugal , 

Nodal Forces, 

Pressure , 

Temperature . 

Each loading type has its own program option card and card groups. Along with 
the program option cards and card groups associated with the different loadings , 
additional card groups are needed for program control. The following table 
summarizes the different loadings and their program option cards and card 
groups . 


LOADING TYPE 

PROGRAM OPTION 
CARD 

CARD GROUP 

Centrifugal 

ANGULAR 

35 

Nodal Forces 

FORCE 

15 



36 



16 



33 

Pressure 

PRESSURE 

34 

Temperature 

TEMPERATURE 

17 



33 



34 

Program 


10 

Control 


32 
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CENTRIFUGAL LOADING 

This option defines a centrifugal distributed loading. 

Program option card: ANGU LAR 

Card group: 35 

Variable 

Columns Format Name Entry 

1st card 

1-30 3F10.4 GRIDP1 The first of two points required to define 

the axis of rotation. See Figure 3.16, for 
the direction of this axis. GRIDP1 is an 
array of three real numbers . The first real 
number is the x coordinate, the second is the 
y coordinate, and the third real number is 
the z coordinate. 

31-60 3F10.4 GRIDP2 The second of two points required to define 

the axis of rotation. GRIDP2 is an array of 
three real numbers. The first real number 
is the x coordinate, the second is the y 
coordinate, and the third real number is the 
z coordinate. 

2nd card 

1-80 8F10.4 ANGVEL Rotational speed in revolutions per second 

at each time step. 
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Figure 3.16: Direction Of Axis For Centrifugal Loading 
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I 


Card Group 35 

1 2 3 4 5 6 

1 0 0 0 0 0 0 

0.0 0.0 0.0 1.0 0.0 0.0 

100.0 


Columns 

Field Name 

Value 

1-10 

GRIDPl(l) 

0.0 

11-20 

GRIDP1(2) 

0.0 

21-30 

GRIDP1(3) 

0.0 

31-40 

GRIDP2 (1) 

1.0 

41-50 

GRIDP2(2) 

0.0 

51-60 

GRIDP2 (3) 

0.0 

1-10 

ANGVEL(l) 

100.0 


Description 

X- coordinate 
Y- coordinate 
Z- coordinate 
X- coordinate 
Y- coordinate 
Z- coordinate 

The rotational velocity in rad/sec. 
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NODAL FORCE LOADING 

This option allows the user to input static concentrated nodal ’loads. 
Program option card: FORC E 

Card groups: 15 and 36 


Variable 

Columns Format Name 

1st card in card group 15 

1-4 14 NCFOR 

1st card in card group 36 

1-4 14 NCFNOD 

5-8 14 NCFDIR 

2nd card in card group 36 

1-80 8F10.4 CFVAL 


Entry 


Number of nodal loads 


Node number 
Degree -of- freedom 


Value of the load at each load step 


Card group 36 is repeated NCFOR times. 
EXAMPLE: 


Card Group 15 


1 


2 

3 

4 

5 

6 

1 . .0. 


.0 

0 

. .0. 

0 

. . .0 

1 







Card Group 36 






1 


2 

3 

4 

5 

6 

1 0. 


.0 

0 

. .0. 

0 

...0 

1 1 







10. 







Columns 

Field Name 

Value 


Description 


1-4 

NCFOR 


i 


Number of nodal forces is 1. 

1-4 

NCFNOD 


i 


The force acts 

at node #1. 

5-8 

NCFDIR 


i 


The force acts 

along the X-a: 

1-10 

CFVAL 


10.0 


The intensity 

of the load. 
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( 

PRESSURE LOADING 


This card set allows the user to input the pressure loading. For either solid or 
hollow structures using the plate element, the user can specify both a lateral 
pressure and an uniform edge pressure. For the plane stress, plane strain, and 
3D solid elements only an uniform edge pressure can be specified. Note that if 
PL and PU are used to input the pressure loading the variables TL and TU must 
also be specified. If temperature effects are not desired, then TL and TU must 
be set to 0.0 

The pressure is entered using one of the following three options: 

A. When the SPLATE model option is specified. 

B. When the HPLATE model option is specified. 

C. When the S3DS0LID model or the READ IN MODEL options are 
specified. 


Program option card: PRES SURE 

Card groups: 16, 33, and 34 


< 


Variable 

Col tunns Format Name Entry 

1st card in card group 16 

1-4 14 NPRES Number of pressure input data sets. These 

data sets specify those elements which have 
a normal surface pressure (positive into 
the element) . 

Option A 


card la in card group 33 


1-10 

F10.4 

TL 

Temperature 

on the lower surface. 

11-20 

F10.4 

TU 

Temperature 

on the upper surface. 

21-30 

F10.4 

PL 

Pressure on 

the lower surface. 

31-40 

F10.4 

PU 

Pressure on 

the upper surface. 


This card is repeated (MSECT(l) + MSECT(2) +. . .+ MSECT(I) )*NTISTP times, where 
I-l.NSECT. 
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Columns 

Format 

Name 

Entrv 

Option B 

card lb in 

card group 

33 


1-10 

F10.4 

TL 

Temperature on the lower surface. 

11-20 

F10.4 

TU 

Temperature on the upper surface. 

21-30 

F10.4 

PU 

Pressure on the lower surface. 

31-40 

F10.4 

PL 

Pressure on the upper surface. 

Note that each card corresponds 

to an input point, thus, this card is repeated 

(MSECT(1,1) 

+ MSECT(2 , 1) +...+ MSECT(I ,1) )*NTISTP times, where 1-1 ,LSECT(1) . 

card 2b in 

card group 

33 


1-10 

F10.4 

TL 

Temperature on the lower surface. 

11-20 

F10.4 

TU 

Temperature on the upper surface. 

21-30 

F10.4 

PU 

Pressure on the lower surface. 

31-40 

F10.4 

PL 

Pressure on the upper surface. 


Note that each card corresponds to an input point, thus, this card is repeated 
(MSECT(1 , 2) + MSECT(2 , 2) +...+ MSECT(I,2))*NTISTP times, where 1-1 ,LSECT(2) . 


card 3 in card group 33 


1-4 

14 

IBEG 

5-8 

14 

I END 

9-12 

14 

INCR 

card 4 

in card group 

33 

1-80 

F10.4 

PREW 


The first node number in this set of nodes. 
The last node in this set of nodes. 
Increment between node numbers. 


An array containing the intensity of the edge 
pressure at each load step. 


Cards 3 and 4 are repeated NPRES times. 


Option C 

card 3 in card group 34 


1-4 

14 

IBEG 

Beginning element number. 

5-8 

14 

I END 

Ending element number. 

9-12 

14 

INCR 

Increment between element numbers. 


card 4 in card group 34 

1-80 F10.4 PREVAL An array containing the intensity of the edge 

pressure at each load step. 


Cards 3 and 4 are repeated NPRES times. 
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Cards 3 and 4 

are repeated NPRES 

times . 



EXAMPLE: 






Card Grout) 16 
1 

2 

3 

4 

5 

6 

1 0 

0 

0. . . 

0 

0 

0 

1 






Card Grout) 33 
1 

2 

3 

4 

5 

6 

1 0 

0 

0. . . 

0 

0 

0 

70.0 

70.0 

0.0 

10.0 



70.0 

70.0 

0.0 

10.0 



70.0 

70.0 

0.0 

20.0 



70.0 

70.0 

0.0 

20.0 



Card Grout) 34 
1 

2 

3 

4 

5 

6 

1 0 

0 

0. . . 

0 

0 

0 


14 1 

10.0 


Columns 

Field Name 

Value 

Description 

1-4 

NPRES 

i 

There is 1 set of data to describe the 
edge loads. 

1-10 

XL 

70.0 

Temperature on the lower surface at 
input point #1 of input section #1. 

11-20 

TU 

70.0 

Temperature on the lower surface at 
input point #1 of input section #1. 

21-30 

PL 

0.0 

Pressure on the lower surface at input 
point #1 of input section #1. 

31-40 

PU 

10.0 

Pressure on the upper surface at input 
point if 1 of input section ifl. 

1-10 

TL 

70.0 

Temperature on the lower surface at 
input point if 2 of input section ifl. 

11-20 

TU 

70.0 

Temperature on the lower surface at 
input point #2 of input section ifl. 

21-30 

PL 

0.0 

Pressure on the lower surface at input 
point ifl of input section if 1. 

31-40 

PU 

10.0 

Pressure on the upper surface at input 
point ifl of input section if 1. 

1-10 

TL 

70.0 

Temperature on the lower surface at 
input point ifl of input section ifl. 

11-20 

TU 

70.0 

Temperature on the lower surface at 
input point ifl of input section ifl. 

21-30 
Chapter 3 

PL 

0.0 

Pressure on the lower surface at input 
point ifl of input section ifl. 
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Columns 

Field Name 

Value 

Descriotion 

31-40 

PU 

20.0 

Pressure on the upper surface at input 

1-10 

TL 

70.0 

point #1 of input section on #2. 
Temperature on the lower surface at 

11-20 

TU 

70.0 

input point #2 of section #2. 
Temperature on the lower surface at 

21-30 

PL 

0.0 

input point #2 of input section #2 . 
Pressure on the lower surface at input 

31-40 

PU 

20.0 

point n of input section #2 . 
Pressure on the upper surface at input 

1-4 

IBEG 

1 

point n of input section #2. 

The edge load is applied to elements 

1 

5-8 

I END 

4 

through 4. 

9-12 

INCR 

1 


1-10 

PREVAL 

10.0 

The intensity of the edge load. 
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TEMPERATURE LOADING 


This card set allows the user to input a distributed temperature loading. Note 
that if TL and TU are used to input the temperature loading the variables PL and 
PU must also be specified. If there is no lateral pressure loading, then PL and 
PU should be set 0.0. The ply temperatures generated by TL and TU will be 
overridden by the ply temperatures contained in the variable TEPLY. The 
temperature is entered using one of the following three options: 

A. When the SPLATE model option is specified. 

B. When the HPLATE model option is specified. 

C. When the S3DS0LID model or the READ IN MODEL options are 
specified. 


Program option card: TEMP ERATURE 

Card groups: 17, 33, and 34 


Variable 

Columns Format Name Entry 

1st card in card group 17 


1-4 14 NTEMP Number of temperature input data sets . These 

data sets be can used either to specify the 
temperature at the nodes (READ IN MODEL or 
S3DS0LID model options) or at the plies at 
each node (HPLATE or SPLATE model options) . 

5-8 14 NUMPLY The maximum number of plies in a data set. 
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Variable 

Columns Format Name Entry 

Option A 


card la in card group 33 


1-10 

F10.4 

TL 

Temperature 

on the lower surface. 

11-20 

F10.4 

TU 

Temperature 

on the upper surface. 

21-30 

F10.4 

PL 

Pressure on 

the lower surface. 

31-40 

F10.4 

PU 

Pressure on 

the upper surface. 


Note that each card corresponds to an input point, thus, this card is repeated 
(MSECT(l) + MSECT(2) + . . .+ MSECT(I) )*NTISTP times, where I-1,NSECT. 

Option B 


card lb in card group 33 


1-10 

F10.4 

TL 

Temperature 

on the lower surface. 

11-20 

F10.4 

TU 

Temperature 

on the upper surface. 

21-30 

F10.4 

PU 

Pressure on 

the lower surface. 

31-40 

F10.4 

PL 

Pressure on 

the upper surface. 


This card is repeated (MSECT(1,1) + MSECT(2,1) + . . .+ MSECT(I,1))* NTISTP times, 
where I-1,LSECT(1) . 


card 2b in card group 33 


1-10 

F10.4 

TL 

Temperature 

on the lower surface. 

11-20 

F10.4 

TU 

Temperature 

on the upper surface. 

21-30 

F10.4 

PU 

Pressure on 

the lower surface. 

31-40 

F10.4 

PL 

Pressure on 

the upper surface. 


This card is repeated (MSECT(1,2) + MSECT(2,2) + ... + MSECT(1 , 2) )*NTISTP times, 
where I-1,LSECT(2) . 


card 3 in card group 33 


1-4 

14 

IBEG 

5-8 

14 

I END 

9-12 

14 

INCR 


The first node number in this set of nodes. 
The last node in this set of nodes. 
Increment between node numbers. 


card 4 in card group 33 

1-80 F10.4 TEPLY Array of ply temperatures at each ply for 

each time step. 


Cards 3 and 4 are repeated NTEMP times. 
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card 

3 

in 

card group 

34 


1-4 



14 


IBEG 

5-8 



14 


I END 

9-12 



14 


INCR 

card 

4 

in 

card group 

34 


1-80 



F10.4 


TENOD 


I 


Beginning node number . 

Ending node number. 

Increment between node numbers . 


Array of nodal temperatures at each 
node in the data set for each load 
step . 


Cards 3 and 4 are repeated NTEMP times. 


EXAMPLE: 


Card Group 33 


1 

2 

3 

4 

1 0 . 

0 

0 

. . .0. 

70.0 

70.0 

0.0 

0.0 

70.0 

70.0 

0.0 

0.0 

100.0 

100.0 

0.0 

0.0 

100.0 

100.0 

0.0 

0.0 

Columns 

Field Name 

Value 


1-10 

TL 

70.0 


11-20 

TU 

70.0 


21-30 

PL 

0.0 


31-40 

PU 

0.0 


1-10 

TL 

70.0 


11-20 

TU 

70.0 


21-30 

PL 

0.0 


31-40 

PU 

0.0 



5 6 

0 0 


Description 

Temperature on the lower surface at 
input point #1 of input point section 
# 1 . 

Temperature on the lower surface at 
input point #1 of input section #1. 
Pressure on the lower surface at input 
point #1 of input section #1. 
Pressure on the upper surface at input 
point #1 of input section #1. 
Temperature on the lower surface at 
input point #2 of input section #1. 
Temperature on the lower surface at 
input point #2 of input section #1. 
Pressure on the lower surface at input 
point #2 of input section #1. 
Pressure on the upper surface at input 
point #2 of input section #1. 
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Columns Field Name Value Description 


1-10 

TL 

100.0 

11-20 

TU 

100.0 

21-30 

PL 

0.0 

31-40 

PU 

0.0 

1-10 

TL 

100.0 

11-20 

TU 

100.0 

21-30 

PL 

0.0 

31-40 

PU 

0.0 


Temperature on the lower surface at 
input point # 1 of input section #2. 
Temperature on the lower surface at 
input point #1 of input section #2. 
Pressure on the lower surface at intip 
point #1 of input section # 2 . 

Pressure on the upper surface at input 
point #1 of input section #2 . 
Temperature on the lower surface at 
input point #2 of input section #2. 
Temperature on the lower surface at 
input point #2 of input section #2. 
Pressure on the lower surface at input 
point #2 of input section #2 . 
Pressure on the upper surface at input 
point #2 of input section #2. 
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PROGRAM CONTROL CARDS 


These cards are required to increment the loading. 


Program option card: 

None 


Card group: 


10 and 32 




Variable 


Columns 


Format 

Name 

Entry 

1st card 

in 

card group 

10 


1-4 


14 

NTISTP 

Number of load steps. 

5-8 


14 

NMECH 

Number of mechanical cycles used in METCAN 
to account for cyclic damage (used for 
fatigue analysis) . 

9-12 


14 

NTHER 

Number of thermal cycles used in METCAN to 
account for cyclic damage (used for fatigue 
analysis) . 

13-16 


14 

LINC 

Number of load increments between load steps , 
see Figure 3.17. 

17-20 


14 

MSTART 

Write a restart file after this many load 
increments . 

21-24 


14 

MITER 

Maximum number of iterations allowed for 
global convergence per load increment. 

2nd card 





1-8 


F8.4 

TOL 

Allowable global tolerance on convergence. 

1st card 

in 

card group 

32 


1-80 


8F10.4 

TISTPS 

Time at load steps. 
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Figure 3.17: Definition Of The Load Increment 
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Card Group 10 


1 

2 

3 

4 

5 

6 

1 0 

, .0 

. . .0 

0 

0 

0 

1111 

1.0 

2 10 






Card Group 32 



i 

2 

3 

4 5 6 

1 0. 

0 

0 

. .0 0 0 

0.0 




Columns 

Field Name 

Value 

Description 

1-4 

NTISTP 

1 

Number of load steps is 1. 

5-8 

NMECH 

1 

Number of mechanical cycles 
is 1 

9-12 

NTHER 

1 

Number of thermal cycles is 1. 

13-16 

LINC 

1 

Number of load increments is 1. 

17-20 

MSTART 

2 

Write a restart file after 2 load 
increments . 

21-24 

MITER 

10 

Number of allowable iterations for 
global convergence. 

1-8 

TOL 

1.0 

Tolerance on the global convergence is 
1.0. 

1-10 

TISTPS(l) 

0.0 

Time at load step #1 is 0.0. 
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3.6 BOUNDARY CONDITIONS 


In this section, the necessary program option cards and card groups used to 
describe the boundary conditions are given. To enter the boundary conditions, 
card groups 19 and 41 are used. If the skewed boundary conditions are required 
in the analysis, coordinate transformations are available. This option requires 
the program option card TRANSFORMATION and the card groups 18 and 40. 


Chapter 3 


March, 1992 


INPUT85 


BOUNDARY CONDITIONS 


HITCAN User's Manual - Version 1.0 


I 


These cards are required to enter the boundary conditions . 


Program option card: 


None 

Card group: 



19 

and 41 




Variable 


Columns 

Format 

Name 


Entry 

1st card 

in card 

group 

19 



1-4 

14 


NBC 


Number of boundary condition data sets. 

1st card 

in card 

group 

41 



1-4 

14 


IBEG 


Beginning node number. 

5-8 

14 


IEND 


Ending node number. 

9-12 

14 


INCR 


Increment at which the nodes in this series 
are incremented. 

13-16 

14 


I DOF 


Degree-of-freedom which is fixed. 


This card is repeated NBC times. 

EXAMPLE: 

Card Group 19 

1 2 3 4 5 6 

1 0 0 0 0 0 0 

1 


Card Group 41 

1 2 3 4 5 6 

1 0 0 0 0 0 0 

1 10 1 3 


Columns 

Field Name 

Value 

Description 

1-4 

NBC 

1 

One data set is chosen to describe the 
boundary conditions. 

1-4 

IBEG 

1 

The 1st node in this series of nodes 
is 1 . 

5-8 

IEND 

10 

The last node in this series of nodes 
is 10. 
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Columns 

Field Name 

Value 

9-12 

INCR 

1 

13-16 

IDOF 

3 


Description 

The increment between node numbers 
is 1. 

The 3rd degree -of- freedom is fixed. 
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COORDINATE TRANSFORMATION 

Allows the user to specify a coordinate transformation of the global coordinate 
system into a local coordinate system at specified nodes. If more than one 
rotation is applied at a node, HITCAN executes all the entered rotations 
successively. This feature can be used to obtain transformations around an 
arbitrary axis . Each subsequent rotation acts upon the last previous coordinate 
system defined at the node. Note that this option follows the right hand 
coordinate system. 

Program option card: TRAN SFORMATION 

Card group: 18 and 40 

Variable 

Columns Format Name Entry 

1st card in card group 18 

1-4 14 NTR Number of coordinate transformation data 

sets . 

1st card in card group 40 

1-4 14 IBEG First node in this set of nodes. 

5-8 14 IEND Last node in this set of nodes. 

9-12 14 INCR Increment at which nodes in this set of nodes 

are to be incremented. 

13-16 14 IDIR Global axis about which the coordinates are 

rotated. 

2nd card in card group 40 

1-10 F10.4 TRANG Angle over which the coordinate system is 

rotated. 

The cards in card group 40 are repeated NTR times. 
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Figure 3.18 : The Different Levels of Analysis Results Available In The Primary Output File 
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OUTPUT CONTROL 


Program option card: 

None 


Card group: 



42 






Variable 


Columns 


Format 

Name 

Entrv 

1st card 

in 

card 

group 

42 


1-4 


14 


NPRT(l.l) 

Initial node for nodal displacements , 
velocities, and accelerations of set 1. 

5-8 


14 


NPRT(2 , 1) 

Final node for nodal displacements , 

velocities, and accelerations of set 1. 

9-12 


14 


NPRT(1,2) 

Initial node for nodal displacements, 
velocities, and accelerations of set 2. 

13-16 


14 


NPRT(2 , 2) 

Final node for nodal displacements , 
reactions, velocities, and accelerations 
of set 2. 

Etc. 





Etc. in 14 format. 

2nd card 

in 

card 

group 

42 


1-4 


14 


NPRTS(l.l) 

Initial node for ply stress output of 
set 1. 

5-8 


14 


NPRTS(2, 1) 

Final node for ply stress output of set 1. 

9-12 


14 


NPRTS(1,2) 

Initial node for ply stress output of set 2 , 

13-16 


14 


NPRTS(2,2) 

Final node for ply stress output of set 2. 

Etc. 





Etc. in 14 format. 

3rd card 

in 

card 

group 

42 


1-4 


14 


NPRTP(l.l) 

Initial node for constituent properties 
and stresses of set 1. 

5-8 


14 


NPRTP(2 , 1) 

Final node for constituent properties 
and stresses of set 1. 

9-12 


14 


NPRTP(1,2) 

Initial node for constituent properties 
and stresses of set 2. 

13-16 


14 


NPRTP(2 , 2) 

Final node for constituent properties 
and stresses of set 2. 

Etc . 





Etc. in 14 format. 


Chapter 3 


INPUT92 


March, 1992 ^ 



HITCAN User's Manual - Version 1.0 


Columns 

Field Name 

Value 

1-4 

NPRTP (1,1) 

41 

5-8 

NPRTP(2 , 1) 

41 

1-4 

NPPLY (1,1) 

1 

5-8 

NPPLY(2 , 1) 

4 

1-10 

TIMEPN 

10.0 


Description 

First node in the 1st series of nodal 
constituent properties and stresses is 
41. 

Last node in the 1st series of nodal 
constituent properties and stresses is 
41. 

First ply in the 1st series of ply 
constituent properties and stresses 
is 1. 

Last ply in the 1st series of ply 
constituent properties and stresses 
is 4. 

A PATRAN results file is desired at 
10.0 sec. 
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3.8 PROGRAM OPTION CARDS 

There are 29 program option cards. The option cards either control the flow 
through the program or activate various card groups. These cards are listed and 
summarized below in alphabetical order. 

ACCE LERATION 

ANGU LAR 

BRICK 

BUCK LE 

DAMP 

DISP LACEMENT 

DYNAM IC 

ECHO 

ENDO PTION 

FE M ODEL ONLY 

FORCE 

HPLA TE 

INTER PHASE 

MODA L 

PLAT E 

PLYOR DER 

PRES SURE 

PROF ILE 

READ IN MODEL 

RESTA RT 
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S3DSO LID 
S PLA TE 
STRA IN 
STRE SS 
TEMP ERATURE 
TITL E 

TRAN SFORMATION 
UNSYMMETRICAL PLYORDER 
VELOCITY 


The HPLATE, S3DSOLID, and the S PLATE option cards determine the type of finite 
element model to be generated by HITCAN. The READ IN MODEL option enables the 
user to input into HITCAN a finite element model consisting of eight node solid 
elements. One of these four cards must be included in the "Program Option Cards" 
block of the input deck. 

The ENDOPTION card is required to designate the end of the "Program Option 
Cards". Note that only the first four characters, underlined in the list of the 
option cards above , are required to be input . A brief description of each program 
option card follows. 

ACCEL ERATION 

This option enables the user to specify initial acceleration, at user- selected 
nodal points, for the direct time integration analysis. If this option is not 
specified, the initial accelerations are set to 0. 

ANGU LAR 

This option defines a centrifugal distributed loading. The variables required 
for this option can be found in card group number 35. 

BRICK 

When this option is specified, a 8 -node isoparametric 3D solid element will be 
used. This option card is to be used in conjunction with the S3DS0LID and the 
READ IN MODEL option cards. This element has 3 translational degrees of freedom 
identified by u*, Uy, and Uj, as shown in Figure 3.11. 
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4th card 

in card group 

42 


1-4 

14 

NPPLY1 , 1) 

Initial ply for constituent properties 
and stresses of set 1. 

5-8 

14 

NPPLY (2,1) 

Final ply for constituent properties 
and stresses of set 1. 

9-12 

14 

NPPLY (1,2) 

Initial ply for constituent properties 
and stresses of set 2. 

13-16 

14 

NPPLY (2,2) 

Final ply for constituent properties 
and stresses of set 2. 

Etc . 



Etc. in 14 format. 

5th card 

in card group 

42 


1-10 

F10.6 

TIMEPN(l) 

1st time at which a patran results 
file is written. 

11-20 

F10.6 

TIMEPN ( 2 ) 

2nd time at which a patran results 
file is written. 

Etc . 



Etc. in F10.6 format. 


EXAMPLE: 






Card Group 42 







1 

2 

3 

4 

5 

6 

i 

0 

. .0 

0 

0 

0 

0 

i 

81 






i 

9 73 81 






41 

41 






1 

4 







10.0 







Columns 

Field Name 

Value 

Description 




1-4 

NPRT (1,1) 

i 

First node in the 
displacements is 

1st 

1. 

series 

of nodal 

5-8 

NPRT (2,1) 

81 

Last node in the 
displacements is 

1st 

81. 

series 

of nodal 

1-4 

NPRTS(1, 1) 

1 

First node in the 
ply stress is 1. 

1st 

series 

of nodal 

5-8 

NPRTS (2 , 1) 

9 

Last node in the 
ply stress is 9. 

1st 

series 

of nodal 

9-12 

NPRTS (1 , 2) 

73 

First node in the 
ply stress is 73. 

2nd series 

of nodal 

13-16 

NPRTS (2, 2) 

81 

Last node in the 
ply stress is 81. 

2nd 

series 

of nodal 
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BUCKL E 

This option activates the buckling analysis. To determine the critical buckling 
load, MHOST uses the subspace iteration method. For a description of this method, 
see Reference 7. When this option is used, a buckling analysis is performed at 
the initial load and the times specified in the array TIMEMB. 

DAMP 


This option defines damping for the direct time integration analysis. At the 
present time, only Rayleigh damping is available in HITCAN. Thus damping matrix 
[C] is of the form 


[Cl - a[M] + B[K], 

where a is the damping coefficient of the mass matrix [M] and £ is the damping 
coefficient of the stiffness matrix (K] . 

DISPLACEMENT 


This option enables the user to specify initial displacements, at user-selected 
nodal points, for the direct time integration analysis. If this option is not 
specified, the initial displacements are set to 0. 

DYNAMIC 


This option activates the direct time integration dynamic analysis. The 
integration scheme used is the Newmark-Beta method. For a description of this 
method, see Reference 7. 

ECHO 

This option activates the print out of a reflection of the input data, prior to 
interpretation by HITCAN. No other cards are required to be input in conjunction 
with this option. 

ENDOPTION 


This option represents the end of the "Program Option Cards" block. It must be 
the last card in this block. No other cards are required to be input in 
conjunction with this option. 

FE MO DEL ONLY 

This option enables the user to generate a finite element mesh of the structure 
without conducting any analysis. The mesh is output in a PATRAN neutral file 
format. No other cards are required to be input in conjunction with this option. 
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FORC E 

This option allows the user to input static concentrated nodal loads. 

HPLA TE 

This option activates the modeling of a hollow structure using the plate element. 
This option can only be used with the plate element, thus the program option card 
PLATE must also be specified. 


INTERPHASE 


This option enables the user to specify the existence of a discrete interphase 
between the fiber and the matrix. The thickness of the interphase is specified 
as a fraction of the fiber diameter. 

MODA L 

When this option is used HITCAN will perform a free vibration dynamic analysis 
to determine the frequencies and the mode shapes . To determine the frequencies 
and the mode shapes, MHOST uses the subspace iteration method. For a description 
of this method, see Reference 7. 

PANE L 

If this option is specified along with the HPLATE program option card a finite 
element model of a hollow panel can be generated. With this option the vertical 
members of the panel are evenly spaced. Note that the spars represent the sides 
of the panel. A typical panel is illustrated in Figure 3.7. No other cards are 
required to be input in conjunction with this option. 

PLAT E 

This option selects the four-node isoparametic plate element both for the solid 
(SPLATE model option) and the hollow (HPLATE model option) structures. This 
element, which was derived from the Reissner-Mindlin theory for plates and 
shells, has 6 degrees -of- freedom at each node. These degrees -of -freedom are 
identified by u^, Uy, u z , H*, Hy, and H z . This element is shown in Figure 3.12. 

PLYORDER 


This option allows the user to build a composite model by defined plies over the 
surface of the structure and through the thickness producing an integrated 
laminated model of the entire structure. The user has the flexibility of 
selecting the number of plies, the ply thickness, the fiber volume ratio, the 
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void volume ratio, and the fiber orientation in building either a symmetric or 
an unsymmetric ply layup. When PLYORDER is used in conjunction with the program 
option card, UNSYMMETRICAL, the user can generate an unsymmetric ply layup. 
However, when the HPLATE model option is specified, only symmetric ply layups can 
be used. 

PRES SURE 

This option allows the user to input the pressure loading. For either solid or 
hollow structures using the plate element, the user can specify both a lateral 
pressure and an uniform edge pressure. For the plane stress, plane strain, and 
3D solid elements only an uniform edge pressure can be specified. Note that if 
PL and PU are used to input the pressure loading the variables TL and TU must 
also be specified. If temperature effects are not desired, then TL and TU must 
be set to the reference temperature. 


PROF ILE 

Without the presence of this card the nodal input points are assumed to be on the 
outer surface of the shell. HITCAN will then correct the grid point positions 
by moving them to the mid-wall poistion in a direction normal to the surface. 

This option will suppress this mid-wall correction and will retain the grid 
points on the external profile of the shell. This option can be used only with 
the HPLATE model option. 

READ IN MODEL 

This option allows the user to input a finite element model using eight-node 
solid elements . The user must provide the element connectivities and the nodal 
coordinates. When this option is chosen, the user must also specify the BRICK 
program option card. 

RESTA RT 

This option enables the user to conduct an analysis in several runs, an useful 
option for large problems. When the RESTART option is specified, a restart file 
created in a previous run is input. This file contains the necessary information 
to continue the analysis, including, the load step number, increment number , ply 
stresses, microstresses, microstress rates, material failure flags, and nodal 
displacements. For a restart run, the input is the same as before, except that 
RESTART is now specified in the program option cards. 

A restart file is created by specifying the variable MSTART. MSTART is the number 
of increments to be preformed in a particular run. After the analysis has 
progressed through MSTART increments, a restart file is created and the run is 
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terminated. Note that a restart file will also be created when the maximum 
n umb er of allowable iterations ( the variable MITER ) in HITCAN is exceeded. 

S3DS 0LID 

This option activates the modeling of a solid structure using the 3D solid 
element. 

SPLAT E 

This option activates the modeling of a solid structure using the plate, plane 
stress, or plane strain elements. 

STRA IN 

This option specifies, that a four-node isoparametic plane strain element will 
be used in the analysis. Each node within this element has 2 degrees -of- freedom, 
identified by u* and Uy. This element is shown in Figure 3.12. 

STRE SS 

This option specifies, that a four-node isoparametic plane stress element will 
be used in the analysis. Each node within this element has 2 degrees-of -freedom, 
identified by u* and Uy. This element is shown in Figure 3.13. 

TEMP ERATURE 

This option allows the user to input a distributed temperature loading. Note 
that if TL and TU are used to input the temperature loading the variables PL and 
PU must also be specified. If there is no lateral pressure loading, then PL and 
PU should be set 0. The ply temperatures generated by TL and TU will be 
overridden by the ply temperatures contained in the variable TEPLY. 

TITLE - 

A maximum of 5 title lines are allowed. Each title line may have a maximum of 
74 characters following the equal sign. 

TRAN SFORMATION 

This option allows the user to specify a coordinate transformation of the global 
coordinate system into a local coordinate system at specified nodes. If more 
than one rotation is applied at a node, HITCAN executes all the entered rotations 
successively. This feature can be used to obtain transformations around an 
arbitrary axis. Each subsequent rotation acts upon the last previous coordinate 
system defined at the node. Note that this option follows the right hand 
coordinate system. 
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UNSYMMETRICAL PLYORDER 


This option specifies an unsymmetric ply layup. This option is to be used in 
conjunction with the PLYORDER program option card. 

VELOCITY 


This option enables the user to specify initial velocities, at user-selected 
nodal points, for the direct time integration analysis. If this option is not 
specified, the initial velocities are set to 0. 
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3.9 DATABANK 

As discussed previously, a HITCAN input file contains all of the necessary input 
to execute HITCAN, except for the thermal -mechanical properties of the composite 
material. These properties reside in a file labeled "Data Bank". The user may 
chose any other name for this file. The data bank file supplied with the code 
contains the properties for several aerospace composite materials. This file can 
be modified by the user for other materials using the input described below. 

The organization of the data bank is shown in Figure 3.19. As can be seen in the 
figure, the first section of the data bank contains the fiber constituent 
properties, the second section contains the matrix constituent properties, and 
the third section contains the interphase properties . The data bank file 
separates the input for the different constituent materials (fiber, matrix, 
interphase) by a material identifier code as the first line of input for that 
material. This identifier is followed by the properties data, at the end of 
which, starts a new material identifier code. The material identifier codes for 
the fiber and the matrix consists of 4 unique characters. The interphase code 
consists of 9 characters. The first 4 of the 9 characters, are the code for the 
fiber type used In the interphase. The fiber type is then followed by a slash. 
The last 4 characters of the interphase code are the code for the matrix type 
used in the Interphase. Note that the interphase code must be a combination of 
the same fiber and matrix codes that are used. 

The format of variables used in the data bank is shown on the following pages. 
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FIBER PROPERTIES 




Variable 


Columns 

Format 

Name 

Entrv 

1st card 

1-2 

12 

NFTYP 

Number of fiber data sets. 

Cards 2 to 

13 are 

repeated NFTYP 

times . 

2nd card 

1-4 

A4 

DUN1 

Name of fiber. 

3rd card 

1-6 

16 

NDUN1 

Number of fibers per bundle. 

7-6 

F10.3 

DUN 2 

Fiber diameter in inches. 

4th card 

1-5 

15 

NUMFP0 

Number of fiber properties. 

6-10 

15 

NUMFS0 

Number of fiber microstresses. 

11-15 

15 

NUMFQ0 

Number of fiber microstress rates. 

16-20 

15 

NUMFT 

Number fiber exponents. 

5th card 

1-10 

F10.3 

TEMPF 

Reference temperature in °F. 

11-20 

F10.3 

TEMPMF 

Melting temperature in °F. 

21-30 

F10.3 

DOTHF 

Limit- state value of stress rate 
lb/in 2 sec. 

31-40 

F10.3 

RHOF 

Fiber density in lb/in 3 . 

41-50 

F10.3 

EF11 

Modulus longitudinal in lb/in 2 . 

51-60 

F10.3 

EF22 

Modulus transverse in lb/in 2 . 

61-70 

F10.3 

GF12 

Shear modulus in lb/in 2 . 

71-80 

F10.3 

GF23 

Shear modulus in lb/in 2 . 

6th card 

1-10 

F10.3 

NUF12 

Poisson's ratio. 

11-20 

F10.3 

NUF23 

Poisson's ratio. 

21-30 

F10.3 

CPF 

Heat capacity in Btu/lb°F. 

31-40 

F10.3 

KF11 

Thermal conductivity longitudinal in 
Btu/hr-ft 2o F/in. 

41-50 

F10.3 

KF22 

Thermal conductivity transverse in 
Btu/hr-ft 2o F/in. 
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Variable 


Columns 

Format 

Name 

Entrx 

51-60 

F10.3 

AF11 

Thermal expansion coefficient longitudinal 
in in/in/°F. 

61-70 

F10.3 

AF22 

Thermal expansion coefficient transverse 
in in/in/°F. 

71-80 

F10.3 

SFllt 

Longitudinal strength - tension in lb/in 2 . 

7th card 

1-10 

F10.3 

SFllc 

Longitudinal strength - compression in 
lb/ in 2 . 

11-20 

F10.3 

SF22t 

Transverse strength - tension in lb/in 2 . 

21-30 

F10.3 

SF22c 

Transverse strength - compression in 
lb/in 2 . 

31-40 

F10.3 

SF12s 

Shear strength in lb/in . 

41-50 

F10.3 

SF23s 

Shear strength in lb/in 2 . 

8th card 

1-10 

F10.3 

RF11 

Reference longitudinal stress in lb/in 2 . 

11-20 

F10.3 

RF22 

Reference transverse stress in lb/in 2 . 

21-30 

F10.3 

RF12 

Reference shear stress in lb/in 2 . 

31-40 

F10.3 

RF23 

Reference shear stress in lb/in 2 . 

41-50 

F10.3 

RF13 

Reference shear stress in lb/in 2 . 

9th card 

1-10 

F10.3 

QF11 

Reference longitudinal stress rate in 
lb/in 2 sec. 

11-20 

F10.3 

QF22 

Reference transverse stress rate in 
lb/in 2 sec . 

21-30 

F10.3 

QF12 

Reference shear stress rate in 
lb/in 2 sec. 

31-40 

F10.3 

QF23 

Reference shear stress rate in 
lb/in 2 sec. 

41-50 

F10.3 

QF13 

Reference shear stress rate in 
lb/in 2 sec. 

10th card 

1-10 

F10.3 

FTVCI(l) 

Modulus - exponent on stress rate. 

11-20 

F10.3 

FTVCI(2) 

Modulus - exponent on stress. 

21-30 

F10.3 

FTVCI(3) 

Modulus - exponent on temperature. 

31-40 

F10.3 

FTVCI(4) 

Strength - exponent on stress rate. 

41-50 

F10.3 

FTVCI(5) 

Strength - exponent on stress. 
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Variable 


Columns 

Format 

Name 

Entrv 

51-60 

F10.3 

FTVCI(6) 

Strength - exponent on temperature. 

61-70 

F10.3 

FTVCI(7) 

Poisson's ratio - exponent on stress rate. 

71-80 

F10.3 

FTVCI(8) 

Poisson's ratio - exponent on stress. 

11th card 

1-10 

F10.3 

FTVCI(9) 

Poisson's ratio - exponent on temperature. 

11-20 

F10.3 

FTVCI(IO) 

Not used. 

21-30 

F10.3 

FTVCI(ll) 

Not used. 

31-40 

F10.3 

FTVCI(12) 

Not used. 

41-50 

F10.3 

FTVC1(13) 

Heat conductivity - exponent on stress rate. 

51-60 

F10.3 

FTVCI(14) 

Heat conductivity - exponent on stress. 

61-70 

F10.3 

FTVCI(15) 

Heat conductivity - exponent on temperature. 

71-80 

F10.3 

FTVCI(16) 

Thermal expansion coefficient - exponent on 
stress rate. 

12th card 

1-10 

F10.3 

FTVCI(17) 

Thermal expansion coefficient - exponent on 
stress . 

11-20 

F10.3 

FTVCI(18) 

Thermal expansion coefficient - exponent on 
temperature . 

21-30 

F10.3 

FTVCI(19) 

Not used. 

31-40 

Etc. 

F10.3 

FTVCI(20) 

Not used. 

Etc. in F10.3 format 

13th card 

1-10 

F10.3 

FTVCI(25) 

Not used. 

11-20 

Etc. 

F10.3 

FTVCI(26) 

Not used. 

Etc. in F10.3 format 
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EXAMPLE: Silicon Carbide on Aluminum 


Card Group 19 

1 2 
1 0 0 

> 

3 

.0 1 

4 5 6 7 

3 0 0 0 

1 

SICA 

1 

21 

70. 

.0056 

5 5 30 

4870. 

1000000. 

.11 

62000000. 62000000. 23800000. 23800000. 

0.3 

.3 

.29 

.75 

.75 .0000018 .0000018 500000. 

650000. 

500000. 

650000. 

300000. 

300000. 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.25 

.25 

.25 

.0 

.0 .25 .25 .25 

.25 

.25 

Columns 

.25 

Field Name 

Value 

.25 .25 .25 .25 

Description 

1-2 

NFTYP 


i 

1 fiber data set. 

1-4 

DUN1 


SICA 

The fiber material identification code 

1-6 

NDUN1 


1 

is SICA. 

Number of fibers per bundle is 1. 

7-16 

DUN 2 


.0056 

The fiber diameter is .0056 inches. 

1-5 

NUMFP0 


21 

There are 21 fiber properties. 

6-10 

NUMFS0 


5 

The number of reference fiber 

11-15 

NUMFQO 


5 

microstresses is 5 . 

The number of reference fiber 

16-20 

NUMFT 


30 

microstress rates is 5. 

The number of exponents is 30. 

1-10 

TEMPF 


70 

The reference temperature is 70 °F. 

11-20 

TEMPMF 


4870 

The melting temperature is 4870 °F. 

21-30 

DOTHF 


1000000 

The limit- state value of stress rate 

31-40 

RHOF 


.11 

is 1000000. lb/in 2 sec. 

Fiber density is .11 lb/in 3 . 

41-50 

EF11 


62000000 

The longitudinal modulus is 62000000 

51-60 

EF22 


62000000 

lb/ in 2 . 

The transverse modulus is 62000000 

61-70 

GF12 


23800000 

lb/ in 2 . 

The shear modulus is 23800000 

71-80 

GF23 


23800000 

lb/in 2 . 

The shear modulus is 23800000 

1-10 

NUF12 


.3 

lb/in 2 . 

Poisson's ratio is .3. 

11-20 

NUF23 


.3 

Poisson's ratio is .3. 

21-30 

CPF 


.29 

Heat capacity is .29 Btu/lb°F. 
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Columns 

Field Name 

Value 

31-40 

KF11 

.75 

41-50 

KF22 

.75 

51-60 

AF11 

.0000018 

61-70 

AF22 

.0000018 

71-80 

SFllt 

500000 

1-10 

SFllc 

650000 

11-20 

SF22t 

500000 

21-30 

SF22c 

650000 

31-40 

SF12s 

300000 

41-50 

SF23s 

300000 

1-10 

RF11 

.0 

11-20 

RF22 

.0 

21-30 

RF12 

.0 

31-40 

RF23 

.0 

41-50 

RF13 

.0 

1-10 

QF11 

.0 

11-20 

QF22 

.0 

21-30 

QF12 

.0 

31-40 

QF23 

.0 

41-50 

QF13 

.0 

1-10 

FTVCI(l) 

.25 

11-20 

FTVCI(2) 

.25 

21-30 

FTVCI(3) 

.25 

31-40 

FTVCI(4) 

.0 

41-50 

FTVCI(5) 

.0 

51-60 

FTVCI(6) 

.25 

61-70 

FTVCI(7) 

.25 


Description 

Thermal conductivity longitudinal is 
.75 Btu/hr-ft 2o F/in. 

Thermal conductivity transverse is .75 
Btu/hr - f t 2o F/in . 

Thermal expansion coefficient 
longitudinal is .0000018 in/in/°F. 
Thermal expansion coefficient 
transverse is. 0000018 in/in/°F. 
Longitudinal strength - tension is 
500000 lb/in 2 . 

Longitudinal strength - compression is 
650000 lb/in 2 . 

Transverse strength - tension is 
500000 lb/in 2 . 

Transverse strength - compression is 
650000 lb/in 2 . 

Shear strength is 300000 lb/in . 
Shear strength is 300000 lb/in 2 . 
Reference longitudinal stress is .0 
lb/in 2 . 

Reference transverse stress is .0 
lb/in 2 . 

Reference shear stress is .0 lb/in 2 . 
Reference shear stress is .0 lb/in 2 . 
Reference shear stress is .0 lb/in 2 . 
Reference longitudinal stress rate 
is .0 lb/in 2 sec. 

Reference transverse stress rate is 
.0 lb/in 2 sec. 

Reference shear stress rate is .0 
lb/in 2 sec . 

Reference shear stress rate is .0 
lb/in 2 sec. 

Reference shear stress rate is .0 
lb/in 2 sec . 

Modulus - exponent on stress rate is 
.25. 

Modulus - exponent on stress is .25. 
Modulus - exponent on temperature is 
.25. 

Strength - exponent on stress rate 
is .0. 

Strength - exponent on stress is .0. 
Strength - exponent on temperature is 
.25. 

Poisson's ratio - exponent on stress 
rate is . 25 . 
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71-80 

FTVCI(8) 

.25 

Poisson's ratio - exponent on stress 
is .25. 

1-10 

FTVCI (9) 

.25 

Poisson's ratio - exponent on 

temperature is .25. 

11-20 

FTVCI(IO) 

Blank 


21-30 

FTVCI (11) 

Blank 


31-40 

FTVCI (12) 

Blank 


41-50 

FTVCI (13) 

.25 

Heat conductivity - exponent on stress 
rate is .25. 

51-60 

FTVCI (14) 

.25 

Heat conductivity - exponent on stress 
is .25. 

61-70 

FTVCI (15) 

.25 

Heat conductivity - exponent on stress 
temperature is . 25 . 

71-80 

FTVCI (16) 

.25 

Thermal expansion coefficient 
exponent on stress rate is .25. 

1-10 

FTVCI (17) 

.25 

Thermal expansion coefficient 
exponent on stress is .25. 

11-20 

FTVCI (18) 

.25 

Thermal expansion coefficient - 
exponent on temperature is .25. 

21-30 

FTVCI(19) 

Blank 


31-40 

FTVCI(20) 

Blank 


41-50 

FTVCI(21) 

Blank 


51-60 

FTVCI (22) 

Blank 


61-70 

FTVCI(23) 

Blank 


71-80 

FTVCI (24) 

Blank 


1-10 

FTVCI(25) 

Blank 


11-20 

FTVCI(26) 

Blank 


21-30 

FTVCI (27) 

Blank 


31-40 

FTVCI(28) 

Blank 


41-50 

FTVCI (29) 

Blank 


51-60 

FTVCI (30) 

Blank 
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MATRIX PROPERTIES 

Variable 

Columns Format Name Entry 

1st card 


1-2 

12 

NMTYP 

Number of matrix data sets. 

Cards 2 to 

14 are 

repeated NMTYP 

times . 

2nd card 

1-4 

A4 

DUN1 

Name of matrix material. 

3rd card 

1-5 

15 

NUMMP0 

Number of matrix properties. 

6-10 

15 

NUMMS0 

Number of matrix microstresses. 

11-15 

15 

NUMMQ0 

Number of matrix microstress rates. 

16-20 

15 

NUMMT 

Number of matrix exponents. 

4th card 

1-10 

F10.3 

TEMPM , 

Reference temperature in °F. 

11-20 

F10.3 

TEMPMM 

Melting temperature in °F. 

21-30 

F10.3 

DOTHM 

Limit- state value of stress rate in 
lb/in 2 sec . 

31-40 

F10.3 

RHOM 

Matrix density in lb/in 3 . 

41-50 

F10.3 

EMU 

Modulus longitudinal in lb/in 2 . 

51-60 

F10.3 

EM22 

Modulus transverse in lb/in 2 . 

61-70 

F10.3 

GM12 

Shear modulus in lb/in 2 . 

71-80 

F10.3 

GM23 

Shear modulus in lb/in 2 . 

5th card 

1-10 

F10.3 

NUM12 

Poisson's ratio. 

11-20 

F10.3 

NUM23 

Poisson's ratio. 

21-30 

F10.3 

CPM 

Heat capacity in Btu/lb°F. 

31-40 

F10.3 

KM11 

Thermal conductivity longitudinal in 
Btu/hr - f t 2o F/in . 

41-50 

F10.3 

KM22 

Thermal conductivity transverse in 
Btu/hr-ft 2o F/in. 

51-60 - 

F10.3 

AMll 

Thermal expansion coefficient 
longitudinal in in/in/°F. 

61-70 

F10.3 

AM22 

Thermal expansion coefficient transverse 
in in/in/°F. 

71-80 

F10.3 

SMllt 

Longitudinal strength - tension in 
lb/in 2 . 
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Variable 


Columns 
6th card 

Format 

Name 

Entry 

1-10 

F10.3 

SMllc 

Longitudinal strength - compression in 
lb/ in 2 . 

11-20 

F10.3 

SM22t 

Transverse strength - tension in lb/in 2 . 

11-30 

F10.3 

SM22c 

Transverse strength - compression in 
lb/in 2 . 

31-40 

F10.3 

SMI 2 s 

Shear strength in lb/in . 

41-50 

F10.3 

SM23s 

Shear strength in lb/in 2 . 

7th card 

1-10 

F10.3 

RM11 

Reference longitudinal stress in lb/in 2 . 

11-20 

F10.3 

RM22a 

Reference transverse stress - region a 
in lb/in 2 . 

21-30 

F10.3 

RM22b 

Reference transverse stress - region b 
lb/in 2 . 

31-40 

F10.3 

RM22c 

Reference transverse stress - region c 
in lb/in 2 . 

41-50 

F10.3 

RM12a 

Reference shear stress - region a in 
lb/in 2 . 

51-60 

F10.3 

RM12b 

Reference shear stress - region b in 
lb/in 2 . 

61-70 

F10.3 

RM12c 

Reference shear stress - region c in 
lb/in 2 . 

71-80 
8th card 

F10.3 

RM23a 

Reference shear stress - region a in 
lb/in 2 . 

1-10 

F10.3 

RM23b 

Reference shear stress - region b in 
lb/in 2 . 

11-20 

F10.3 

RM23c 

Reference shear stress - region c in 
lb/in 2 . 

21-30 

F10.3 

RM13a 

Reference shear stress - region a in 
lb/ in 2 . 

31-40 

F10.3 

RM13b 

Reference shear stress - region b in 
lb/in 2 . 

41-50 

F10.3 

RM13c 

Reference shear stress - region c in 
lb/in 2 . 

9th card 

1-10 

F10.3 

QM11 

Reference longitudinal stress rate 
lb/in 2 sec. 
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I 




Variable 


Columns 

Format 

Name 

Entry 

11-20 

F10.3 

QM22a 

Reference transverse stress rate - 
region a in lb/in 2 sec. 

21-30 

F10.3 

QM22b 

Reference transverse stress rate - 
region b in lb/in 2 sec. 

31-40 

F10.3 

QM22c 

Reference transverse stress rate - 
region c in lb/in 2 sec. 

41-50 

F10.3 

QM12a 

Reference shear stress rate - region a 
in lb/in 2 sec. 

51-60 

F10.3 

QM12b 

Reference shear stress rate - region b 
in lb/in 2 sec. 

61-70 

F10.3 

QM12c 

Reference shear stress rate - region c 
in lb/in 2 sec. 

71-80 

F10.3 

QM23a 

Reference shear stress rate - region a 
in lb/in 2 sec. 

10th card 




1-10 

F10.3 

QM23b 

Reference shear stress rate - region b 
in lb/in 2 sec. 

11-20 

F10.3 

QM23c 

Reference shear stress rate - region c 
in lb/in 2 sec. 

21-30 

F10.3 

QM13a 

Reference shear stress rate - region a 
in lb/in 2 sec. 

1-40 

F10.3 

QM13b 

Reference shear stress rate - region b 
in lb/in 2 sec. 

41-50 

F10.3 

QM13c 

Reference shear stress rate - region c 
in lb/in 2 sec. 

11th card 




1-10 

F10.3 

MTVCI(l) 

Modulus - exponent on stress rate. 

11-20 

F10.3 

MTVCI(2) 

Modulus - exponent on stress. 

21-30 

F10.3 

MTVCI(3) 

Modulus - exponent on temperature. 

31-40 

F10.3 

MTVCI(4) 

Strength - exponent on stress rate. 

41-50 

F10.3 

MTVCI(5) 

Strength - exponent on stress. 

51-60 

F10.3 

MTVCI(6) 

Strength - exponent on temperature. 

61-70 

F10.3 

MTVCI(7) 

Poisson's ratio - exponent on stress rate. 

71-80 

F10.3 

MTVCI(8) 

Poisson's ratio - exponent on stress. 

12th card 




1-10 

F10.3 

MTVCI(9) 

Poisson's ratio - exponent on temperature. 

11-20 

F10.3 

MTVCI(IO) 

Not used. 

21-30 

F10.3 

MTVCI(ll) 

Not used. 

31-40 

F10.3 

MTVCI(12) 

Not used. 

41-50 

F10.3 

MTVCI(13) 

Heat conductivity - exponent on stress rate. 
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Columns 

Format 

Name 

Entrv 

51-60 

F10.3 

MTVCI(14) 

Heat conductivity - exponent on stress. 

61-70 

F10.3 

MTVCI(15) 

Heat conductivity - exponent on temperature . 

71-80 

F10.3 

MTVCI(16) 

Thermal expansion coefficient - exponent on 
stress rate. 

13th card 

1-10 

F10.3 

MTVCI(17) 

Thermal expansion coefficient - exponent on 
stress. 

11-20 

F10.3 

MTVCI(18) 

Thermal expansion coefficient - exponent on 
temperature . 

21-30 

F10.3 

MTVCI(19) 

Not used. 

31-40 

Etc. 

F10.3 

MTVCI(20) 

Not used. 

Etc. in F10.3 format 

14th card 

1-10 

F10.3 

MTVCI(25) 

Not used. 

11-20 

Etc. 

F10.3 

MTVCI(26) 

Not used. 

Etc. in F10.3 format 

EXAMPLE: 

Titanium Aluminum 



1 

1 

0 

2 

.0 

3 

.0 

4 5 

.0 0. 

0 

6 

0 

7 

1 








TI15 








21 

13 13 

30 






70. 

1800. 

1000000. 

. .172 

12300000. 

12300000. 

4659091. 

4659091. 

.32 

.32 

.12 

.39 

.39 

.0000045 

. 0000045 

130000. 

130000. 

130000. 

130000. 

91000. 

91000. 




.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 




.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 




.25 

.25 

.25 

.0 

.0 

.25 

.25 

.25 

.25 




.25 

.25 

.25 

.25 

.25 

.25 
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Columns 

Field Name 

Value 

Description 

1-2 

NMTYP 

1 

1 matrix data set. 

1-4 

DUN1 

TI15 

The matrix material identification code is 
TI15 . 

1-5 

NUMMPO 

21 

There are 21 matrix properties. 

6-10 

NUMMSO 

5 

The number of reference matrix microstresses 
is 5 . 

11-15 

NUMMQO 

5 

The number of reference matrix microstress 
rates is 5 . 

16-20 

NUMMT 

30 

The number of exponents is 30. 

1-10 

TEMPM 

70 

The reference temperature is 70 °F. 

11-20 

TEMPMM 

1800 

The melting temperature is 1800 °F. 

21-30 

DOTHM 

1000000 

The limit- state value of stress rate 
is 1000000. lb/in 2 sec. 

31-40 

RHOM 

.172 

Matrix density is .172 lb/in 3 . 

41-50 

EMU 

12300000 

The longitudinal modulus is 12300000 
lb/ in 2 . 

51-60 

EM22 

12300000 

The transverse modulus is 12300000 
lb/ in 2 . 

61-70 

GM12 

45659091 

The shear modulus is 4659091 lb/in 2 . 

71-80 

GM23 

4659091 

The shear modulus is 4659091 lb/in 2 . 

1-10 

NUM12 

.32 

Poisson's ratio is .32. 

11-20 

NUM23 

.32 

Poisson's ratio is .32. 

21-30 

CPM 

.29 

Heat capacity is .12 Btu/lb°F. 

31-40 

KM11 

.75 

Thermal conductivity longitudinal is 
.39 Btu/hr-ft 2o F/in. 

41-50 

KM22 

.75 

Thermal conductivity transverse is . 39 
Btu/hr- f t 2o F/in. 

51-60 

AM11 

.0000045 

Thermal expansion coefficient 

longitudinal is .0000045 in/in/°F. 

61-70 

AM22 

. 0000045 

Thermal expansion coefficient 
transverse is. 0000045 in/in/°F. 

71-80 

SMI It 

130000 

Longitudinal strength - tension is 
130000 lb/in 2 . 

1-10 

SMllc 

130000 

Longitudinal strength - compression 
is 130000 lb/in 2 . 

11-20 

SM22t 

130000 

Transverse strength - tension is 
130000 lb/in 2 . 

21-30 

SM22c 

130000 

Transverse strength - compression is 
130000 lb/in 2 . 

31-40 

SM12s 

91000 

Shear strength is 91000 lb/in . 

41-50 

SM23s 

91000 

Shear strength is 91000 lb/in 2 . 

1-10 

RM11 

.0 

Reference longitudinal stress is .0 
lb/ in 2 . 

11-20 

RM22a 

.0 

Reference transverse stress - region 
a is .0 lb/in 2 . 

21-30 
Chapter 3 

RM22b 

.0 

Reference transverse stress - region 
b is .0 lb/in 2 . 
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Columns 

Field Name 

Value 

Description 

31-40 

RM22c 

.0 

Reference transverse stress - region 
c is .0 lb/in 2 . 

41-50 

RM12a 

.0 

Reference shear stress - region a is 
.0 lb/in 2 . 

51-60 

RM12b 

.0 

Reference shear stress - region b is 
.0 lb/in 2 . 

61-70 

RM12c 

.0 

Reference shear stress - region c is 
.0 lb/in 2 . 

71-80 

RM23a 

.0 

Reference shear stress - region a is 
.0 lb/in 2 . 

1-10 

RM23b 

.0 

Reference shear stress - region b is 
.0 lb/in 2 . 

11-20 

RM23c 

.0 

Reference shear stress - region c is 
.0 lb/ in 2 . 

21-30 

RM13a 

.0 

Reference shear stress - region a is 
.0 lb/in 2 . 

31-40 

RM13b 

.0 

Reference shear stress - region b is 
.0 lb/in 2 . 

41-50 

RM13c 

.0 

Reference shear stress - region c is 
.0 lb/in 2 . 

1-10 

QM11 

.0 

Reference longitudinal stress rate 
is .0 lb/in 2 sec. 

11-20 

QM22a 

.0 

Reference transverse stress rate - 
region a is .0 lb/in 2 sec. 

21-30 

QM22b 

.0 

Reference transverse stress rate - 
region b is .0 lb/in 2 sec. 

31-40 

QM22c 

.0 

Reference transverse stress rate - 
region c is .0 lb/in 2 sec. 

1-50 

QM12a 

.0 

Reference shear stress rate - region 
a is .0 lb/in 2 sec. 

51-60 

QM12b 

.0 

Reference shear stress rate - region 
b is .0 lb/in 2 sec. 

61-70 

QM12c 

.0 

Reference shear stress rate - region 
c is .0 lb/in 2 sec. 

71-80 

QM23a 

.0 

Reference shear stress rate - region 
a is .0 lb/in 2 sec. 

1-10 

QM23b 

.0 

Reference shear stress rate - region 
b is .0 lb/in 2 sec. 

11-20 

QM23c 

.0 

Reference shear stress rate - region 
c is .0 lb/in 2 sec. 

21-30 

QM13a 

.0 

Reference shear stress rate - region 
a is .0 lb/in 2 sec. 

31-40 

QM13b 

.0 

Reference shear stress rate - region 
b is .0 lb/in 2 sec. 

41-50 

QM13c 

.0 

Reference shear stress rate - region 
c is .0 lb/in 2 sec. 

1-10 

MTVCI(l) 

.25 

Modulus - exponent on stress rate 
.25. 
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Field Name 

Value 

Description 

11-20 

MTVC1 (2) 

.25 

Modulus - exponent on stress is .25. 

21-30 

MTVCI(3) 

.25 

Modulus - exponent on temperature is 
.25. 

31-40 

MTVCI(4) 

.0 

Strength - exponent on stress rate 
is .0. 

41-50 

MTVCI(5) 

.0 

Strength - exponent on stress is .0. 

51-60 

MTVCI ( 6 ) 

.25 

Strength - exponent on temperature is 
.25. 

61-70 

MTVCI ( 7 ) 

.25 

Poisson's ratio - exponent on stress 
rate is . 25 . 

71-80 

MTVCI (8) 

.25 

Poisson's ratio - exponent on stress 
is .25. 

1-10 

MTVCI (9) 

.25 

Poisson's ratio - exponent on 
temperature is .25. 

11-20 

MTVCI (10) 

Blank 


21-30 

MTVCI (11) 

Blank 


31-40 

MTVCI (12) 

Blank 


41-50 

MTVCI (13) 

.25 

Heat conductivity - exponent on stress 
rate is .25. 

51-60 

MTVCI (14) 

.25 

Heat conductivity - exponent on stress 
is .25. 

61-70 

MTVCI (15) 

.25 

Heat conductivity - exponent on 

temperature is .25. 

71-80 

MTVCI (16) 

.25 

Thermal expansion coefficient - 
exponent on stress rate is .25. 

1-10 

MTVCI (17) 

.25 

Thermal expansion coefficient - 
exponent on stress is .25. 

11-20 

MTVCI(18) 

.25 

Thermal expansion coefficient - 
exponent on temperature is . 25 . 

21-30 

MTVCI (19) 

Blank 


31-40 

MTVCI (20) 

Blank 


41-50 

MTVCI (21) 

Blank 


51-60 

MTVCI(22) 

Blank 


61-70 

MTVCI (23) 

Blank 


71-80 

MTVCI (24) 

Blank 


1-10 

MTVCI (25) 

Blank 


11-20 

MTVCI(26) 

Blank 


21-30 

MTVCI (27) 

Blank 


31-40 

MTVCI(28) 

Blank 


41-50 

MTVCI (29) 

Blank 


51-60 

MTVCI (30) 

Blank 
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INTERPHASE PROPERTIES 




Variable 


Columns 
1st card 

Format 

Name 

Entrv 


1-2 

12 

NDTYP 

Number of interphase data sets. 

Cards 2 to 

15 are 

repeated NDTYP 

times . 

2nd card 




1-4 

A4 

DUN1 

Name of fiber. 

5 

IX 



6-9 

A4 

DUN2 

Name of matrix. 

3rd card 




1-5 

15 

NUMDP0 

Number of interphase properties . 

6-10 

15 

NUMDS0 

Number of interphase microstresses . 

11-15 

15 

NUMDQ0 

Number of interphase microstress rates. 

16-20 

15 

NUMDT 

Number of interphase exponents. 

4th card 




1-10 

F10.3 

TEMPD 

Reference temperature in °F. 

11-20 

F10.3 

TEMPMD 

Melting temperature in °F. 

21-30 

F10.3 

DOTHD 

Limit- state value of stress rate in 




lb/in 2 sec . 

31-40 

F10.3 

RHOD 

Interphase density in lb/in 3 . 

41-50 

F10.3 

ED11 

Modulus longitudinal in lb/in 2 . 

51-60 

F10.3 

ED22 

Modulus transverse in lb/in 2 . 

61-70 

F10.3 

GD12 

Shear modulus in lb/in 2 . 

71-80 

F10.3 

GD23 

Shear modulus in lb/in 2 . 


5th card 




1-10 

F10.3 

NUD12 

Poisson's ratio. 

11-20 

F10.3 

NUD23 

Poisson's ratio. 

21-30 

F10.3 

CPD 

Heat capacity in Btu/lb°F. 

31-40 

F10.3 

KD11 

Thermal conductivity longitudinal in 
Btu/hr- f t 2o F/in. 

41-50 

F10.3 

KD22 

Thermal conductivity transverse in 
Btu/hr- ft^F/ in. 

51-60 

F10.3 

AD11 

Thermal expansion coefficient 
longitudinal in in/in/°F. 

61-70 

F10.3 

AD22 

Thermal expansion coefficient transvi 
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I 




Variable 


Columns 

Format 

Name 

Entry 

in in/in/°F. 

71-80 

F10.3 

SDllt 

Longitudinal strength - tension in 
lb/in 2 . 

6th card 

1-10 

F10.3 

SDllc 

Longitudinal strength - compression in 
lb/ in 2 . 

11-20 

F10.3 

SD22t 

Transverse strength - tension in lb/in 2 

21-30 

F10.3 

SD22c 

Transverse strength - compression in 
lb/in 2 . 

31-40 

F10.3 

SD12s 

Shear strength in lb/ in . 

41-50 

F10.3 

SD23s 

Shear strength in lb/in 2 . 

7th card 

1-10 

F10.3 

RD11 

Reference longitudinal stress in lb/in 2 

11-20 

F10.3 

RD22b 

Reference transverse stress - region b 
lb/in 2 . 

21-30 

F10.3 

RD22c 

Reference transverse stress - region c 
in lb/ in 2 . 

31-40 

F10.3 

RD12b 

Reference shear stress - region b in 
lb/in 2 . 

41-50 

F10.3 

RD12c 

Reference shear stress • region c in 
lb/in 2 . 

51-60 

F10.3 

RD23b 

Reference shear stress • region b in 
lb/in 2 . 

61-70 

F10.3 

RD23C 

Reference shear stress - region c in 
lb/in 2 . 

71-80 
8th card 

F10.3 

RD13b 

Reference shear stress - region b in 
lb/in 2 . 

1-10 

9th card 

F10.3 

RD13c 

Reference shear stress - region c in 
lb /in 2 . 

1-10 

F10.3 

QD11 

Reference longitudinal stress rate 
lb/in 2 sec . 

11-20 

F10.3 

QD22b 

Reference transverse stress rate - 
region b in lb/in 2 sec. 

21-30 

F10.3 

QD22c 

Reference transverse stress rate - 
region c in lb/in 2 sec. 

31-40 

F10.3 

QD12b 

Reference shear stress rate - region b 
in lb/in 2 sec. 
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Variable 


Columns 

Format 

Name 

Entrv 

41-50 

F10.3 

QD12c 

Reference shear stress rate - region c 
in lb/in 2 sec. 

51-60 

F10.3 

QD23b 

Reference shear stress rate - region b 
in lb/in 2 sec. 

61-70 

F10.3 

QD23c 

Reference shear stress rate - region c 
in lb/in 2 sec. 

71-80 

F10.3 

QM13b 

Reference shear stress rate - region b 
in lb/in 2 sec. 

10th card 




1-10 

F10.3 

QD13c 

Reference shear stress rate - region c 
in lb/in 2 sec. 

11th card 




1-10 

F10.3 

DTVCI(l) 

Modulus - exponent on stress rate. 

11-20 

F10.3 

DTVCI(2) 

Modulus - exponent on stress. 

21-30 

F10.3 

DTVCI(3) 

Modulus - exponent on temperature. 

31-40 

F10.3 

DTVCI(4) 

Strength - exponent on stress rate. 

41-50 

F10.3 

DTVCI(5) 

Strength - exponent on stress . 

51-60 

F10.3 

DTVCI ( 6 ) 

Strength - exponent on temperature . 

61-70 

F10.3 

DTVCI(7) 

Poisson's ratio - exponent on stress 




rate . 

71-80 

F10.3 

DTVCI (8) 

Poisson's ratio - exponent on stress. 

12th card 




1-10 

F10.3 

DTVCI (9) 

Poisson's ratio - exponent on temperature. 

11-20 

F10.3 

DTVCI (10) 

Not used. 

21-30 

F10.3 

DTVCI (11) 

Not used. 

31-40 

F10.3 

DTVCI (12) 

Not used. 

41-50 

F10.3 

DTVCI (13) 

Heat conductivity - exponent on stress rate. 

51-60 

F10.3 

DTVCI (14) 

Heat conductivity - exponent on stress. 

61-70 

F10.3 

DTVCI (15) 

Heat conductivity - exponent on temperature. 

71-80 

F10.3 

DTVCI (16) 

Thermal expansion coefficient - exponent on 




stress rate. 

13th card 




1-10 

F10.3 

DTVCI (17) 

Thermal expansion coefficient - exponent on 




stress. 

11-20 

F10.3 

DTVCI (18) 

Thermal expansion coefficient - exponent on 
temperature . 

21-30 

F10.3 

DTVCI (19) 

Not used. 

31-40 

F10.3 

DTVCI (20) 

Not used. 

Etc. 



Etc . in F10 . 3 format 
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Columns 

Format 

Variable 

Name 

Entry 


14th card 
1-10 

F10.3 

DTVCI(25) 

Not used. 


11-20 

F10.3 

DTVCI(26) 

Not used. 


Etc . 



Etc. in F10.3 format 


15th card 

1-10 

F10.3 

TREF 

Reference temperature 

in °F 

EXAMPLE: 

1 

SICA/TI15 

2 

Interphase 

3 

4 5 

6 

1 0 

0 . 

0 


.0 


1 

SICA/TI15 


21 

9 9 30 



70. 

3335. 

1000000. 

.141 

.31 

.31 

.205 

.57 

390000. 

315000. 

390000. 

195500. 

.0 

.0 

.0 

.0 

.0 




.0 

.0 

.0 

.0 

.0 




.5 

.5 

.5 

.0 

.5 




.5 

.5 



70. 




Columns 

Field Name 

Value 

1-2 

NDTYP 


1 

1-4 

DUN1 


SICA 

6-9 

DUN 2 


TI15 

1-5 

NUMDP0 


21 

6-10 

NUMDS0 


5 

11-15 

NUMDQ0 


5 

16-20 

NUMDT 


30 

1-10 

TEMPD 


70 

11-20 

TEMPMD 


3335 

21-30 

DOTHD 


1000000 
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37150000. 

37150000. 

.14229546. 

14229546 

.57 

195500. 

.0000032 

.0000032 

315000. 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.5 

.5 

.5 

.5 

.5 

.5 

.5 


Description 


1 interphase data set. 

The fiber material identification code 


is SICA. 

The matrix material identification code 


is TI15 . 

There are 21 interphase properties. 
The number of reference interphase 
microstresses is 5. 

The number of reference interphase 
microstress rates is 5. 

The n umb er of exponents is 30. 

The reference temperature is 70 °F. 
The melting temperature is 3335 °F. 
The limit- state value of stress rate 
is 1000000. lb/in 2 sec. 
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Columns 

Field Name 

Value 

Description 

31-40 

RHOD 

.141 

Interphase density is .172 lb/in 3 . 

41-50 

ED11 

37150000 

The longitudinal modulus is 37150000 
lb/in 2 . 

51-60 

ED22 

37150000 

The transverse modulus is 37150000 
lb/in 2 . 

61-70 

GD12 

14229546 

The shear modulus is 14229546 lb/in 2 . 

71-80 

GD23 

14229546 

The shear modulus is 14229546 lb/in 2 . 

1-10 

NUD12 

.31 

Poisson's ratio is .31. 

11-20 

NUD23 

.31 

Poisson's ratio is .31. 

21-30 

CPD 

.205 

Heat capacity is .205 Btu/lb°F. 

31-40 

KD11 

.57 

Thermal conductivity longitudinal is 
.57 B tu/hr - f t 2o F/ in . 

41-50 

KD22 

.57 

Thermal conductivity transverse is .57 
Btu/hr-ft 2o F/in. 

51-60 

AD11 

.0000032 

Thermal expansion coefficient 
longitudinal is .0000032 in/in/°F. 

61-70 

AD22 

.0000032 

Thermal expansion coefficient 
transverse is. 0000032 in/in/°F. 

71-80 

SDllt 

315000 

Longitudinal strength - tension is 
315000 lb/in 2 . 

1-10 

SDllc 

390000 

Longitudinal strength - compression 
is 390000 lb/ in 2 . 

11-20 

SD22t 

315000 

Transverse strength - tension is 
315000 lb/in 2 . 

21-30 

SD22c 

390000 

Transverse strength - compression is 
390000 lb/in 2 . 

31-40 

SD12s 

195500 

Shear strength is 195500 lb/in . 

41-50 

SD23s 

195500 

Shear strength is 195500 lb/in 2 . 

1-10 

RD11 

.0 

Reference longitudinal stress is .0 
lb/ in 2 . 

11-20 

RD22b 

.0 

Reference transverse stress - region 
b is .0 lb/in 2 . 

21-30 

RD22c 

.0 

Reference transverse stress - region 
c is .0 lb/in 2 . 

31-40 

RD12b 

.0 

Reference shear stress - region b is 
.0 lb/in 2 . 

41-50 

RD12c 

.0 

Reference shear stress - region c is 
.0 lb/in 2 . 

51-60 

RD23b 

.0 

Reference shear stress - region b is 
.0 lb/in 2 . 

61-70 

RD23c 

.0 

Reference shear stress - region c is 
.0 lb/in 2 . 

71-80 

RD13b 

.0 

Reference shear stress - region b is 
.0 lb/in 2 . 

1-10 

RD13c 

.0 

Reference shear stress - region c is 
.0 lb/in 2 . 

1-10 

QD11 

.0 

Reference longitudinal stress rate 
is .0 lb/in 2 sec. 
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Columns 

Field Name 

Value 

Description 

11-20 

QD22b 

.0 

Reference transverse stress rate - 
region b is .0 lb/in 2 sec. 

21-30 

QD22c 

.0 

Reference transverse stress rate - 
region c is .0 lb/in 2 sec. 

31-40 

QD12b 

.0 

Reference shear stress rate - region 
b is .0 lb/in 2 sec. 

41-50 

QD12c 

.0 

Reference shear stress rate - region 
c is .0 lb/in 2 sec. 

51-60 

QD23b 

.0 

Reference shear stress rate - region 
b is .0 lb/in 2 sec. 

61-70 

QD23c 

.0 

Reference shear stress rate - region 
c is .0 lb/in 2 sec. 

71-80 

QD13b 

.0 

Reference shear stress rate - region 
b is .0 lb/in 2 sec. 

1-10 

QD13c 

.0 

Reference shear stress rate - region 
c is .0 lb/in 2 sec. 

1-10 

DTVCI(l) 

.5 

Modulus - exponent on stress rate is 

11-20 

DTVCI(2) 

.5 

. J • 

Modulus - exponent on stress is .5. 

21-30 

DTVCI(3) 

.5 

Modulus - exponent on temperature is 

31-40 

DTVCI(4) 

.0 

• j . 

Strength - exponent on stress rate 
is .0. 

41-50 

DTVCI ( 5 ) 

.0 

Strength - exponent on stress is .0. 

51-60 

DTVCI ( 6 ) 

.5 

Strength - exponent on tempeature is 
s 

61-70 

DTVCI (7) 

.5 

• J 9 

Poisson's ratio - exponent on stress 
rate is . 5 . 

71-80 

DTVCI (8) 

.5 

Poisson's ratio - exponent on stress 
rate is .5. 

1-10 

DTVCI (9) 

.5 

Poisson's ratio - exponent on 
temperature is .5. 

11-20 

DTVCI (10) 

Blank 


21-30 

DTVCI (11) 

Blank 


31-40 

DTVCI (12) 

Blank 


41-50 

DTVCI (13) 

.5 

Heat conductivity - exponent on stress 
rate is . 5 . 

51-60 

DTVCI (14) 

.5 

Heat conductivity - exponent on stress 
is . 5 . 

61-70 

DTVCI (15) 

.5 

Heat conductivity - exponent on 
temperature is .5. 

71-80 

DTVCI (16) 

.5 

Thermal expansion coefficient 
exponent on stress rate is .5. 

1-10 

DTVCI (17) 

.5 

Thermal expansion coefficient - 
exponent on stress is .5. 

11-20 

DTVCI (18) 

.5 

Thermal expansion coefficient - 
exponent on temperature is .5. 

21-30 

DTVCI (19) 

Blank 
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Columns 

Field Name 

Value 

31-40 

DTVCI(20) 

Blank 

41-50 

DTVCI(21) 

Blank 

51-60 

DTVCI(22) 

Blank 

61-70 

DTVCI(23) 

Blank 

71-80 

DTVCI(24) 

Blank 

1-10 

DTVCI(25) 

Blank 

11-20 

DTVCI(26) 

Blank 

21-30 

DTVCI(27) 

Blank 

31-40 

DTVCI(28) 

Blank 

41-50 

DTVCI(29) 

Blank 

51-60 

DTVCI (30) 

Blank 

1-10 

TREF 

70. 
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CHAPTER 4 

EXECUTION PROCEDURE ON THE LeRC CRAY X-MP AND Y-MP 

HITCAN is presently being executed in the batch mode on the CRAY X-MP and the 
CRAY Y-MP here at the NASA LeRC. Figure 4.1 shows both the input and output 
files for HITCAN. The primary input file is input . data . The file data. bank 
contains the data bank, restart. in is the restart file. This file is only read 
in if the program option card RESTART is specified. The file restart. out is the 
restart file that is created by HITCAN. This file is generated by HITCAN when 
there is a lack of global convergence or when the variable MSTART is used. 
pat.neut is a PATRAN neutral file. This file contains the necessary data to view 
the finite element model using PATRAN. The file pat. disn contains the nodal 
displacements. This file can be used for post-processing on PATRAN. 

To compile and load HITCAN on the CRAY X-MP, the following nqs script can be 
used: 


# USER- userid PW- userowd 

# QSUB-r iob id 

# QSUB-1T cnutime 

# QSUB -1M memory 

# QSUB-eo 
set -x 

eft -d p /aerospace2 /userid/hitcan. f 

segldr -o /aerospace2 / userid/hitcan hitcan.o 

exit 

The variable userid is the user's id on the CRAY X-MP. The variable userowd is 
the user's password, iob id is the name of the job on the X-MP. eputime is the 
epu time limit and memory 

is the maximum memory size allowed, hitcan. f is the compiler 
input file containing the Fortran source code to be compiled. 

hitcan.o is the object file created by the compiler and hitcan is the executable 
file. Note that the all of the files are assumed to reside in the user's home 
directory. 

To execute HITCAN on the X-MP, the following script can be used: 

# USER- userid PW- userpwd 

# QSUB-r iob id 

# QSUB -IT eputime 

# QSUB -1M memory 

# QSUB-eo 
set -x 
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cp /aerospace2 / userid/restart . in fort. 68 
cp /aerospace2 / userid/data . bank fort. 70 

/aerospace2 / user id/hitcan < /aerospace / userid/input . data 

cp fort. 18 /aerospace2 / userid/pat . neut 
cp fort. 68 /aerospace2 / userid/restart . out 
cp fort. 76 /aerospace2 / userid/pat . disp 
exit 

To compile and load HITCAN on the CRAY Y-MP, the script below can be used: 


# USER -userid PW- userpwd 

# QSUB-r iob id 

# QSUB - IT cputime 

# QSUB -1M memory 

# QSUB-eo 
set -x 

eft -dp /wrk/ user id/hitcan. f 

segldr -o /wrk/ user id/hitcan /wrk/ user id/hitcan. o 
exit 

Here userid is the user's id on the Y-MP and userpwd is the user's password on 
the Y-MP. Note all of the files are assumed to reside in the user's workspace 
on the Y-MP. To execute HITCAN on the Y-MP, the following script can be used: 

# USER- userid PW -userowd 

# QSUB-r iob id 

# QSUB - IT cputime 

# QSUB -1M memory 

# QSUB-eo 
set -x 
cd $W 

mkdir tmp $$ 
cd tmp$$ 

cp /wrk/ userid/restart. in fort. 68 
cp /wrk/ user id/data. bank fort. 70 
/wrk/ userid/hitcan < /wrk / userid/inout . data 
cp fort. 18 /wrk/ user id/pat . neut 
cp fort. 68 /wrk/ userid/restart . out 
cp fort. 76 /wrk / user id /pat . disp 
cd /wrk/ user id 
rm -fr tmp$$ 
exit 


Here all of the underlined variables are the same as before. 
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CHAPTERS 


OUTPUT DESCRIPTION 


The following is a description of the HITCAN output. 

1 . Echo print of the input dataset to HITCAN prior to interpretation by 
HITCAN. This is triggered by specifying the program option card 
ECHO. 


2. HITCAN logo and version number. 

3. A list of the program option cards specified. 

4. A summary of the model and load data input. 

5. Number of words required to perform the analysis. If there are more 
words required for memory, a message is printed specifying the 
number of additional words needed. 


6. Description of the finite element model created, including the ply 
layup for each node. Also, the corresponding temperature and 
pressure at each node are listed, if program option cards TEMPERATURE 
or PRESSURE are specified. 

7. The record of execution beginning at the first load step is given for 
each load increment. By setting the variable arrays NPRT, NPRTS, 
NPRTP, and NPPLY the user can specify the output for each increment. 
NPRT(1,I) is the initial node of set I at which displacement output 
is desired, NPRT(2,I) is the final node. NPRTS is an array 
containing the nodes selected for output of the ply stresses. As in 
NPRT, NPRTS (1,1) is the initial node and NPRTS (2,1) is the final 
node. By specifying NPRTP the user can select sets of nodes at which 
ply properties will be output. NPPLY contains the plies selected for 
output. NPPLY(1,J) is the initial ply and NPPLY(2,J) is the final 
ply. For NPRT, NPRTS, NPRTP, and NPPLY a maximum of 10 sets in each 
is allowed. These variables are specified in card group 42. 

8. If a modal analysis is performed, the following information is given: 


a) The eigenvalue number and value. 

b) The corresponding frequency in both radians per time and cycles 
per time. 

c) The corresponding eigenvector, normalized so that the Lnorm is 

1 . 0 . 

d) The generalized mass associate with the normalization of the 
eigenvector. 
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9. If a buckling analysis is performed, the following information is 
given: 

a) The eigenvalue number and value. 

b) The corresponding eigenvector, normalized so that the Lnorm is 

1 . 0 . 
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CHAPTER 6 

EXAMFLEFROBLEM »1 

Example #1 demonstrates the thermal analysis of a composite shell structure using 
the HITCAN code. A curved shell (40 degree segment) has a radius of 120 in., a 
width of 120 in. , and a thickness of 0.8 in. Both straight edges are clamped and 
both curved edges are free. Initially the shell is subjected to a temperature 
gradient of 20 deg. F. After 10 sec., an external pressure load of 0.1 psi is 
applied to the top surface. The shell is made of Sic/Ti- 15- 3-3-3 composite 
material (Silicon Carbide fiber, Titanium matrix with 15Z Vanadium, 3X Aluminum, 
3X Chromium, and 3% Tin, and interphase with average properties of the fiber and 
the matrix) . The composite laminate consists of 4 plies (0/0/90/90) of equal 
thickness with 0.5 fiber volume ratio. The ply layup is such that the 0 degree 
plies are at the top and the 90 degree plies are on the bottom. The geometry, 
boundary conditions, loading, and ply layup are shown in Figure 6.1. The finite 
element model is shown in Figure 6.2. A portion of the output is shown after the 
input deck. 
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I 



Figure 6.1: Example Problem Number 1 
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Figure 6.2: Finite Element Mesh For Example Number 1 




Chapter 6 


March, 1992 


EXAMP3 



HITCAN User's Manual - Version 1.0 


INPUT DECK DESCRIPTION 


PROGRAM OPTION CARDS 


The program option cards chosen are: 

TITLE 

SPLATE 

PLATE 

PLYORDER 

UNSYMMETRICAL 

PRESSURE 

TEMPERATURE 

ENDOPTION 


CARD GROUP 1 

The number of material systems (plies) to be described in Card Group 28 of this 
input deck is 2. One material system is required for the ply with a fiber 
orientation angle of 0°, another material system is required for the ply with a 
fiber orientation angle of 90°. 

CARD GROUP 5 

Two cross-sections will be used to define the shell, this is indicated by NSECT 
- 2. Since, the model is described by the input X, Y, Z, and TB, IGRD is set to 
4. The number of nodes along the x-axis is 9, so IU is set to 9. Nine nodes are 
needed along the y-axis, so JU is equal to 9. The finite element model is to 
begin at x-0.0 in. and end at x-120.0 in. These values are the variables XBEG 
and XEND. 

CARD GROUP 6 

The number of plies for one-half of the wall thickness (MAXPLY) is set to 2. 
This means that the number of plies for the bottom half of the shell will be 2. 

CARD GROUP 7 

Since the program option card UNSYMMETRICAL was specified, the number of plies 
for the top half of the shell is required. Because the ply layup is to consist 
of 4 plies, and because MAXPLY is equal to 2, the variable LMAX is set to 2. 
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CARD GROUP 10 

In the this card group the number of load steps (NTISTP) is set to 2. Since the 
degradation of material properties due to cumulative mechanical/thermal load 
cycles is not to be included in this analysis, the number of mechanical cycles 
(NMECHC) is set to 1, and the number of thermal cycles (NTHERC) is also set to 
1. The number of load increments (LINC) between load steps is 2. Since a 
restart file is not desired, the variable MSTART is set to 20. The variable 
MITER, the maximum number of iterations allowed for global convergence, is set 
to 10. The tolerance on global convergence (TOL) is set to 1.0. 

CARD GROUP 16 

The number of pressure data sets (NPRES) is 0, since there are no edge loads in 
this analysis. 

CARD GROUP 17 

The number of temperature data sets (NTEMP) is 0, since the temperature is input 
at the input points. 

CARD GROUP 19 

The number of boundary condition data sets (NBC) is 11. 

CARD GROUP 27 

The array MSECT is set equal to (5, 5), thus each 

cross-section will be defined by 5 points. The next card is the number of 
elements between cross-section (NXSPC) , in this example it is 8. Input for the 
two cross-sections is provided in cross-section mid-plane and thickness form 
(specified by IGRD - 4). The next 5 records are the points of the first cross- 
section. The seventh record is the number of elements between the next 2 cross- 
sections. Since there is only 1 cross-section left the variable (NXSPC) is set 
to 0. The final 5 records in this card group represent the geometry of the 
second cross-section. 
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CARD GROUP 28 

The n umb er of different material systems is 2. Each material system is 
represented by two records, thus four lines are required. The first record 
defines the location of the ply on the model by percent of thickness, length, and 
width. The first two values on this record refer to the initial and final 
percent thickness of the thickest grid point of the model. The next four values 
are the initial and final percent width and initial and final percent length. 
Since both plies are to be used over the entire shell, this record has the values 
0.0, 100.0, 0.0, 100.0, 0.0, and 100.0. The second record specifies fiber type, 
matrix type, ply thickness, void volume ratio, fiber volume ratio, and 
orientation angle in degrees relative to the structural x-axis. Here both plies 
are identical, except for the fiber orientation angle. The fibers of the bottom 
plies are along the x-axis, while the top plies are transverse to the x-axis. 

CARD GROUP 29 

The ply stack-up order for the bottom half of the shell is designated by 2 plies 
starting at the bottom surface. This layup is in the array MPLY. MPLY is equal 
to (2, 2), where these values are the material system identification numbers. 

CARD GROUP 30 

The ply stack-up order for the top half of the shell is designated by 2 plies 
starting at the top surface. This layup is in the array NPLY. NPLY is equal to 
(1, 1), where these values are the material system identification numbers. 

CARD GROUP 32 

The time at each load step is given in this card group. 

CARD GROUP 33 

Since the program option cards PRESSURE and TEMPERATURE were used, this group 
must be given. In this card group the temperature and pressure are given using 
the variables TL, TU, PL, and PU. Each line consist of the variables TL, TU, PL, 
and PU. The number of lines in this card group is (no. of load steps) x (number 
of input sections) x (number of points per input section); i. e. , 2x2x5- 20. 

CARD GROUP 41 

This CARD GROUP contains the boundary condition data. Each record consist of the 
beginning node number, the ending node number, the node numbering increment, and 
the degree -of- freedom which is constrained. 
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CARD GROUP 42 

This card group controls the output. The nodal displacements are desired at 
nodes 5 through 32, 37 through 45, and 50 through 77. The ply stresses are 
desired at node 41. Also, a PATRAN results file containing nodal displacements 
is desired at the end of the analysis. 
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INPUT DECK FOR PROBLEM #1 


Note that the Card Group Nos. are not part of the input file. 

BLOCK #1 PROGRAM OPTION CARDS 

TITLE-CYLINDRICAL SHELL WITH PRESSURE AND TEMPERATURE LOADING 

S PLATE MODEL OPTION 

PLATE 

PLYORDER 

UNSYMMETRICAL 

PRESSURE 

TEMPERATURE 

ENDOPTION 

BLOCK #2 CARD GROUPS 


CARD 

GROUP 

NO. 




1 


2 

3 


1 . . . . 

0 


. .0 

. . .0. . 

1- 

2 





5- 

2 





5- 

4 

9 

9 



5- 


.0 

120. 



6- 

2 





7- 

2 





10- 

2 

1 

1 2 

20 10 


10- 


1 . 




16- 

0 





17- 

0 





19- 

10 





27- 

5 

5 




27- 

4 





27- 


0 . 

■41.68 

236.4 

.8 

27- 


0 . ■ 

■20.92 

239.1 

.8 

27- 


0 . 

0 . 

240. 

.8 

27- 


0 . 

20.92 

239.1 

.8 

27- 


0 . 

41.68 

236.4 

.8 

27- 

0 





27- 

120. ■ 

•41.68 

236.4 

.8 

27- 

120. - 

■20.92 

239.1 

.8 

27- 

120. 

0 . 

240. 

.8 

27- 

120. 

20.92 

239.1 

.8 

27- 

120. 

41.68 

236.4 

.8 

28- 


.0 

100. 

.0 

28- 

SICA 

TI15 


.2 

.0 

28- 


.0 

i 

00. 

.0 

28- 

SICA 

TI15 


.2 

.0 

Chapter 6 
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0 0 0 0 . 


100. 

.0 

100 

.5 

0 . 


100. 

.0 

100 

.5 
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CARD 

GROUP 

NO. 




1 

2 

3 

4 


1 

0 . . . 

0 

0 

0 

29- 

2 

2 




30- 

1 

1 




32- 


0 . 

10. 



33- 


70. 

90. 

0 . 

0 . 

33- 


70. 

90. 

0 . 

0 . 

33- 


70. 

90. 

0 . 

0 . 

33- 


70. 

90. 

0 . 

0 . 

33- 


70. 

90. 

0 . 

0 . 

33- 


70. 

90. 

0 . 

0 . 

33- 


70. 

90. 

0 . 

0 . 

33- 


70. 

90. 

0 . 

0 . 

33- 


70. 

90. 

0 . 

0 . 

33- 


70. 

90. 

0 . 

0 . 

33- 


70. 

90. 

0 . 

.1 

33- 


70. 

90. 

0 . 

.1 

33- 


70. 

90. 

0 . 

.1 

33- 


70. 

90. 

0 . 

.1 

33- 


70. 

90. 

0 . 

.1 

33- 


70. 

90. 

0 . 

.1 

33- 


70. 

90. 

0 . 

.1 

33- 


70. 

90. 

0 . 

.1 

33- 


70. 

90. 

0 . 

.1 

33- 


70. 

90. 

0 . 

.1 

41- 

9 

81 9 

1 



41- 

9 

81 9 

2 



41- 

9 

81 9 

3 



41- 

9 

81 9 

4 



41- 

9 

81 9 

5 



41- 

1 

73 9 

1 



41- 

1 

73 9 

2 



41- 

1 

73 9 

3 



41- 

1 

73 9 

4 



41- 

1 

73 9 

5 



42- 

5 

32 37 

45 50 77 



42- 

41 

41 




42- 






42- 






42- 


10. 





6 

0 


7 

0 . 
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I 

HITCAN D U T P U T SUMMARY 
CYLINDRICAL SHELL WITH PRESSURE AND TQtPERATURE L CABINS 


INITIAL LOAD 


PLY STRESSES 

(in psx. units J 

IN THE MATERIAL 

COORDINATE SYSTEM 

FOR NODE 

41 


PLY NO. 

SIBL'U 

SIBL-22 

SIGL-33 

SX6L-12 

SI6L-23 

SSL* 

1 

-0. !<«<*£♦•* 

- 9 .26<»€+0* 

3.0005466 

0.1335-64 

-4.509E-49 

-4.5645 

2 

•0.X56E+4* 

-0.1SM+O* 

3.0005460 

0. 6785-65 

-4.1095-40 

-4.1215 

? 

3 .285E+43 

-0.LS1E+0* 

3.0005460 

-6 .2155-66 

-4.1215-43 

0.1095* 

- 

3. !?«♦«♦ 

-9.149E*** 

0.0005460 

3.6345-65 

-6.5665-49 

0.5095* 

DISPLACEMENTS AFTER THE INITIAL LOAD AT SELECTED NODES 




NODE NO. 

X 

r 

Z 

TMETA-X 

THETA -Y 

TWTA- 


(in.) 

(In.) 

(in.) 

(r-M.i 

(ran. j 

(rti 

S 

-4.1205-01 

-6.2155-10 

0.7355-01 

3.7295-11 

-6.2325-02 

0.99451 

6 

-0.114E-01 

3.2205-62 

0.6615-41 

-6.1475-62 

-6.2275-62 

-0.2275* 

7 

-6.961E-02 

0.3115-62 

0.4475-41 

-4.274E-62 

-4.1645-62 

-0.6535- 

a 

*4 .667E-0Z 

0.1735-62 

0.1495-41 

-6.2985-62 

-6.1315-42 

-0.1375- 

a 

-0 .5S7E-26 

-6.1225-25 

0.2395-26 

-6.5195-26 

-6.6515-27 

-0.4135- 

10 

-0. 3625-25 

3.1035-2* 

0.1635-25 

3.7905-26 

-6.2025-27 

0.1695- 

11 

-4.3425-42 

-6.1445-03 

0 .2175*41 

3.3395-62 

3.2355-03 

0.2175- 

1 

-4.6685-02 

-0.2215-02 

3.6035-01 

3.3525-02 

-4.4185-43 

0.1125- 

12 

-o.aiiE-oa 

-0.1965-02 

3.9045-41 

3.1995-42 

-4. 1Q75-42 

3 .2535- 

:* 

- 4 .347E-02 

-0.2635-10 

0.1015400 

3.5535-11 

-4.1275-42 

0.4925- 

15 

-•3 .3115-42 

3.1965-62 

3.9045-41 

-0.1995-62 

-4.1075-42 

-4.2935- 

16 

-4.6485-02 

3.2215-02 

0.6035-41 

-4.3525-02 

-4.4185-43 

-4. 1125- 

17 

-4.3425-02 

3.1*45-03 

0.2175-41 

-4.3395-42 

3.2355-63 

-6.2175- 

IS 

*4.5625-25 

-6.1035-24 

0. 1635-25 

-6.7905-26 

-4.2025-27 

-4.1695- 

19 

*0 .2235-25 

0.U7E-24 

0. 1335-25 

3 .4335-26 

-0.3525-27 

0.7265- 

23 

-4.2065-02 

3.2625-63 

0.1945-41 

3.3795-42 

-0.5365-44 

0.1065- 

21 

-4 .3325-02 

-6.1715-02 

0 .6295-41 

0.4525-62 

3.1705-43 

0.3665- 

22 

-4 .5096-42 

-6.1345-62 

0.1005440 

3.2795-42 

-4.3035-43 

0.3985— 

23 

-4.5515-02 

-6.2495-10 

0. 1155400 

3.3702-11 

-4.5135-43 

0.9465— : 

2* 

•4 .5095*42 

0.1345-62 

0.1005440 

-4.2795-42 

-4.5035-63 

-4.3985H 

25 

-4. 3325- 02 

3. 1715-62 

0.6295-41 

-0.4525-02 

3.1705-43 

-4.3665H 

26 

-4.2065-02 

-6.2625-63 

0.1945-41 

-0.3795-02 

-4.5365-46 

-0.1065H 

27 

• 4 .2235-25 

-9.1175-24 

0.1835-25 

-4.4335-26 

-4.3525-27 

-6.726EH 

23 

■J. Lfl 55 -25 

0.1185-24 

0.180E-25 

3.3095-26 

3.3625*23 

0.6375H 

19 

-4.9775-03 

3 .3265-63 

0. 1935-41 

3.3465-02 

3.400F-46 

0.4305H 

23 

-0. 1335-02 

-0.1535-02 

0.620E-41 

3.4645-62 

-0.4365-94 

o.smH 
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21 

-0.2422*02 

-0.1302*02 

0.1032000 

0.3042-02 

*0.5272*04 

0.1792-04 

22 

-0.2652*02 

•0.1342-10 

0. 1192000 

0.1552-11 

r 0.7002-04 

0.1792-11 

27 

•9.4282*33 

3 .1172*24 

0.1532*25 

0.3452-24 

-0.7982-54 

• .2882-11 

25 

0.3212-12 

3.3572*03 

0.1552*01 

0.3472*02 

0.4042-12 

-0.7832-12 

29 

-0.1562-11 

-0.1552-02 

0.6252-01 

0.4492-02 

-0.7802-12 

• .1492—11 

40 

•0.1652*11 

-0.1792*02 

0*1032000 

0.3242-02 

-0.9012-12 

-0.2572-12 

41 

0.310E-14 

3.2312*14 

0.1282000 

0.1202-14 

0.5612*14 

0.1452-11 

42 

9 .1662-11 

3.1792-02 

0.1032000 

-0*3242-02 

0.5942-12 

*0 .2562-12 

43 

0.15 62-11 

3.1552*02 

0*4252*01 

-0.4492-92 

0.7772-12 

0*1492-11 

44 

-0.3192*12 

-0.3572*03 

0*1552*01 

-0.3472*02 

-0.4012-12 

-0*7842-12 

45 

0.4252*33 

-0.1172*24 

0*1532-25 

-0.3452-26 

0.7982-54 

0*2fK-U 

50 

9 .2652*02 

9.1342-10 

0*1192000 

-0.1542-11 

0.7002-04 

0*1792-11 

51 

0.2422-02 

9.1502-02 

0 *1032000 

-0.3042-02 

0.5272-04 

0*1792-04 

52 

0.1532-02 

3.1532-02 

0.4202*01 

-0.4442-02 

0.4842-04 

0*3782-04 

53 

0.977E-93 

•0.3242*03 

0.1932-01 

•0.3442*02 

-0.4006-04 

0.4352-04 

54 

3.1052-25 

-0.1152-24 

0.1582*25 

-0.3092-26 

-0.3422*28 

0.4372-04 

55 

0.2232-25 

3 .1172*24 

0.1532*25 

0.4532-26 

0.3522-27 

-0.7242-04 

54 

9 .2062*02 

3.2422-03 

0*1942-01 

0.3792-02 

0.5842-04 

-0.1542-03 

57 

0.3522*02 

-0.1712-02 

0.4292*01 

0.4522*02 

-0.1702-03 

-0.5442-04 

55 

9*5092*02 

•0.1542*02 

0*1082000 

0.2792-02 

0.3082-03 

-0.3982-04 

59 

9 .5512*92 

3.2492-10 

0 *1152000 

-0.3782-11 

0.5132*03 

0*9442-1 2 

60 

9.5092*02 

3.1542*02 

0*1802000 

-0.2792-02 

0.3082-03 

0*3982-04 

61 

9*3522*02 

3.1712*02 

0*4292-01 

-0.4522*02 

-0.1702-03 

0.0412-04 

62 

9*2042*02 

•0.2422*03 

0*1942-01 

-0.3792-02 

0.5842-04 

0*1042-05 

63 

9*2252*25 

•0.1172*24 

0.1532-25 

-0.4532-26 

0.3522-27 

0*7242-04 

64 

0.3422*25 

9.1032*24 

0.1452-25 

0.7902*26 

0. 2022-27 

•0.1492-03 

65 

0*3422*02 

-0.1442-03 

0.2172-01 

0.3592-02 

-0.2852-03 

*0.2172*05 

64 

0.4452*02 

-0.2212*02 

9.4052*01 

0.35ZE-02 

0.4182*03 

-0.1122*03 

o7 

3.5112*02 

-0.1962*02 

0.904E-01 

0.1992-02 

0.1072-02 

-0.2932-04 

s4 

9.5472*02 

- .2452*10 

9.1012000 

-0.5512-11 

0.1272*02 

0 *4912-11 

o9 

3.5112*02 

3.1962-02 

9.9042-01 

•0.1992-02 

0.1072-02 

0.2932-04 

70 

9.4452*02 

3.2212*02 

0.4052-01 

-0.3522-02 

0.4182-03 

0.1122*03 

71 

9.5422*02 

3.1442-03 

0.2172*01 

-0.3592*02 

-0.255E-03 

0*2172-03 

72 

9. 3422-25 

-0.1032-24 

0.1452-25 

-0.7902*26 

0.2022*27 

0.1492-03 

73 

0.5572*26 

3.1222-25 

0.2392-26 

0.5192-26 

0.6512*27 

*0.4132*03 

74 

9.4472*02 

-0.1752*02 

0 *1496*01 

0.2982*02 

0.1312*02 

*0.1572*05 

75 

9.9412-02 

-0.3112*02 

0.4472-01 

0.2742*02 

0.1642-02 

*0.4452-04 

76 

9.1142*01 

-0.2202*02 

0.6612*01 

0.1472*0 2 

0.2272*02 

•0.2272*04 

77 

9 .1202*01 

3 .2152*10 

0.7352-01 

-0.7382-11 

0.2322*02 

••9942—11 

REQUIRED T3 

: LOOP TV1U HETUM 


16.911 SEC. 




DETEt. THE OISPL. IN HHOST 


2.910 SEC. 




AM4LYSZZ2 THE 

FIRST TOTE STEF 


37.424 SEC. 
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HITCAN OUTPUT SUMMARY 
CYLINDRICAL SHELL WITH PRESSURE AND TEMPERATURE LOADING 



LOAD INCREMENT NUMBED 1 
TIME STEP Z 
CYCLE NUMBER 1 
TIME S.flOBMM 


PLY STRESSES (in pss. units) 


IN THE MATERIAL C88BOZNATE SYSTEM FOR NODE 


41 


?L r MO. 

SISL-11 

SISL-22 

SXCL-33 

SICL-12 

SXCL-25 

sxa*; 

1 

•0.14*E*** 

-0.2652+04 

t.OHEHt 

3.155E-04 

-9.5002-09 

-a. saw 

z 

-O.ISAE+O* 

-0.1572+04 


3.4012-05 

-0.1092-98 

-0.1222- 

T 

a .iasE+93 

-0.LS12+O* 

0 . 0002+00 

-0.52SE-05 

-9.12X2-90 

0.109* 

4 

9.XT7E*** 

-0.1492+04 

O.0IOE4H 

3.4142-05 


uw 


HITCAN OUTPUT SUMMARY 
CYLINDRICAL SHELL WITH PRESSURE AND TEMPERATURE LOADING 



LOAD INCREMENT NUMBER I 
TIME STEP 2 
CYCLE NUMBER I 
jXME 5.0008800 

TOTAL OISfLACEWEKTS FOR SELECTED NOTES 
NODE NO. X T 

2 

THETA -X 

THETA-Y 

THETA- 


(in. J 

(in. ) 

Cin.i 

(rad.i 

<rN.J 

( 1 - 00 . 

s 

-0.1202-01 

-0. 2132-10 

0.7322-01 

3.7422-11 

-0.2522-02 

0.99SE- 

» 

-O.114E-01 

0 .2282-02 

0.4502-01 

•9.147E-02 

-0.2272-02 



-0.9412-02 

0 .2X12-02 

0.4442-01 

-0.273E-02 

-0.1442-02 

-0.4S9E— 

s 

-0.4472-02 

0.1752-02 

9. 1492-01 

-0.2982-02 

-0.1512-02 

-0.1072- 

a 

-0.5902-24 

*0.1212-25 

3.2372-26 

-0.5192-24 

-0.4512-27 

-0.414EH 

13 

-0.5432-25 

0.1852-24 

9.1472-25 

3.7092-26 

-0.2052-27 

0.149CH 

* . 

-0.5422-02 

-9. 1432-05 

3.2142-01 

3.5092-02 

3.2042-05 

0.217EH 

12 

-0.4442-02 

*9.2222-02 

3.4072-01 

3 .5512*02 

-0.4102-05 

0.1122H 

11 

-O.O11E-02 

-0. 1952-02 

3.9022-01 

3.1982-02 

-0.107E-OT 


14 

-O. *442-02 

-0.2472-10 

9.1012000 

3.5452-11 

-0.127E-02 

1.401 
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IS 

-0.0115-02 

0. 1956*02 

0.9026*01 

•0.1986*02 

-0.1075-02 

-0.2946 >4 

lb 

-0.6406*02 

0.2286-02 

0.6075*01 

-0.3516*02 

-0.4186*03 

-0.1122*05 

17 

-4. 3425*02 

0.1436*03 

0.2166*01 

-0.3896*02 

4.2846*05 

-0.2175*05 

18 

-0.3635*25 

•0.1836*24 

0.1676*25 

-0.7896*26 

-0.2036-27 

-0.1496*05 

19 

*9.2235*25 

0.1175-24 

0.1836*25 

0.4826*26 

-0.3526-27 

0.7276*04 

20 

*0.2888*02 

0 .2625*03 

0.1946*01 

0.3786*02 

•0.5856*04 

0.1846*05 

21 

*0.3822*02 

-0.1706*02 

0.6275*01 

0.4526-02 

9.1706*05 

0.8475*04 

22 

*0.5095*02 

*0.1836*02 

0.1866000 

0.2796-02 

-0.3086*05 

0.3906*04 

23 

*0.5515*02 

-0.2496*10 

0.1155000 

0.3766*11 

-0.5135-03 

0.9526*12 

2$ 

•0.5095*02 

0.1836*02 

0.1086000 

-0.2796*02 

-0.3086*05 

-0.5986 04 

25 

*0.3825*02 

0.1706*02 

0.6275*01 

•0.4526*02 

0.1706-03 

-O *8475*04 

2b 

•0.2865*02 

*0.2625*03 

0.1946*01 

-0.3786*02 

-0.5856*04 

-0.1046*05 

27 

-0.2235*25 

-0.1175*24 

0.1835-25 

-0.4826*26 

-0.3S26-27 

-0.7276*04 

28 

-0.1855*25 

0.1175-24 

0.1806*25 

0.3076*26 

9.3616*28 

0.4376*04 

29 

-0.9786*03 

0.3275-03 

0.1935*01 

0.3456*02 

0. 4006*04 

0.43X6-04 

20 

*0.1835*02 

-0.1525*02 

9.6195*01 

0.4636*02 

-0.4866*04 

0.3786*04 

21 

-0.2425*02 

-0.1795-02 

0.1035000 

0.3056*02 

-0.5286*04 

0.1796*04 

22 

•0.2656*02 

-0.1365*10 

0.1195000 

0.1586*11 

-0.7006*04 

0.1806*11 

37 

•0.4285*33 

0.1175*24 

0.1825*25 

0.3676*26 

-0.7975-34 

0. 2816-11 

28 

9.3285*12 

0.3886*03 

0.1886*01 

0.3666*02 

0.3996-12 

-0.7956*12 

39 

-0.1566*11 

-0.1555-42 

0.6246*01 

0.4686*02 

-0.7955-12 

0.1506*11 

<40 

-0.1656-11 

-0.1796*02 

0.1036000 

0.3236*02 

-0.9125-12 

*0.8446*12 

41 

9.3896*14 

0.2315-14 

0.1285000 

0.1216*14 

0.5606*14 

0.1486*11 

42 

9.1666*11 

0.1795-02 

0.1035000 

-0.3225*02 

0.9075-12 

*0.8466*12 

43 

9.1565*11 

0.1555*02 

0.6246-01 

-0.4686*02 

0.7935-12 

0.1506*11 

44 

-0.3266*12 

-0.3886*03 

0.1886-01 

-0.3666*02 

-0.3966-12 

-0.7946*12 

45 

9.4256*33 

-0.1175-24 

0.182E-25 

-0.3675*26 

0.7975-34 

0. 2815-11 

50 

9.2656*02 

0.1366-10 

0.1195000 

-0.1575*11 

0.7696*04 

0.1806*11 

51 

9.2425*02 

0.1796-02 

9.1035000 

-0.3056*02 

0 *5232-04 

0.1795*04 

52 

9.1535*02 

0.152E-02 

9.6196-01 

-0.4636*02 

0.4466*04 

0.3785*04 

53 

3.9786-03 

-0.3275-03 

9.1936-01 

-0.3455*02 

-0.4006*04 

0.4315*04 

54 

3.1056*25 

-0.1175-24 

9.1806*25 

-0.3875-26 

-0.3612-28 

0.4375*04 

55 

3.2225*25 

0.1175-24 

9.1836*25 

0.4826*26 

0.3522-27 

-0.7275-04 

5b 

3.2866*02 

0.2625*03 

9.1946*01 

0.3786*02 

0.3852*04 

-0.1046-03 

57 

9.3826*02 

-0.1705*02 

9.6275*01 

0.4526*02 

-0.1782-03 

-0.8476*04 

58 

9.5096*02 

-0.1835-02 

9.1006000 

0.2796*02 

0.3036*03 

-0.8986*04 

59 

9.5515*02 

0.2496*10 

9.1156000 

-0.3766*11 

0 .5132*03 

0.9536*12 


Cha P ter 6 March, 1992 


OF POO ft QUAU7Y 


EXAMP13 



HITCAN User's Manual - Version 1.0 


60 

3 .5096-42 

3. 1036-02 

0.1006400 

-0.2796-02 

0.3086-03 

••3986-04 

ol 

3.3826-92 

3.1706-92 

9.6276-01 

-0.6526-02 

-0.1706-93 

0*0676-04 

0 2 

3.2006*92 

-9.2026-03 

0.1946-01 

-0.3786-02 

0.5856-94 

•0066-03 

o3 

3.222E-2S 

-9.1176-24 

9.1336-25 

-0.0026-26 

0 .3526-27 

1.7276-04 

64 

3 .3036-25 

3.1036-24 

0.1676-25 

0.7896-26 

0.2036-27 

-00696-09 

65 

3 .3426*92 

-9.1436-03 

0.2106-01 

0.3096-02 

-0.280E-93 

*0.2176-09 

66 

9.0406-92 

-O.220E-02 

9.6076-01 

0. 3516-02 

0.0106-93 

-00126-43 

67 

3.8116-92 

-0.1956-02 

0. 9926*01 

0.1986-02 

0.1076-92 

-0*2906-04 

60 

9.3006-92 

0. 2676-10 

0.1016400 

-0.5 616- 11 

0.1276-92 

0.0926-13 

69 

9 .SUE-92 

0. 1956-02 

0.9926-01 

-0.1986-02 

0.1076-92 

00906-04 

70 

0.0006-92 

0.2206-02 

0.6076-01 

-0.3516-02 

0.0106-93 

10126-03 

71 

9.3026-92 

9.1436-03 

0. 2106-01 

-0.3896-02 

2806-93 

10176-03 

72 

3.3636-25 

-9.1036-24 

0.1676-25 

-0.7896-26 

O.203E-27 

1. 1696-03 

73 

3.5906-26 

9.1216-25 

0.2376-26 

0. 5196-26 

0.6516-27 

-0.6106-03 

74 

3.6676-92 

-9.1736-02 

0.1496-01 

0 .2986-02 

0.1316-92 

*9.1076-03 

-s 

3.9616-92 

-9.3116-02 

9.6066-01 

0.2736-02 

0.1646-92 

•9.6596-00 

76 

3.U4E-91 

-9.2206-02 

9.0506-01 

O .1476-02 

0 .2276-92 

•00276-00 

77 

3.120E-91 

3.2126-10 

0.7326-01 

-0.7436-11 

0.2326-92 

0.9956-11 

REQUIRED 

T: EVALUATE THIS 

LOAD INCREMENT 


71.790 SEC. 


i 



Hire 

AH OUT 

PUT SUHMA 

R Y 



CTLXHDItClL SHELL WITH PRESSURE AND TEMPERATURE LOADING 


LOAD INCREMENT NUMBER 
TINE STEP 

cycle number 

TIME 


2 

2 

1 

1O.0OOO0O* 


PLY STRESSES tin P*x. -nxtsi IN THE MATERIAL COORDINATE SYSTEM FOR MODE 41 


PLY NO. 

SIGL-I1 

SXGL-Z2 

SZ6L-33 

SZGL-12 

SXGL-23 

sat-: 

1 

-90436+96 

•9 .2666+9* 

O.8O0E4OO 

0.1336-04 

-0.5086-00 

-•.SMC 

2 

•0. 1566400 

•4.1586404 

O.000E4O0 

0.6056-05 

-0.1096-08 

-•.12ZE- 

3 

0.1666403 

•4.1S2E+94 

9. 0006499 

-0.0356-06 

-0.1225-98 

(.11* 

4 

O.195E+00 

-9.149E+96 

O.000E4M 

0 .5976-05 

-0.5666-99 

•.SMC 


Chapter 6 


March, 1991 I 


EXAMP14 


QRSGlWi- FAC 4 ?. K» 

AC t»no« QUALITY 



HITCAN User's Manual - Version 1.0 

NXTCAH OUTPUT 

CYUNBKXQU. SHELL VZTH PRESSURE AND TEMPERATURE L0AO2MB 


SUMMARY 


LOAD INCREMENT NUMBER 
TOTE STEP 
CYCLE NUWER 
TIME 


2 

2 

I 


IS. 


TOTAL DISPLACEMENTS FOR SELECTED NODES 


D £ HO. 

X 

Y 

Z 

TH5TA-X 

TH5TA-Y 

TWTA-2 


( in.J 

(in. i 

(in.J 

(rao.J 

(rstf.l 

(rai») 

5 

-0.1285-01 

•0.2185*10 

0.7305-01 

0.7565-11 

-0.2535*02 

0.9965*11 

» 

-♦.ll*E-0i 

0 .2195*02 

0.6565-01 

-0.1465*02 

-0.2275-02 

-0.2265-04 

7 

-0.9625*02 

0.5105*02 

0.4445*01 

-0.2725*02 

-0.1645*02 

-0.4595-04 

a 

-0.4675-02 

0.1785-02 

0.1465*01 

-0.2975-02 

*0 .13X5-02 

-0.1875-OS 

9 

-I.S9ZC-26 

-0.1285*25 

0.2355*26 

-0.5185-26 

*0.6505*27 

-0.4145*03 

10 

-0.3635-25 

0.1055-24 

0.1675*25 

0.7885-26 

-0.2835-27 

•0695*03 

u 

-♦•5425-02 

-0.1455-05 

0.2165*01 

0.5885-02 

0.2645-03 

•.2175-03 

12 

-0.4495-02 

-0.2205-02 

0.6065*01 

6.5505-02 

-0.4165-03 

0.2125*03 

13 

-♦. 8125-02 

-0.195E-02 

0.9005*01 

0.1985-02 

-0.1875-02 

••2995*04 

14 

-0.6645-02 

-0.2655-10 

0.1015000 

0.5735-11 

-0.1275*02 

•.4935*11 

IS 

*♦•8125-02 

0.195E-02 

0.9005*01 

-0.1985-02 

-0.1875-02 

-•.2995-04 

16 

-0.6496-02 

0.2205-02 

0.6065-01 

-0.5505-02 

-0.4165*03 

-0.1125-03 

17 

-0 .3425-02 

0.1435*03 

0.216E-01 

-0.5885*02 

0.2645-03 

-0.2175-03 

IS 

-0 .5635-25 

-0.1035*24 

0.1675*25 

-0.7885-26 

-0.2835-27 

-0.1615*03 

19 

-♦•2225-25 

0.1175-24 

0.1625*25 

0 .4605*26 

-0.3525*27 

0.7215*04 

20 

•♦ .2865-02 

0.2655*03 

0.1945*01 

0.5775-02 

-0.5835*04 

0.1I6C-03 

21 

-♦•5535-02 

-0.1705*02 

0.6265*01 

0.4515*02 

0.1785-03 

•■ 1665*04 

22 

•♦.5098*02 

-0.1855-02 

0.1005000 

0.2785*02 

-0.3885*03 

0.3985*04 

23 

-♦.5525-02 

-0.2485-10 

0.114EO00 

0.5825-11 

-0.5135*03 

0.9585*12 

2* 

♦;sm 4 ?2 

0.1855*02 

0 .1005006 

-0.2785*02 

-0.3885*03 

-0.59K-04 

25 

•♦. 5855*02 

0.1705-02 

0.6265*01 

-0.4515-0 2 

0.1785-05 

-0.8685*04 

26 

-♦.2065*02 

-0.2655*03 

0.1948*01 

-0.3775-02 

-0.5855-04 

-0.1165*08 

27 

*0.2225*25 

-0.1175-24 

0.1825*25 

-0.4605*26 

-0.5525-27 

-0.7215*04 

28 

-♦.1055*25 

0.U7E-24 

0.1805*25 

0.5865-26 

0.5615*28 

0.6375*04 

29 

-♦. 9785*03 

0.5275-03 

0.1925-01 

0.5445-02 

0.3995-04 

0.43X5*04 

20 

-0.1855*02 

-0.1525-02 

0.6185-01 

0.4625-02 

-0.4865*04 

0.3795*04 

31 

-0.2425*02 

-0.1795-02 

0.1055000 

0.3055*02 

-0.5285*04 

0.1795-04 

32 

•0.2655*02 

-0.1565-10 

0*1X95000 

0.1605*11 

-0.7815-04 

0. 1805*11 

37 

-0.4275*33 

0.U7E-24 

0.1825*25 

0.3665*26 

-0.7965*54 

• .2125-11 

38 

:. 5355*12 

0.5895-03 

0.1885-01 

0.5655-02 

0.3945*12 

*••8015*12 
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< 

$9 

156E-11 

-9.1545*42 

9.622-01 

0. 4475*92 

-4.3125*12 

USU'U 


•Uitf’U 


0. 1020 ft 

0 .3235*92 

•4.9245*12 


<*1 

O.SOSE-l** 

0.2222*14 

8.UJE+O0 

8.1225*34 

965575-14 

O.M90-U 

42 

0.166E-U 

0.1795-42 

0.102006 

-9.2235*92 

0.9195-12 

-0.6760-13 

45 

o.xsu-u 

9.1545*42 

0.622-61 

-9.4475*95 

0.3895*12 

0.1SK-11 

4* 

-9.3S3E-12 

-«.S89f-#5 

0.1S2-61 

-4.3455-42 

-4.3915*12 

>0.0660-13 

45 

0.6270-23 

-4.U75-24 

0.1S2-2 

•4. 3545*24 

9.7945*34 

9.2020-11 

51 

0.265E-0Z 

0.1252*10 

0.112000 

-4 *1595*31 

0.7815-94 

0.1000611 

51 

!.»{£•«: 

0. 1795-02 

0.102066 

-4*3055*42 

9.5285*94 

0 .010-04 

52 

I.183&K 

9*1522*02 

0.612-61 

-4*4425*42 

ft illf ■ 66 

1.0710-64 

52 

a.»ntc-«5 

•0.2272*05 

0.192-61 

-4*5445*42 

*4.3995*94 

0.600-04 

54 

0.10SE-25 

•0. 1172-24 

0.102-25 

*4*5045*24 

•4.3415*23 

0.6070-04 

55 

8.222E-2S 

9*1172-24 

0.102-2 

0*4505*24 

0.3525*27 

>6.7200-04 

55 

9.Z06E-0Z 

0*2522-05 

O.104E-01 

0.577E-92 

0.5335*94 

>6.1060-01 

57 


•0.1782*02 

0.622-61 

0.4515*92 

-4.1706-43 

-«♦ 0100-6* 

55 

9.500E-OZ 

*0.1532-02 

O.1O0E000 

9.2785*95 

9.3835*43 


59 


0.2452*10 

9 .112000 

-9. 5335*31 

0.5135*43 

I.MM>L 

65 

o.sooe-oz 

9.1552-02 

O.1OOE006 

•9.2785*92 

9*3836*43 

I4MM< 

61 

O.SUE-92 

9.1782*02 

0.622-61 

*9.4515*95 

-4.1705-43 

6 ■ >100-0* 

62 

0.2O6£-«Z 

•0.2532*05 

0.196E-61 

•4.5775*92 

9*5335*94 

1 .1660-6! 

65 

0.222E-ZS 

-0.1172-24 

0.162-2 

*4*4395*24 

0 *3525- 27 

0.7220-0J 

64 

8. SUE-25 

9.1822*24 

0.162-2 

4*7835*24 

0*2035*27 

*0.1600-0! 

65 

9.S62E-02 

•0.1452*05 

0.22-61 

4*5835*42 

*4*2345*93 

•0.200-0! 

66 

O.660C-O2 

-0.2282-02 

0.602-61 

4*5506*92 

4.4135*43 

>6.1120-0! 

67 

8.SX2E-92 

-0.1952-02 

O.90OE-OI 

0*1986*92 

0*1075-42 

>6.2060-0* 

64 

8. SUE-02 

9.2552-10 

O.1O1E0O0 

-4.5735-11 

0.1275-92 

0 .60SE-1! 

o9 

9.SX2E-92 

0.1552-02 

9 .9006-01 

-4.1985-92 

0.1075-42 

«.206C-0< 

:: 

9.660E-OZ 

0.2252-02 

0.602-01 

-4.3505*42 

0.4135*43 

0.1120-6. 

-• 

3.S42E-0Z 

0.1432-05 

0.22-01 

-4. 3835*42 

-4.2345*43 

0.200-0 

72 

8. SUE-25 

-0.1832-24 

0.162-2 

-9.7335*24 

0.2935-27 

0.1600-6 

73 

9.50ZE-26 

0.1282-25 

0.222-26 

0 .5135*24 

0.4545*27 

•0.6160-0 

74 

8.U7E-02 

-0.1722*02 

0.1UC-O1 

4.2975-42 

0.1325*42 

>0.1070-0: 

75 

9.062E-OZ 

-0. 3132-02 

0.642-01 

9.2725-42 

0.1445*42 

-6. 6800-6 

76 

9.114E-01 

-0. 2192-02 

0.652-01 

0.1445-42 

0.2275-92 

*0.2260-0 

77 

9 .1202-01 

0. 2182-10 

t.nti-n 

-4.7545*11 

0.2335-42 

0.0060-1 
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EXAMPLE PROBLEM #2 

Example #2 demonstrates the basic characteristics of a dataset required to run 
the S3DS0LID model option in the HITCAN code. This second example is of a thick 
ring subjected to centrifugal loading at an angular velocity of 20,000 rpm. The 
ring has an inside radius of 2.875 in. and an outside radius of 3.475 in. The 
ring is made of Sic/Ti-15-3-3-3 composite material (Silicon Carbide fiber, 
Titanium matrix with 15% Vanadium, 3% Aluminum , 3% Chromium, and 3% Tin, and 
interphase with average properties of the fiber and the matrix) . The composite 
laminate consists of 7 plies each with a fiber volume ratio of 0.5. All of the 
fibers have an orientation angle of 90 deg. w.r.t. the x-axis. Since this is an 
axisymmetric problem, only a small sector (3 degrees) of the ring was modeled. 
The geometry and the boundary conditions are shown in Figure 6.3. The finite 
element model is shown in Figures 6.4 and 6.5. A portion of the output is shown 
after the input deck. 
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Figure 6.3: Example Problem Number 2 
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Figure 6 . 4 : Element Numbers For Example Number 2 
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INPUT DECK DESCRIPTION 


PROGRAM OPTION CARDS 

The program option cards chosen are: 

TITLE 

S3DSOLID 

BRICK 

PLYORDER 

UNSYMMETRICAL 

TRANSFORMATION 

ANGULAR 

ENDOPTION 


CARD GROUP 1 

The number of material systems to be described in Card Group 28 of this input 
deck is 2. The two material systems are identical, except for the ply 

thicknesses. The first material system has a ply thickness of 0.05 in., the 
second has a ply thickness of 0.1 in. 

CARD GROUP 4 

Two input planes will be used to define the ring, this is indicated by NIPL - 2. 
The number of output sections (NOSC) was chosen to be 3. The number of elements 
along the y-axis is 3, so NEYY is set to 3. Six elements are needed through the 
thickness , so NETT is equal to 6 . 

CARD GROUP 6 

The number of plies for one-half of the wall thickness (MAXPLY) is set to 4. 
This means that the number of plies for the bottom half of the ring will be 4. 

CARD GROUP 7 

Since the program option card UNSYMMETRICAL was specified, the number of plies 
for the top half of the ring needs to be given. Because the ply layup is to 
consist of 7 plies the variable LMAX is set to 3. 

CARD GROUP 10 

In the this card group the number of load steps (NTISTP) is set to 1, the number 
of mechanical cycles (NMECHC) is set to 1, and the number of thermal cycles 
(NTHERC) is also set to 1. The number of load increments (LINC) is 1. Since a 
restart file is not desired, the variable MSTART is set to 20. The variable 
MITER, the maximum number of iterations allowed for global convergence, is set 
to 10. The tolerance on global convergence (TOL) is set to 1.0. 
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CARD GROUP 18 

The number of transformation data sets (NTR) is 1. Coordinate transformation is 
needed to apply skewed boundary conditions on the right edge of the sector. 

CARD GROUP 19 

The number of boundary condition data sets (NBC) is 11. 

CARD GROUP 26 

In this card group, the first card contains the number of elements between output 
sections and the x-coordinate of the output sections are specified. Between 
output sections 1 and 2, 2 elements are desired. Between sections 2 and 3, 4 
elements are needed. Finally, between sections 3 and 4, 2 elements are needed. 
The array NXSPC will be equal to (2, 4, 2). The x-coordinate of each output 
section is in the array X. This array has the values 

(.0, .1, .5, .6). The next card is the array LSECT. This variable contains the 
number of sets of input points for each input plane. For the ring, 2 sets of 
input points for each input plane was chosen. Thus LSECT is equal to (2, 2). 
The next two cards are the array MSECT. This array contains the number of input 
points in each set of input points. The first line is for the bottom input 
plane. Here each set of input points will contain 4 points. The second line is 
for the top input plane. For this input plane, each set of input points will 
again have 4 input points. The last card contains the coordinates of an input 
point. The first 4 lines is for the first set of input points of the bottom 
input plane. Each line contains the coordinates of 1 input point. The next 4 
lines are the coordinates of the input points of the second set of input points. 
The next 8 records are the input points of the top input plane. 

CARD GROUP 28 

The number of different material systems is 2. Each material system is 
represented by two records, thus four lines are required. The first record 
defines the location of the ply on the model by percent of thickness, length, and 
width. The first two values on this record refer to the initial and final 
percent thickness of the thickest grid point of the model. The next four values 
are the initial and final percent width and initial and final percent length. 
Since both plies are to be used over the entire shell, this record has the values 
0.0, 100.0, 0.0, 100.0, 0.0, and 100.0. The second record specifies fiber type, 
matrix type, ply thickness, void volume ratio, fiber volume ratio, and 
orientation angle in degrees relative to the structural x-axis. Here both plies 
are identical, except for the ply thickness. The top and bottom plies are 0.05 
in. thick, while all the other plies are 0.1 in. thick. 
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CARD GROUP 29 

The ply stack-up order for the bottom half of the ring is designated by 4 plies 
starting at the bottom surface. This layup is in the array MPLY. MPLY is equal 
to (1, 2, 2, 2). 

CARD GROUP 30 

The ply stack-up order for the top half of the ring is designated by 3 plies 
starting at the top surface. This layup is in the array NPLY. NPLY is equal to 
( 2 , 2 , 1 ). 

CARD GROUP 32 

The time at each load step is given in this card group. 

CARD GROUP 39 

This card group contains the information required for centrifugal loading. The 
first line defines the axis about which the structure is rotating. The second 
line contains the rotational velocity at each load step in revolutions per sec. 

CARD GROUP 40 

Since the TRANSFORMATION program option card was used, this group must be given. 
This card group provides for local coordinate transformation at various nodes. 
The first line consists of the variables IBEG, IEND, INCR, and IAXIS. Here IBEG 
is set to 4, IEND is set to 252, INCR is set to 4, and IAXIS (the axis about 
which the coordinates are to be transformed) is set to 1. IAXIS corresponds to 
the x-axis. The second line has the angle of rotation, which is 87.0 degrees. 

CARD GROUP 41 

This CARD GROUP contains the boundary condition data. Each record consist of the 
beginning node number, the ending node number, the node numbering increment, and 
the degree-of-freedom which is constrained. Note that for those nodes, which had 
their coordinate systems transformed, the degree-of-freedoms which are 
constrained will be in the transformed coordinate system. 

CARD GROUP 42 

This card group controls the output. Here the nodal displacements are desired at 
nodes 1 through 252. The ply stresses are desired at nodes 17, 89, 125, and 161. 
The constituent material properties and the constituent stresses of ply number 
1 are needed at nodes 53 and 197. 
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INPUT DECK FOR PROBLEM #2 

Note that the Card Group Nos. are not part of the input file 

BLOCK #1 PROGRAM OPTION CARDS 

TITLE-THICK RING SUBJECTED TO A CENTRIFUGAL LOADING 

S3DSOLID MODEL OPTION 

BRICK 

PLYORDER 

UNSYMMETRICAL 

TRANSFORMATION 

ANGULAR 

ENDOPTION 


BLOCK #2 CARD GROUPS 

CARD 

GROUP 

NO. 

1 2 3 4 5 6 7 

1 0 0 0 0 0 0 0 . 


1- 

2 





4- 

2 

3 

3 

6 


6- 

4 





7- 

3 





10- 

1 

1 

1 

1 20 

10 

10- 


1. 




18- 

1 





19- 

11 





26- 

2 

.0 




26- 

4 

.1 




26- 

2 

.5 




26- 


.6 




26- 

2 

2 




26- 

4 

4 




26- 

4 

4 




26- 


.0 


.0 

2.875 

26- 


.0 


.05 

2.875 

26- 


.0 


.1 

2.873 

26- 


.0 


.151 

2.871 

26- 


.6 


.0 

2.875 

26- 


.6 


.05 

2.875 

26- 


.6 


.1 

2.873 

26- 


.6 


.151 

2.871 
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CARD 

GROUP 

NO. 






1 

2 

3 

4 

5 

6 


1 . ... 


0 . . . 

0 

0 

0 

0 . . . 

0 

26- 



0 

.0 

3.475 




26- 



0 

.061 

3.475 




26- 



0 

.121 

3.473 




26- 



0 

.182 

3.47 




26- 



6 

.0 

3.475 




26- 



6 

.061 

3.475 




26- 



6 

.121 

3.473 




26- 



6 

.182 

3.47 




28- 



0 

100. 

.0 

100. 

.0 

100. 

28- 

SICA 

TI15 

.05 

.0 

.35 

90. 


28- 



0 

100. 

.0 

100. 

.0 

100. 

28- 

SICA 

TI15 

.1 

.0 

.35 

90. 


29- 

1 

2 

2 

2 





30- 

2 

2 

1 






32- 



0 






39- 



0 

.0 

.0 

.6 

.0 

.0 

39- 


333. 

3 






40- 

4 

252 

4 

1 





40- 


87. 






41- 

17 

233 

36 

1 





41- 

1 

217 

36 

2 





41- 

5 

221 

36 

2 





41- 

9 

225 

36 

2 





41- 

13 

229 

36 

2 





41- 

17 

233 

36 

2 





41- 

21 

237 

36 

2 





41- 

25 

241 

36 

2 





41- 

29 

245 

36 

2 





41- 

33 

249 

36 

2 





41- 

4 

252 

4 

3 





42- 

1 

252 







42- 

17 

17 

89 

89 125 125 

161 161 




42- 

53 

53 

197 

197 






42- 1 1 


42- .0 


7 

0 . 



Chapter 6 


March, 1992 


EXAMP25 


HITCAN User's Manual - Version 1.0 


HIT CAN OUTPUT SUMMARY 
THICK RING USING SRIS ELEMENT SUBJECTED TO A CENTRIFUGAL LOADING 

INITIAL LOAD 

PLY STRESSES (in p*i. unit si IN THE MATERIAL COORDINATE SYSTEM FOR NODE 17 

PLY HO. SICL-11 SXCL-22 SICL-33 SICL-12 SICL-23 SXSL 

I 0.189E+05 8.158E+92 9.O8OE+00 -8.9Z0E-8I -8.4S0E+00 0.2791 

’LY STRESSES (in psi. units! IN THE MATERIAL COORDINATE SYSTEM FOR NODE 39 

PLY NO. SICL-11 SICL-2Z SICL-33 SXCL-12 SICL-23 SXSL' 

1 9.170E+OS 8.530E+98 8.880EPO# 9.320E+00 -4.995E4-0O 0.2551 

PLY STRESSES (in psi. units) IN THE MATERIAL COORDINATE SYSTEM FOR NODE 12S I 

PLY NO. SICL-11 SXCL-22 SICL-33 SICL-12 SICL-23 SICL- 

1 0.170E+OS 0.968E-O1 9.08OEPO0 0.3SOE-01 -4.195E+00 9.2551 

PLY STRESSES (in pti. units) IN THE MATERIAL COORDINATE SYSTEM FOR NODE 161 

PLY NO. SICL-11 SISL-22 SICL-33 SICL-12 SICL-23 SXGL- 

: 3 . 160E+0S -0.21SE+01 3.8OOEP0O 9.330E+00 3.126E+O0 9.2S5E 
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Chapter 6 


HITCAN OUTPUT SUMMARY 
THICK RING USING BRICK ELEMENT SUBJECTED TO A CENTRIFUGAL LOADING 


DISPLACEMENTS AFTER THE INITIAL LOAD AT SELECTED NODES 


NODE NO. 

X 

(in.) 

1 

:.685E-04 

Z 

3.683E-04 

3 

3.6858*04 

4 

3.6648-04 

5 

5.5718-04 

6 

5.5698-04 

7 

0.5718-04 

8 

9.569E-04 

9 

3.4576-04 

10 

0.4558-04 

11 

9.4576-04 

12 

9.4558-04 

13 

3.2276*04 

14 

3.2276-04 

15 

3.2298-04 

16 

3.2278-04 

17 

3.X666-5I 

18 

-0.498E-07 

19 

3.7148-0 7 

ZO 

-0. 1628-06 

21 

-0.2278-04 

22 

-0.2298*04 

23 

-0.2278-04 

24 

-0. 2308-04 

25 

•0.4558*04 

26 

-0.4588*04 

27 

•0.4546*04 

28 

-0. 4578*04 

29 

-0. 5698-04 

30 

-0. 5738-04 


r 

Z 

(in.) 

(in.) 

0.2286-28 

0.3378-02 

0.583E-O4 

0.3378-02 

0.1188*05 

0.3368-02 

0*1766-03 

0.336E-02 

0.4578*28 

0.3378-02 

0.5838*04 

0.3378-02 

0.1188*03 

0.3378-02 

0.1768*03 

0.3368-02 

0.6868*28 

0.3378-02 

0.5838*04 

0.3378-02 

0.1186*03 

0.3378-02 

0.1768*03 

0 .3366*02 

0.9158*28 

0 .3378*02 

0.5836*04 

0.3378-02 

0.1186-05 

0.3378-02 

0.176E-03 

0.3365-02 

0.914E-2S 

0.3378-02 

0.5848*04 

0.3378-02 

0.1188*03 

0.3376-02 

0.176E-03 

0.3368-02 

0.9098*28 

0.3378-02 

0.5858-04 

0.3376-02 

0.1188*03 

0.3376-02 

0.1768*03 

0.336E-02 

0.6818*28 

0.3378-02 

0.5856*04 

^0. 3378-02 

0.1186-03 

0.3378-02 

0.176E-O3 

0 .3366*02 

0.4546*28 

0.3378-02 

0.5858-04 

0.3378-02 
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3X 

32 

33 
24 

35 

36 

37 
35 
3* 

40 

41 

42 

43 

44 

45 
*6 
47 
45 

49 

50 

51 

52 

53 

54 

55 

56 

57 
55 

59 
do 
*i 
oZ 
oS 
»4 
o5 
66 
67 
65 

60 

70 

71 

72 

73 

74 

75 

Chapter 6 


-0 .5652*04 

0.1152*03 

-0.5702*04 

0.1762*03 

-9.6542*04 

0.2272*25 

-0.6572-04 

0.5552*04 

-0.6522*04 

0.1152-03 

-0.6532*04 

0.1762-03 

0.6602*04 

0.4452*26 

0.6622-04 

0.5552*04 

0.6602*04 

0.1172*03 

0.6622*04 

0.1752-03 

0.550E-04 

0.5952*25 

0.5522*04 

0.5532-04 

0.559E-04 

0.1172*03 

0. 5522-04 

0.1752-03 

0.440E-04 

0.1342-27 

3 .4422-04 

0.5532*04 

0.4402*04 

0.1172-03 

0.4422-04 

0.1752-03 

0 .2282-04 

0.179E-27 

0.2212-04 

0.5532-04 

9.2292-04 

0.1172-03 

0.2212-04 

0*. 1752* 93 

-0.1962*31 

0.1792-27 

0.6452-07 

0.5542*04 

-0.1192-06 

0.1172-03 

0.5062-07 

0.1752-03 

-0.2222-04 

0.1792-27 

-0.219E-04 

0.5532-04 

-0.2222-04 

0.1172-03 

-0.2192-04 

0.1752-03 

-0.4432-04 

0.1342-27 

-0.4392-04 

0.5532-04 

-0.4432-04 

0.1172-03 

-0.4392*04 

0.1752-03 

-0.5532-04 

0.5962-25 

-0.5492*04 

0.5522*04 

-0.5532-04 

0.1172-03 

-0.5502-04 

0.1752*03 

-0.6632-04 

0.4452*25 

-0.6592*04 

0.5522*04 

-0.6632-04 

0.1162*03 

-0.6602-04 

0.1752-03 

9.6402-04 

0.4312-25 

0.6352-04 

0.5632-04 

0.6412-04 

0.1172-03 


0.3372-02 
0.3362-02 
0.3372-02 
0 .3372*02 
0.3362*02 
0.3362*02 
0.3352*02 
0.3352*02 
0.3342*02 
0.3342*02 
0.3352-02 
0.3352*02 
0.3342*02 
0.3342*02 
0.3352*02 
0.3352*02 
0.3342*02 
0.3342*02 
0.3352-02 
0.3352-02 
0.3352*02 

0.3342-02 

0.3352-02 { 

0.3352*02 

0.3352*02 

0.3342-02 

0.3352-02 

0.3352-02 

0.3352-02 

0.3342*02 

0.3352*02 

0.3352*02 

0.3342-02 

0.334E-02 

0.3352-02 

0.3352*02 

0.334E-02 

0.334E-02 

0.3352-02 

0.3352*02 

0.334E-02 

0.334E-02 

0.3332-02 

0.5332-02 

0.3322*02 
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76 

3.6385-04 

0.1746-03 

0.3326-02 

77 

3 .3345-04 

0.8635-28 

0.3336-02 

78 

0.5325-04 

0.5836-04 

0.3336-02 

79 

3.5345-04 

0.1175-03 

0.3326-02 

80 

3.5325-04 

0.1746-03 

0.3326-02 

81 

3.4276-04 

3.1296-27 

0.3336-02 

82 

0. 4255-04 

0.5836-04 

0.3336-02 

83 

3.4275-04 

0.1176-03 

0.3326-02 

86 

3.4255-04 

0.1745-03 

0.3326-02 

as 

0.2146-04 

0.1736-27 

0.3336-02 

86 

0.2125-04 

0.5845-04 

0.3336-02 

87 

0.2146-04 

0.1176-03 

0.3336-02 

88 

0.2125-04 

0.1746-03 

0.3326-02 

89 

0.2326-31 

0.1736-27 

0.3336-02 

90 

-0.7526-07 

0.5846-04 

0.3336-02 

91 

3.9525-07 

0.1176-05 

0.3336-02 

92 

-0.1266-06 

0.1746-03 

0.3326-02 

93 

-0.2125-04 

0.1736-27 

0.3336-02 

94 

-0.2145-04 

0.5845-04 

0.3336-02 

95 

-0.2126-04 

0.1176-03 

0.3336-02 

96 

-0.2155-04 

0.174E-03 

0.3326-02 

97 

-0.4255-04 

0.1295-27 

9.3336-02 

98 

-0.4286-04 

0.5836-04 

0.3336-02 

99 

-0.4256-04 

0.1176-03 

0.3326-02 

100 

-0.4286-04 

0.1746-03 

0.3326-02 

101 

-0.5316-04 

0.8636-28 

0.3336-02 

102 

-0.5345-04 

0.5836-04 

0.3336-02 

103 

-0.5315-04 

0.1176-03 

0.3326-02 

104 

-0.53S6-04 

0.1746-03 

0.3326-02 

105 

-0.4376-04 

3.4316-28 

0.3336-02 

106 

- 0 . 4415*04 

3 . 5836-04 

0.3336-02 

107 

-0.4385-04 

0.1176-03 

0.3326-02 

108 

-0.4415-04 

0 .1746-03 

0.3326-02 

109 

3.4156-04 

0.4165-28 

0.3316-02 

110 

0.4185-04 

0.5825-04 

0.3316-02 

1 U 

3.4165-04 

0.1155-03 

0.3305-02 

112 

3.6195-04 

0.1736-03 

0.3306-02 

113 

3.5135-04 

0.8316-28 

0.3316-02 

114 

3.5155-04 

9.5826-04 

0.3316-02 

115 

0.513E-04 

0.1166-03 

0.3306-02 
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0.516E-94 

O. 1735- 03 

0.3306-02 

:i7 

a.410E-«4 

3.1255-27 

0.3315-02 

;:s 

a.*12£-0* 

3.5435-0* 

0.3315-02 

:io 

a.*ioE-** 

0 • 1165*05 

0.3315-02 

123 

e.W2£-»4 

0.1736-05 

0.3305-02 

:zi 

0.20SE«H 

0.1665-27 

0.3315-02 

: 22 

0.20EE-9* 

0.5435-0* 

0.3315-02 

: 22 

a .209E-9* 

0.1165-03 

0.3315-02 

: 2 * 

0.207E-94 

0.1735-05 

0 *3306-02 

125 

•9.309E-31 

0.1665-27 

0.3315-02 

125 

0.*4«E~«7 

0.5435-0* 

0.3315-02 

127 

•«.XX5E-«fc 

0.1165-03 

0.3315-02 

123 

UUE*» 

0.1735-05 

0.3306-02 

125 

-9.207E-0* 

0.1665-27 

0.3315-02 

113 

-4.209E-9* 

0.5435-0* 

0.3315-02 

1ZI 

-9.207E-9* 

0.1166-03 

0.3315-02 

122 

-9.209E-9* 

0.1735-03 

0.3305-02 

122 

-a.413E-0<* 

0.1256-27 

0.3315-02 

12* 

-4..09E-0. 

0.543E-04 

0.3315-02 

125 

-4.*X3E“4* 

0.1165-03 

0.3315-02 

125 

-4.*X0E-49 

0.1735-03 

0.3305-02 

127 

-4.5X6E-4* 

0.4315-28 

0.3315-02 

12S 

-4.5X2E-4. 

0.5425-0* 

0.3315-02 

129 

-4.5X6E-44 

0.1165-03 

0.3305-02 

;*o 

-4.SX 3E-4* 

0.1736-03 

0.3305-02 

1*1 

•4.4X8E-44 

0.*166-2S 

0.3315-02 

1*2 

-4.4X3 E-44 

0.5426-0* 

0.3316-02 

1*2 

-4.419E-44 

0.1156-03 

0 .3306-02 

1 ** 

-4.414E-44 

0.1736-03 

0.3306-02 

1*5 

3.S94E-49 

3. * 005-28 

0.3295-02 

1*5 

3 .S9ZE-4* 

3.570E-0* 

0.3295-02 

1*7 

3 .5946-49 

3.11*5-03 

0.3235-02 

1*4 

4.59XE-44 

3.1725-03 

0.3245-02 

1*9 

4.495E-0* 

0.79*5-23 

0.3295-02 

121 

0.493E-44 

0.5705-0* 

0.3295-02 

121 

0 .495E-0. 

0.11*5-03 

0.3296-02 

122 

fr/TOE-** 

0.1725-03 

0.3235-02 

122 

0.396E-44 

0.1205-27 

0.3295-02 

12 * 

a.S94E-44 

0.5705-0* 

0.3295-02 

125 

t.SHE*M 

0.11*5-03 

0.3295-02 

126 

8 .3996-99 

0.1725-03 

0.3285-02 

127 

3.19BE-44 

0.1605-27 

0.3295-02 

124 

4.X97E-44 

0.5705-0* 

0.3295-02 

129 

0.X99E-44 

0 .11*5-03 

0.3295-02 

lii 

3.196E-44 

0.1725-03 

0.3235-02 
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161 

0.2292-31 

0.1602-27 

0.3292-62 

162 

-0. 7542-97 

0.5702-44 

0.3292-62 

163 

0. 1082-46 

0.1142-43 

0.3292-62 

164 

-4.1332-46 

0.1722-63 

0.3282-42 

165 

-0.1972-44 

0.1602-27 

0.3292-42 

166 

-6.1992-44 

0.5702-44 

0.3292-42 

167 

-6.1962-44 

0.1142-63 

0.3292-42 

166 

-4.1992-44 

0.1722-63 

0.3282-62 

169 

-6.3942-44 

0.1202-27 

0.3292-42 

170 

-6.3972-44 

0.5702-44 

9.3292-62 

171 

-6.3942*64 

0.1142-63 

0.3292-42 

172 

-4.3972-44 

0.1722-63 

6.3282-62 

173 

-6.4922-44 

0.7992-28 

6.3292-62 

174 

-6.4962-64 

0.5702-44 

6.3292-62 

17S 

-6.4922-44 

0.1142-63 

6.3292-42 

176 

-6.4962-44 

0.1722-63 

0.3282-62 

177 

-6.5912-64 

0.4642-25 

0.3292-42 

178 

-6.5952-44 

0.5702-44 

0.3292-42 

179 

-6.591E-64 

0.1142-43 

6.3282-62 


180 

-6.5948-64 

0.1722-43 

0.3282-02 

181 

0.5682-64 

0.3862-28 

0.327E-02 

182 

0.5702-04 

0.5682-64 

0.3272-02 

183 

0.5682-44 

9.1142-43 

0.3272-02 

184 

0.5702-64 

0.1712-63 

0.3262-02 

135 

0.4732-64 

0.7712-28 

0.3272-02 

136 

0.4752-64 

0.5692-64 

0.3272-02 

137 

0.4732-04 

0.114E-63 

0.3272-02 

138 

0.4762-04 

0.1712-43 

0.3262-02 

139 

0.3792-44 

0.1162-27 

0.3272-02 

190 

0.33 02-44 

0.5692-44 

0.3272-02 

191 

0 . 3782-64 

0.1142*63 

0.3272-02 

192 

0.3312-64 

0.171E-63 

0.3262-62 

193 

0.1392-64 

0.1542-27 

0.3282-02 

194 

0.1912-04 

0.5692-64 

0.3272-02 

195 

0.1392-44 

0.1142-03 

0.3272-02 

196 

0.1912-64 

0.1712-43 

0.3262-02 

197 

-6.23 42-31 

0.1542-27 

0.3282-62 

198 

0.6812-67 

0.5692-64 

0.3272-02 

199 

-6.1232-66 

0.1142-63 

0.3272-02 

204 

0.1282-66 

0.1712-63 

0.3262-62 

201 

-6.1912-64 

0.1542-27 

0.3282-62 

202 

-6.1392-64 

0.5692-64 

0.3272-02 

203 

-6.1512-04 

0.1142-43 

4.3272-02 

204 

-4.1382-64 

0.171E-63 

0.3262-02 

205 

-6.3312-64 

0.1162-27 

0.3272-02 
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20* 

'#.3782-04 

9.5*96-04 

0.3275-02 

207 

-0 .3815-0* 

9.11*5-03 

0.3275-02 

208 

•0.3785-04 

9.1715-03 

0.32*5-02 

20* 

•0.476E-04 

9.7715-28 

0.3275-02 

210 

-«.*73E-0* 

9.5*05-04 

0.3275-02 

211 

-0.47*5-04 

9.11*5-03 

0.3275-02 

212 

-0.4725-04 

9.1715-03 

0.32*5-02 

213 

-0.S715-04 

9.38*5-28 

0.3275-02 

21* 

-0.5*75-04 

9.5*85*04 

0.3275-02 

215 

-0. 5715-0* 

9.11*5-03 

0.3275-02 

21* 

-*.<5*78-44- 

0.1716-03 

0.32*5-02 

217 

9.S49E-04 

0.1906-28 

0.3255-02 

218 

9.5455-04 

9.5735-04 

0.32S5-02 

21* 

9.54*5-04 

0.1125-03 

0.3255-02 

220 

9.543E-04 

9.1786-03 

0.32*5-02 

221 

9.457E-04 

0.3815-28 

0.3255-02 

222 

9.4545-04 

9.5735-04 

0.3255-02 

222 

9.45*5-04 

9. 1125-03 

0 .3255-02 

22* 

9.4525-04 

0.1785-03 

0.3245-02 

225 

9.36*5-04 

9.5715-28 

0.326E-02 

226 

0.3*35-04 

9.5735-04 

0.3255-02 

227 

9.3*55-04 

0.1125-03 

0.3255-02 

228 

9.362E-04 

9.1786-03 

0.32*5-02 

22* 

9.1835-04 

9.7*26-28 

0.32*5-02 

230 

9.1815-04 

9.57*5-04 

0.32*5-02 

231 

9.1835-04 

9.1135-03 

0.3255-02 

232 

0. 1805-04 

9.1785-03 

0.324E-02 

233 

9.1*15-31 

9.7*25-28 

0.32*5-02 

23* 

-0.64*5-07 

9.5745-04 

0.3265-02 

235 

3.1115-0* 

9.1135-03 

0.32SE-0Z 

236 

-0.17*5-0* 

0.1705-03 

0.324E-02 

237 

-0.1815-04 

0.7*25-28 

0.526E-02 

228 

-0.1835-04 

0.57*5-04 

0.326E-02 

23* 

-0. 1815-04 

9.1135-03 

0.32SE-0Z 

2*0 

-0.1845-04 

9.1705-03 

0.324E-02 

2*1 

-0.3*35-04 

9.5715-23 

0.326E-02 

2*2 

-0.36*5-04 

9.5735-04 

0 .3255-02 

2*3 

-0.3*25-04 

9.1135-03 

0 .3255-02 

2*4 

-0.36*5-04 

0.1705-03 

0.32*5-02 

2*5 

-0.4545-04 

0.3815-28 

0 .3235-02 

2*6 

-0.4575-04 

0.5735-04 

0 .3255*02 

2*7 

-0.4535-04 

0 .1125-03 

0.3255-02 

2*8 

-0.4575-04 

0.1785-03 

0.32*5-02 

2*9 

-0.54*5-04 

0.1905-28 

0.325E-02 

250 

-0.5485-04 

0.5735-04 

0.3255-02 

251 

-0.5435-04 

0.1135-03 

0.3255-02 

252 

-0.5475-04 

0.1705-03 

0.32*5-02 
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EXAMPLE PROBLEM #3 


Example #3 demonstrates the basic characteristics of a dataset required to run 
the HPLATE model option in the HITCAN code. This example is of a modal analysis 
of a hollow built-up structure. The structure has 0.5 in. length, 0.2 in. width 
, 0.075 in. total thickness, 0.02 in. thickness at the top plate, 0.01 in. 
thickness of the bottom plate, and 3 spars in the x-z plane equally spaced in the 
y- direction with 0.02 in. thickness. The structure is made of Sic/Ti-15-3-3-3 
composite material (Silicon Carbide fiber, Titanium matrix with 15% Vanadium, 3% 
Aluminum, 3% Chromium, and 3% Tin, and interphase with average properties of the 
fiber and the matrix) . The composite laminate consists of 4 plies (90/0/0/90) 
of equal thickness for the top plate, 2 plies (90/90) of equal thickness for the 
bottom plate, and 4 plies (0/0/0 /0) of equal thickness for the spars. Each ply 
has a fiber volume ratio of 0.5. See Figure 6.6 for a complete description of 
the geometry, boundary conditions, and the loading. Figure 6.7 illustrates the 
ply layup. The finite element model is shown in Figures 6.8 and 6.9. A portion 
of the output is shown after a discussion of the input deck. 
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Figure 6.7: Ply Layup 
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Figure 6 . 9 : Finite Element Showing Node Numbers 
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INPUT DECK DESCRIPTION 


PROGRAM OPTION CARDS 

The program option cards chosen are: 

TITLE 

HPLATE 

PANEL 

PLATE 

PROFILE 

PLYORDER 

MODAL 

FORCE 

ENDOPTION 


CARD GROUP 1 

The number of material systems to be described in Card Group 28 of this input 
deck is 2. One material system is required for the ply with a fiber orientation 
angle of 0°, another material system is required for the ply with a fiber 
orientation angle of 90°. 

CARD GROUP 2 

The number of output sections (NOSC) was chosen to be 2. Both of the output 
sections are to have spars, so the variable NXSPAR is set to 2. The number of 
spars was chosen to be 3. Since NYSPC is equal to 7, the number of elements along 
the 

y-axis will be 8. 

CARD GROUP 6 

The number of plies for one -half of the wall thickness (MAXPLY) is set to 2. 
This means that the maximum number of plies available to fill the wall thickness 
is 4. 

CARD GROUP 9 

The first value on this line is the number of input cross sections on the top 
surface. The second value is the number of input cross sections on the bottom 
surface. Thus, the array LSECT has the values (3, 3). 


Chapter 6 


March, 1992 


EXAMP38 



HITCAN User's Manual - Version 1.0 


CARD GROUP 10 

In the this card group the number of load steps (NTISTP) is set to 1, the number 
of mechanical cycles (NMECHC) is set to 1, and the number of thermal cycles 
(NTHERC) is also set to 1. The number of load increments (LINC) between load 
steps is 1. Since a restart file is not desired, the variable MSTART is set to 
20. The variable MITER, the maximum number of iterations allowed for global 
convergence, is set to 10. The tolerance on global convergence (TOL) is set to 
1 . 0 . 

CARD GROUP 11 

Since a modal analysis is desired this card group must be used. NEIGV is set to 
1. INCREG is set to 0 card and MHITER is set to 20. 

CARD GROUP 19 

The number of boundary condition data sets (NBC) is 5. 

CARD GROUP 24 

In this card group the number of elements between output sections, the x- 
coordinate of the output sections, and the spar descriptions are given. The 
first line in this card group contains the ply designation numbers for the spars. 
This is the array NSPDES. Here this array has the values (1, 1, 1). The next 
line has the number of elements between output sections 1 and 2 and the x- 
coordinate of output section 1. The following line describes the 3 spars. Note 
that since the program option card PANEL was specified, the location of the spars 
is not needed. HITCAN will automatically place the spars at each end of the 
panel. In this example each spar has a wall thickness of 0.02 in. The array 
NXSPC will be equal to (2, 2). The x-coordinate of each output section is in the 
array X. This array has the values (.0, .25, .5). The next card is the number 
of points for each input cross section on the top surface. Each input cross 
section will have 3 input points. The next card contains this same information, 
but for the bottom of the panel. Here again 3 input points will be used for each 
input cross section. The next card contains the coordinates of an input point. 
The first 9 lines is for the input points of the top of the panel. The first 3 
lines are for the first cross section. The next 3 lines are for the second cross 
section, etc. Each line contains the coordinates and wall thickness of 1 input 
point. The next 9 lines are the coordinates of the input points of the bottom 
of the panel. 
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I 


The number of different material systems is 2. Each material system is 
represented by two records, thus four lines are required. The first record 
defines the location of the ply on the model by percent of thickness, length, and 
width. The first two values on this record refer to the initial and final 
percent thickness of the thickest grid point of the model. The next four values 
are the initial and final percent width and initial and final percent length. 
Since both plies are to be used over the entire shell, this record has the values 
0.0, 100.0, 0.0, 100.0, 0.0, and 100.0. The second record specifies fiber type, 
matrix type, ply thickness, void volume ratio, fiber volume ratio, and 
orientation angle in degrees relative to the structural x-axis. Both plies are 
identical, except for the ply orientation. 

CARD GROUP 29 

The ply stack-up order for one-half of the wall thickness is designated by 2 
plies starting at the bottom of the top surface. This layup is in the array 
MPLY. MPLY is equal to 
( 2 , 1 ). 

CARD GROUP 32 

The time at each load step is given in this card group. 

CARD GROUP 36 

This card group contains the information required for concentrated force loading. 
The first line contains the node number and the direction of the applied force. 
The second line contains the value of the force at each load step. 

CARD GROUP 41 

This CARD GROUP contains the boundary conditions. Each record consists of the 
beginning node number, the ending node number, the node numbering increment, and 
the degree-of-freedom which is constrained. 

CARD GROUP 42 

This card group controls the output. Since the output controlled by this card 
is not desired, the cards in this group are blank. 

CARD GROUP 44 

Since a buckling analysis is desired the last card must contain the time at which 
the analysis is to be done. 
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INPUT DECK FOR PROBLEM # 3 

Note that the Card Group Nos. are not part of the Input File. 

BLOCK #1 PROGRAM OPTION CARDS 

TITLE-BUCKLING ANALYSIS OF A PANEL SUBJECTED TO A NODAL FORCES 

HPLATE MODEL OPTION 

PANEL 

PLATE 

PROFILE 

PLYORDER 

MODAL 

FORCE 

ENDOPTION 


BLOCK #2 CARD GROUPS 

CARD 

GROUP 

NO. 




1 

2 

3 


1 . 

. . .0 

0 

0 

1- 

2 




2- 

2 

2 3 

7 


6- 

2 




9- 

3 

3 



10- 

1 

1 1 

1 20 10 


10- 


1 . 



11- 

1 

0 20 



15- 

10 




19- 

5 




24- 

1 

1 1 



24- 

2 

.0 



24- 


.0 

eg 

o 

.0 

24- 

2 

.25 



24- 


.0 

CM 

o 

.0 

24- 


.5 



24- 


.0 

.02 

.0 

24- 

3 

3 3 



24- 

3 

3 3 



24- 


.0 

-.1 

.04 

24- 


.0 

.0 

.04 

24- 


.0 

.1 

.04 

24- 


.25 

-.1 

.04 

24- 


.25 

.0 

.04 
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4 5 6 7 

0 0 0 0 . 


CM 

O 

.0 

CM 

O 

O 

ro 

.0 

CM 

O 

CM 

o 

.0 

CM 

O 


.02 

.02 

.02 

.02 

.02 
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( 


CARD 

GROUP 

NO. 



1 

1 0 

24- 

.25 

24- 

.5 

24- 

.5 

24- 

.5 

24- 

.0 

24- 

.0 

24- 

.0 

24- 

.25 

24- 

.25 

24- 

.25 

24- 

.5 

24- 

.5 

24- 

.5 

28- 

.0 

28- 

SICA TI15 

28- 

.0 

28- 

SICA TI15 

41- 

37 54 1 

41- 

9 81 18 

41- 

10 82 18 

41- 

9 81 18 

41- 

42- 
42- 
42- 
42- 

10 82 18 

44- 

.0 

44- 

.0 


2 

3 

0 . . . 

0 . 

.1 

.04 

-.1 

.04 

.0 

.04 

.1 

.04 

-.1 

.035 

.0 

.035 

.1 

.035 

-.1 

.035 

.0 

.035 

.1 

.035 

-.1 

.035 

.0 

.035 

.1 

.035 

100. 

.0 

.005 

.0 

100. 

.0 

.005 

.0 


1 

2 

2 

3 

3 


4 

5 

o 

0 . 

.02 


.02 


.02 


.02 


.01 


.01 


.01 


.01 


.01 


.01 


.01 


.01 


.01 


100. 

.0 

.4 

.0 

100. 

.0 

.4 

90. 


6 

0 


7 

0 . 


100 . 

100 . 
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EIGENVALUE EXTRACTION 


EIGENVALUE HUMBER I VALUE * 9.*17307E*1I 


FREWENCT IN RADIANS PER TIRE* 
IN CYCLES PER TIRE* 


E.NWEHS 

>.?CT2SEm 


EIGENVECTOR 


Chapter 6 


1 

-0.448442-04 

-0.501332*03 

9.198792-01 

9.440142-01 

9.452242-02 

-0 .143555-92 

z 

-9.443712-04 

-9.474392-93 

9.184132-01 

9. 744222-01 

0.185722-01 

-0.144245-95 

z 

- 9 . 745212-04 

-O .450132-03 

9.141272-01 

9.92X052-01 

9.348372-02 

-0.151435-92 

<♦ 


-0.422922-03 

0.137X92-01 

0. 922442-01 

-0.424252-02 

0*524185-93 

5 

-9.101702-93 

-0. 394712-03 

9.114582-01 

9.732252-01 

9.842352-03 

*0*25/772-92 

* 

-0.132772-03 

-0. 298412-03 

0.904882-02 

9.1X342 

0.105382-01 

0*734455-95 

7 

-0.135732-03 

-0.201452-03 

9.555572-02 

0.15X73 

-0.354092-02 

*9.172825-92 

5 

-0.1X9002-03 

-0.102592-03 

0.189312-02 

0.11079 

-0*183732-01 

0 *249415-92 

4 

-9.740542-04 

0.00000 

0.00000 

0. 551 152-01 

-0.184052-01 

*9.131945-92 

10 

9.422652-04 

3.00000 

0.00000 

-0.409X52-01 

-0.573072-02 

*9.278415*92 

11 

9.142702-03 

9.175282-03 

0.493242-03 

9*15274 

-0*511432-01 

*9.214335-93 

12 

9.194332-03 

9.344342*03 

9.541532-02 

9.19121 

-0*248042-02 

*9*249245-92 

12 

9.1887*2-03 

3.514072-03 

0.104002-01 

9.15454 

0.439702-01 

0 •392035-92 

14 

9.114582-03 

3.487132-03 

9 .114742-01 

-0.254352-01 

-0.340702-02 

-0*294385-92 

15 

3.114842-03 

3.729902-03 

9.120082-01 

9.973972-01 

-0.340802-01 

0. 23X595-92 

16 

3.10741E-93 

3.774072-03 

9.155942-01 

3.13321 

-0.257012-02 

-0.242322-02 

17 

3 . 942322-04 

3.819072-03 

3.190382-01 

3.955332-01 

9.249202-01 

9 .25X412-02 

;a 

3.457932-04 

3.847522-03 

3.198842-01 

•0.247502-03 

-0. 155002-02 

-0.200595-92 

19 

-0.315902-04 

-0. 551532-03 

0.2002X2-01 

3.450432-01 

-0.474022-02 

9.759445-95 

20 

-0.325932-04 

0.521352-03 

9.181812-01 

9.754742-01 

-0.4722S2-02 

-0.550245-95 

21 

-9.383702-94 

-0.491582-03 

9.143092-01 

9.932432-01 

-0.344942-02 

9.702885-95 

22 

-9.451442-04 

-0.442402-03 

9.138532-01 

9.953332-01 

-0.144942-03 

-0.788246-05 

22 

-9.448382-04 

-0.434002-03 

0. U48S2-01 

9.745142-01 

-0.247942-02 

0*122415-02 

2* 

-0. 441472-04 

-0. 324332-03 

0.884572-02 

0.12X37 

-0.444772-02 

-9*747035-03 

25 

-0. 429022-04 

0.214112-03 

0.544212-02 

0.15425 

0.272052-03 

9 *100455-02 

25 

-9.343052-04 

-4.105302-03 

9.273842-02 

9.11491 

0.812572-02 

-0.415132-03 

ZT 

-0.244222-04 

3.00000 

0.00000 

0. S79O82-01 

0.442492-02 

9.101395-92 

22 

3.140852-04 

3.00000 

0.00000 

-0.409572-01 

9.148832-02 

9.339915-93 

2* 

3.197X22-00 

3.143542-03 

9 .222182-02 

0.15770 

0 .250992-01 

-9.113372-92 

10 

C.JIMS&04 

3.332*45-03 

9.540052-02 

9.1987 8 

0.108032-03 

9*935935-93 
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51 

*.395452-44 

9.501402-43 

0.934712-02 

9.14157 

-4. 232352-91 


32 

3 .554722-45 

3.664742-43 

9.117892-91 

*0.23X722-01 

-4.567452-95 

9.12Z7X3H 

53 

0 .514562-45 

9.713792-43 

9.133162-91 

3 .999*72-01 

0*153542-91 

-«.UBX£H 

5* 

9.475442-44 

9.763322-03 

9.157442-41 

9.13592 

-4*969522-95 

I.IHMEH 

55 

9.423535*44 

3.404952-03 

0. 143162-41 

9.971892-01 

-0.144342-41 

•*.420SH 

55 

3.43371E-44 

0.451672-03 

9 .299392-91 

0.441442-93 

•4.113992-92 

1.MU4EH 

57 

9.40949 

•0.524362-03 

9*289992-91 

3*463462-91 

0*364262-99 

-3.941SKH 

51 

9.49949 

•0.501662-03 

9*144242-91 

9.737172-91 

•4*184442-99 

9.337233*1 

39 

9.94499 

•4.473272-93 

9*163722-91 

9*946492-91 

0*117982-99 

8.234413*1 

59 

9.49949 

-4*442622-03 

0*144912-91 

9*914972-91 

-4*155422-99 

3.444733-1 

51 

9.09499 

•4.410242-03 

9 .117522-91 

0*736442-91 

0 *425942-99 

-3.115373-1 

52 

9.00949 

•4.304112-03 

9*93X662-92 

9*12476 

-4*149762-49 


55 

9.00000 

•0*205462-03 

9*573292-92 

0*13433 

0*119972-99 

•«jsbc<] 

54 

3.00000 

•4*103252-03 

9*194492-92 

9*11446 

-4*304442-99 


55 

3.00090 

0*00000 

9. 94499 

0*576462-91 

0.643762-99 

9.1*3793-2 

54 

3.00094 

0*00000 

0*40494 

*9*423442-91 

•4*964472-49 

i.UUS*l 

47 

3.00099 

9. 171502-03 

9*532992-93 

9*13621 

0*632772-49 

a.iQPit-1 

54 

3.00004 

0*341522*03 

9.574452-92 

9.19429 

•4*446166-99 


44 

3.90004 

9*512692-03 

9*184382-91 

0*14474 

0*620232-49 

3.12272*1 

55 

9.00004 

0*647542-03 

9 *117752-91 

-9.234442-91 

-4*432572-49 

-3.493S2E-1 

51 

3.00004 

0*733732-03 

9*123132-91 

9*967222-91 

0 *534196-99 


52 

9.00000 

0.740392-93 

0*154542-91 

0*13451 

-4.358262-99 

I4HUK 

55 

3.00000 

9*427362-93 

0*192542-91 

0.945352-91 

0*443712-99 

8 

54 

3.00009 

0*475742-93 

0*289972-91 

0.399692-93 

•4*494522-49 

U(S4E>1 

55 

3 .3X5902*44 

*0*551532-93 

0*2902X2-91 

0*459632-91 

0.476922-42 


54 

9.326932-04 

-9.521352-93 

0*141412-91 

0*754762-91 

0*672252-42 

UIIMM 

57 

3.344702-04 

0.491542-03 

9*163092-91 

3*932432-01 

9.366966-42 

-3.702843-3 

55 

3.451652-05 

*9.462692-03 

9*134532-91 

9.934442-01 

9.144946-43 

9.7SS24E-3 

54 

3 . 564342-04 

*9.434002-43 

3.116452-01 

9.745142-01 

9.267946-42 

-9.122413-3: 

o* 

3.551472-04 

-9.324332-03 

3.446672-02 
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